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Abstract 


A  model  developed  to  analyze  the  global  emergent  dynamics  of  situational 
awareness  (SA)  dissemination  witiun  a  digitized  brigade  undergoing  hostile 
information  operations  (lO)  stress  events  is  presented.  Networked  battlefield 
platforms  are  modeled  as  interconnected  tactical  internet  nodes  through  the  use 
of  cellular  automata.  The  cellular  automata  form  a  discrete  spatially  extended 
d5mamical  system  consisting  of  a  parallel  networked  lattice  of  computational 
cells  in  two  dimensions.  The  global  impact  of  localized  lO  stress  events  on  SA 
dissemination  can  then  be  analyzed  as  a  function  of  time  by  mapping  the  SA 
dissemination  architecture  onto  the  cellular  automata  lattice.  The  lattice  cells  can 
represent  tactical  internet  client  or  server  platforms.  This  model  demonstrates 
the  fluctuating  patterns  of  situational  "self-awareness"  (i.e.,  a  friendly  platform's 
situational  awareness  of  all  other  friendly  platforms  within  its  area  of  operations) 
which  emerge  across  the  digitized  brigade  on  a  global  scale  and  also  illustrates 
the  emergent  d5mamics  of  SA  dissemination  within  a  "movement-to-contact" 
brigade  structure. 
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Executive  Summary 


In  this  report,  a  model  for  analyzing  the  emergent  global  impact  of  hostile 
information  operations  (lO)  stress  directed  at  the  networked  platforms  making 
up  a  digitized  brigade-level  situational  awareness  (SA)  architecture  is  presented. 
Since  the  model  specifically  addresses  a  platform-  to  network-level  mapping  of 
structure  and  behavior,  the  complexity  of  message  traffic  patterns  (which  arises 
from  intra-platform  processes  embedded  within  the  platform-to-network-level 
mapping)  cannot  be  properly  quantified  and  analyzed.  Thus,  the  current  model 
focuses  on  the  dynamic  robustness  of  file  brigade  SA  network  topology  (i.e.,  the 
potential  ability  of  the  network  structure  to  perpetuate  SA  dissemination  under 
disruptive  stress)  within  a  hostile  lO  enviroiunent. 

A  brigade-level  tactical  conunurucatioiis  network,  as  implemented  tiirough  the 
U.S.  Army's  Force  XXI  Battle  Coiiunand  Brigade-and-Below  (FBCB2),  consists  of 
two  distinct  subarchitectures  for  the  parallel  dissemination  of  both  situational 
awareness  (SA)  and  command  and  control  (C2)  message  traffic  through  the  use 
of  enhanced  position  location  reporting  system  (EPLRS)  data  radios.  The  SA 
subarchitecture  within  a  digitized  brigade  uses  host  platforms  functioning  as 
position  servers  to  collect  and  distribute  the  data  vertically  and  horizontally 
across  a  predetermined  hierarchy  of  radio  nets,  which  include  local  single 
channel  ground  airborne  radio  system  (SINCGARS)  based  nets  as  weU  as  &e 
EPLRS-based  battaUon-level  carrier  sense  multiple  access  (CSMA)  nets  and 
brigade-wide  multisotirce  group  (MSG)  net.  This  SA  dissemination  architecture 
provides  a  framework  within  which  to  construct  a  cellular-automata-based  SA 
network  topology  model. 

In  general,  a  cellular  automata  (CA)  consists  of  a  possibly  infimte,  n-dimensional 
regular  lattice  of  cells.  Each  cell  can  be  in  a  measurable  state  chosen  from  a  finite 
alphabet  (i.e.,  a  set  of  possible  cell  states).  The  states  of  aU  cells  in  the  lattice  are 
updated  simultaneously  in  discrete  time  steps  through  the  use  of  a  transition 
function,  which  takes  as  its  input  the  current  state  of  the  center  cell  and  some 
finite  collection  of  nearby  cells  tiiat  lie  within  a  finite  distance,  collectively  known 
as  a  neighborhood.  CeU  neighborhoods  may  be  either  fixed  or  variable  over  time 
(but  fixed  witihin  a  time  step);  these  neighborhoods  are  usually  defined  to  be 
local  to  a  center  reference  cell,  but  can  extend  beyond  this  conventional  limit  if 
necessary.  This  last  condition  is  necessary  to  model  the  SA  dissemination 
network  topology  within  a  brigade,  where  neighborhoods  must  span  across  all 
lattice  scales  (i.e.,  local  to  global). 


The  first  step  in  building  the  SA  network  is  to  spatially  coiafigure  a  brigade 
structure  relative  to  a  CA  lattice;  such  an  example  network  is  a  generic 
''movement  to  contact"  brigade  structure.  Within  this  generic  structure,  each  of 
the  612  blue  squares  represents  a  networked  platform  node  which  occupies  a 
lattice  ceU  (in  the  current  version  of  the  model,  nodes  are  ground  vehicles  which 
remain  fixed  relative  to  one  another),  where  a  node  is  designated  as  either  a 
SINCG ARS  net  client,  a  self  server,  a  CSMA  server,  or  an  MSG  server. 

The  next  step  in  model  development  is  to  define  the  spatial  scope  of  situational 
awareness  operationally  required  of  a  platform  in  order  to  formulate  one  or  more 
CA  neighborhoods  relative  to  the  lattice  cell  which  tihe  platform  occupies. 
Applying  an  echelon-dependent  "area  of  operations"  (AO)  rectangle  to  each 
platform  within  the  CA  brigade  model  defines  an  AO  neighborhood  relative  to  the 
overall  CA  lattice,  where  the  AO  neighborhood  is  a  sublattice  stirrounding  a 
center  node  cell  which  defines  the  region  of  friendly  neighboring  node  cells 
which  the  center  node  must  be  situationally  aware  of. 

Once  the  spatial  scope  of  situational  awareness  has  been  delineated  and 
implemented  for  each  node  ceU  within  the  brigade  (thus  creating  an  AO-based 
network  topology),  a  quantitative  measure  of  the  degree  to  which  a  node  is 
situationally  aware  of  other  friendly  neighbor  nodes  can  now  be  defined. 
Within  the  context  of  the  CA  model,  the  situational  awareness  of  a  TI  node  ceU  is 
a  collective  measure  of  the  state  of  the  communication  channels  which  exist 
between  the  reference  node  and  all  neighboring  friendly  nodes  within  the 
reference  node's  AO  neighborhood.  The  state  of  the  local  network  topology 
relative  to  the  z*  node  can  be  represented  by  a  metric  called  the  SA  connectivity 
Ci(tn)  that  indicates  the  fraction  of  neighbor  nodes  within  the  z*  node's  AO 
neighborhood  that  can  transmit  SA  messages  to  the  z*^  node  at  time  step  t„. 

There  are  five  data  processing  filters  that  have  been  coded  into  the  CA  model  for 
analyzing  time  series  data. 

(1)  The  first  data  filter  calculates  the  normalized  perturbative  stress 

that  is  introduced  into  the  entire  brigade  system  via  localized  lO  stress. 

(2)  The  second  data  filter  calculates  globally  averaged  SA  connectivity  <C(f„)> 
as  a  fimction  of  scenario  time,  where  an  average  value  is  calculated  at  each 
point  in  discrete  simulation  time  across  all  node  platforms  within  the  CA 
brigade. 

(3)  The  tiiird  data  filter  calculates  the  five  fitness  metrics  Pvery  high{tn),  Phigh{t„), 
Pmcdiumitn),  Pi„u,(t„),  and  Zou,  (fn),  where  Pfitness(fn)  is  the  fraction  of  network 
nodes  occupying  a  particular  "fitness"  state  at  time  step  t„.  The  five  fitness 
states  are  defined  by  the  inteiwals  C(fn)  <  0.1  (veiy  low),  0.1  <  Ci(t„)  <  0.4 
(low),  0.4  <  C,(f„)  <  0.6  (medium),  0.6  <  G(f„)  <  0.9  (high),  and  G(f„)  >  0.9  (very 
high). 
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(4)  The  fourth  data  filter  calculates  the  state  entropy  H{tn),  a  quantitative 
measure  of  the  disorder/volatility  within  the  network  at  time  step  t„,  and 
the  temporal  entropy  H{t„,  a  quantitative  measure  of  the  temporal 
disorder  between  sequential  network  states  measured  at  f„-i  and  t„. 

(5)  The  last  data  filter  presents  fitness  metric  data  in  the  form  of  a  return  map, 
which  is  a  2-D  scatter  plot  of  Pfitness{tn+i)  vs.  Pfitness{t„)  for  each  fitness  state. 

Taken  together,  these  filters  can  be  used  to  analyze  the  emergent  dynamics  of  an 
SA  network  exposed  to  lO  stress. 

The  final  step  in  developing  the  CA  model  is  the  addition  of  some  form  of  lO 
stress  which  directly  interacts  with  a  node  cell  in  order  to  perturb  its 
communication  capabilities.  This  stress  can  arise  either  from  a  stochastic 
parameter  embedded  within  a  nodal  state  vector  or  an  external  neighboring  cell 
(representing  a  localized  lO  threat  soiuce)  that  generates  a  local  "stress  field." 
Examples  of  both  types  of  lO  stress  were  added  to  the  CA  model,  and 
subsequent  simulation  results  were  analyzed. 

The  first  type  of  lO  stress  was  modeled  as  a  stochastic  process  which  acts  to 
perturb  the  functional  commtmication  states  of  each  node  witiun  the  SA 
network.  The  dynamics  of  this  stochastic  process  are  described  by  a  stationary 
Markov  chain  which  is  appUed  uniformly  and  independently  to  each  of  the 
nodes  within  the  digitized  brigade.  The  stochastic  process  then  acts  to  transition 
each  network  node  between  functional,  transient,  or  permanent  dysfunction,  and 
functional  but  "spoofed"  states.  This  type  of  stress  is  similar  to  the  mean  field 
approximation  from  statistical  physics,  where  a  particular  element  of  a  multi¬ 
element  system  is  subjected  to  tiie  average  influence  of  a  global  field  interacting 
with  the  rest  of  the  system.  In  this  sense,  the  "mean  field"  lO  stress  is  at  best  a 
first-order  approximation  to  a  real  multi-tiueat  lO  environment,  serving  to 
demonstrate  the  emergent  behavior  of  the  brigade's  SA  network  topology  imder 
a  coherent  hypothetical  rather  than  realistic  stress. 

Simulation  results  using  the  mean  field  lO  stress  model  show  that  message 
flooding  (resultant  from  Global  Positioning  System  [GPS]  spoofing)  is  the 
dominant  cause  of  SA  connectivity  degradation,  since  extreme  levels  of 
communication  channel  latency  (resulting  from  server  message  queue 
saturation)  are  also  assumed  to  cut  intemodal  commrmication  channels  tmtil 
latency  levels  drop  below  a  threshold  level.  However,  the  mean  field  model  has 
shortcomings  in  two  different  areas.  First,  since  the  C2  network  structure  within 
the  brigade  has  deliberately  been  decoupled  from  the  SA  network  structure,  the 
model  can  significantly  imderpredict  the  impact  of  GPS  spoofing  on  nodal  SA 
connectivity.  Second,  the  model  fails  to  realistically  consider  correlations 
between  spatially  dispersed  stress  sources,  as  well  as  the  resultant  correlations 
between  SA  connectivity  levels  associated  with  impacted  nodes. 
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A  more  realistic  form  of  lO  stress  arises  from  perturbations  to  nodal 
communication  due  to  localized  external  fields  that  are  generated  by 
electromagnetic  (EM)  sources.  In  this  second  stress  model,  the  CA  d5nnamics  are 
added  by  introducing  spatially  dispersed  source  cells  representing  generic  radio 
frequency  (RF)  "jammer  bombs"  that  act  to  locally  jam  platform-mounted  radios 
passing  by.  Thus,  a  node  is  jammed  if  it  passes  within  a  predefined  jamming 
neighborhood  of  a  jammer  bomb;  alternately,  a  node  is  jammed  if  and  only  if 
there  is  at  least  one  jammer  bomb  within  an  identical  jamming  neighborhood 
relative  to  the  target  node  (where  the  target  node  is  at  the  center  of  the  translated 
neighborhood).  Since  the  localized  jammer  bomb/node  interactions  act  to 
perturb  nodal  reception  capability  states,  this  in  turn  perturbs  nodal  SA 
connectivity  levels  throughout  the  CA  brigade. 

A  jammer  bomb  scenario  is  defined  with  two  phases:  (1)  the  digitized  brigade 
moves  through  a  field  of  jammer  bombs  previously  planted  by  hostile  forces,  and 
(2)  hostile  unmanned  aerial  vehicles  flying  over  the  advancing  brigade  deliver  a 
barrage  of  additional  jammer  bombs.  Twenty  simulation  trials  of  the  two-phase 
jammer  bomb  scenario  are  nm  using  the  CA  model,  where  each  trial  is  seeded 
wifli  a  different  random  ntimber  (needed  to  stochastically  vary  the  distribution 
of  jammer  bombs  upon  initial  insertion  into  a  simtdation).  In  addition,  four 
variants  to  the  original  jammer  bomb  scenario  are  also  simulated,  which  include 
(1)  reduced  AO  neighborhoods,  (2)  reduced  jammer  bomb  lifetimes,  (3)  a 
combination  of  reduced  jammer  bomb  neighborhoods  and  increased  average 
bomb  density,  and  (4)  jammer  bomb  susceptibility  gradients. 

Simulation  results  from  all  of  the  jammer  bomb  scenarios  indicate  the  following 
set  of  shared  characteristics. 

•  The  dynamic  fraction  of  jammed  node  within  the  brigade  increases  to  a 
maximum  value,  temporarily  remains  in  a  state  of  quasi-equiHbrium,  and 
then  steadily  decreases  back  down  to  a  level  of  0. 

•  The  dynamic  SA  connectivity  averaged  across  aU  nodes  within  the  brigade 
decreases  down  to  a  minimum  value,  again  temporarily  remains  in  a  state 
of  quasi-equiHbrium,  and  then  steadily  increases  back  up  to  a  level  of  1. 

•  The  SA  fitness  Fvery  iow{tn)  time  series  roughly  track  the  scenario-respective 
</z5'’’«*(f„)>  time  series  in  terms  of  metric  ampHtude. 

•  The  SA  fitness  Pvery  high{tn)  time  series  roughly  track  the  scenario-respective 
<C(f„)>  time  series  in  terms  of  relative  time  series  profile  (with  differences  in 
metric  ampHtude). 
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•  The  SA  fitness  Phigh{tn)  time  series  profile  is  a  mirror-image  reflection  of  the 
respective  SA  fitness  Pveryhigh{tn)  time  series  profile. 

•  The  dynamic  state  entropy  H{t„)  and  temporal  entropy  H{t„,  tn-i) 
demonstrate  separation  between  global  brigade,  CSMA  net,  and  local  net 
node  populations. 

However,  the  jammer  bomb  scenario  variants  are  distinguishable  by  their 
distinctive  emergent  behaviors  as  displayed  within  their  respective  Pveiy  kigh{tn) 
time  series.  Although  variations  in  AO  neighborhoods  and  jammer  bomb 
lifetimes  lead  to  differences  in  collective  structure  as  displayed  within  the  return 
maps,  the  most  noticeable  factor  observed  within  jammer  bomb  simulations 
appears  to  be  the  size  of  the  jammer  bomb  neighborhood  and  its  impact  on  Pvery 
high{tn)-  These  latter  results  suggest  fliat  the  greater  the  jamming  range  of  a  stress 
source  and  tire  more  ubiquitous  this  type  of  source,  the  more  coordinated  the 
overall  network  response  (i.e.,  the  less  resistant  is  the  network  to  associated 
perturbations  to  SA  connectivity). 

The  spatio-temporal  dynamics  of  other  t3q)es  of  dispersed  lO  stress  sources  that 
modify  nodal  transmit/ receive  functionality  are  likely  to  create  similar  modes  of 
emergent  global  behavior  in  situational  awareness  levels  which  cannot  be 
anticipated  by  studying  isolated  platforms  or  even  local  area  networks  (LANs). 
Given  that  the  threat/ target  interaction  dynamics  of  an  lO  stress  source  can  be 
modeled,  it  is  a  fairly  straightforward  process  to  program  the  stress  source  into 
the  CA  code  within  the  context  of  a  scenario.  There  is  also  no  limitation  on  the 
number  of  different  types  of  lO  stress  sources  which  can  be  programmed  into  a 
single  scenario  given  that  the  discrete  spatio-temporal  interaction  d5namics  are 
known. 
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1.  Introduction 


1.1  Purpose 

The  cellular-automata-based  network  model  presented  in  this  report  was 
developed  to  analyze  the  emergent  global  dynamics  of  situational  awareness 
(SA)  dissemination  within  a  digitized  brigade  exposed  to  several  different  types 
of  hostile  information  operations  (lO)  stress  (U.S.  Department  of  the  Army 
1996b). 

1.2  Background 
1.2.1  General 

A  complex  system  can  be  defined  as  a  group  of  many  interacting  parts  functioning 
as  a  unified  whole.  It  is  distinguishable  from  its  environment  by  recognizable 
boundaries,  where  the  global  behavior  of  the  whole  is  generally  not  equal  to  the 
sum  of  the  parts.  Even  if  the  individual  behavior  of  an  isolated  part  is  linear,  the 
collective  behavior  that  emerges  globally  through  networked  interactions  and 
interdependencies  among  the  parts  is  often  nonlinear.  This  collective  behavior 
arises  from  the  detailed  structure,  behavior,  and  interpart  relationships  as  they 
are  defined  on  a  finer  scale,  aU  of  which  dynamically  evolve  within  the  context  of 
an  irreducible  whole. 

To  analyze  the  U.S.  Army's  evolving  tactical  internet  (TI)  as  a  complex  system 
operating  within  an  equally  complex  lO  environment,  the  following  three 
distinct  scales  define  the  structure,  behavior,  and  relationships:  (1)  the  microscale, 
which  addresses  interactions  and  interdependencies  between  information  system 
(IS)  components  (e.g.,  1553b  data  bus,  combat  net  radios,  personal  computers) 
within  a  platform,  (2)  the  mesoscale,  which  addresses  local  interactions/ 
interdependencies  between  networked  battlefield  platforms,  and  (3)  the 
macroscale,  which  addresses  the  holistic  global  behavior  of  collective  "system-of- 
systems"  structures  such  as  the  TI  (zum  Bnmnen  et  al.  2000).  Usually,  lO 
threat/ target  system  interactions  occur  at  the  micro-  and/or  mesoscale  levels, 
affecting  information  flows  between  IS  components  within  a  platform,  and/  or 
flows  between  locally  networked  platforms,  respectively.  In  the  current  model, 
the  emergent  global  or  macroscale  impact  of  lO  stress,  which  is  directed  at  the 
mesoscale  networked  platforms  making  up  a  digitized  brigade-level  situational 
awareness  (SA)  architecture  is  analyzed.  Since  the  model  specifically  addresses  a 
meso-  to  macroscale  mapping  of  structure  and  behavior,  the  complexity  of 
message  traffic  patterns  (which  arises  from  micro-  to  mesoscale  processes 
embedded  within  the  mesoscale/ macroscale  mapping)  cannot  be  properly 
quantified  and  analyzed.  Thus,  this  model  focuses  on  the  d)mamic  robustness  of 
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the  brigade  SA  network  topology  (i.e.,  the  potential  ability  of  the  network 
structure  to  perpetuate  SA  dissemination  under  disruptive  stress)  within  a 
hostile  lO  environment. 

A  brigade-level  tactical  communications  network,  as  implemented  through  the 
U.S.  Army's  Force  XXI  Battle  Command  Brigade-and-Below  (FBCB2),  consists  of 
two  distinct  subarchitectures  for  the  parallel  dissemination  of  both  situational 
awareness  (SA)  and  command  and  control  (C2)  message  traffic  using  enhanced 
position  location  reporting  system  (EPLRS)  data  radios.  The  EPLRS-based  radio 
nets  provide  contention-free  routing  for  both  SA  and  C2  traffic  down  to  local 
single  channel  groimd  airborne  radio  system  (SINCGARS)  nets  which  are  also 
known  as  stub  nets  (since  tihey  function  as  jimction  stubs  for  the  SA  and  C2 
subarchitectures)  (U.S.  Department  of  the  Army  1999). 

The  SA  subarchitecture  within  a  digitized  brigade  uses  host  platforms 
fimctioning  as  position  servers  to  collect  and  distribute  the  data  vertically  and 
horizontally  across  a  predetermined  hierarchy  of  radio  nets.  These  radio  nets 
include  local  SlNCGARS-based  nets,  as  well  as  the  EPLRS-based  battalion-level 
carrier  sense  multiple  access  (CSMA)  nets  and  the  brigade-wide  multisource 
group  (MSG)  net  (see  Figure  1)  (U.S.  Department  of  the  Army  1999;  U.S.  Army 
Communication  Electronics  Command  1999).  Any  digitized  platform  equipped 
with  an  EPLRS  radio  unit  can  function  as  a  server.  These  servers  include  (a) 
CSMA  servers  that  route  SA  traffic  upward  from  their  local  SINCGARS  net  to 
their  CSMA  net  and  downward  from  both  CSMA  and  brigade- wide  MSG  nets  to 
local  nets,  (b)  MSG  servers  that  route  traffic  from  a  CSMA  net  upward  to  the 
MSG  net  and  (c)  self  servers  fliat  broadcast  self-generated  SA  messages  directly 
onto  their  CSMA  community  net  and  are  thus  not  members  of  a  local  SINCGARS 
net.  Whereas  CSMA  server  selection  is  dynamically  updated  amongst  aU 
capable  platforms  within  a  local  SINCGARS  community  (although  each 
community  must  have  at  least  a  primary  and  alternate  server),  primary  and 
alternate  MSG  servers  are  predetermined.  In  addition,  aU  platforms  equipped 
with  EPLRS  radio  units  receive  aU  MSG  net  broadcasts  within  the  brigade,  as 
weU  as  aU  broadcasts  within  the  platform's  respective  CSMA  net  community 
(note  that,  unless  they  are  CSMA  or  self  servers,  these  platforms  stUl  require  a 
local  CSMA  server  to  broadcast  self-generated  SA  messages  outside  of  their  local 
SINCGARS  net).  This  SA  dissemination  architecture  provides  a  framework 
within  which  to  construct  a  ceUular-automata-based  SA  network  topology 
model. 

1.2.2  Threat 

The  threats  addressed  within  this  report  are  information  operation  events  which 
result  in  either  denial  of  service  or  message  flooding  within  a  hierarchical 
digitized  brigade  network. 
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Figure  1.  Situational  awareness  dissemination  architecture  within  a  digitized  brigade. 


1.3  Scope 

The  applicability  of  the  SA  dissemination  model  presented  in  this  report  is 
constrained  by  die  following  assumptions. 

•  Networked  platform  nodes  within  a  digitized  brigade  receive  SA  update 
messages  from  other  platform  nodes  only  (i.e.,  network  control  stations 
[NCS]  are  not  modeled). 

•  lO  stress  will  only  impact  a  node's  ability  to  transmit  and/or  receive  SA 
messages  from  other  neighboring  nodes  and  not  affect  internal  information 
flows  within  a  platform  node. 

•  The  spatial  configuration  of  platform  nodes  relative  to  one  another  remains 
fixed  throughout  a  simulation  (i.e.,  platforms  do  not  move  relative  to  each 
other). 


2.  The  Situational  Awareness  Network  Model 


2.1  Cellular  Automata 

In  general,  a  cellular  automata  (CA)  (von  Neumann  1966)  consists  of  a  possibly 
infinite,  n-dimensional  regular  lattice  of  cells.  Each  cell  can  be  in  a  measurable 
state  chosen  from  a  finite  alphabet  (i.e.,  a  set  of  possible  cell  states).  The  states  of 
aU  cells  in  the  lattice  are  updated  simultaneously  in  discrete  time  steps  using  a 
transition  function,  which  takes  as  its  input  the  current  state  of  the  center  cell  and 
some  finite  collection  of  nearby  cells  that  lie  within  a  finite  distance,  collectively 
known  as  a  neighborhood.  Figure  2  displays  two  well-known  CA  cell 
neighborhoods  in  two  dimensions,  where  C  =  center  cell  and  N  =  neighbor  cell. 
Cell  neighborhoods  may  be  either  fixed  or  variable  over  time  (but  fixed  within  a 
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Figure  2.  CA  cell  neighborhoods  in  two  dimensions:  (a)  von  Neumaim  neighborhood 
and  (b)  Moore  neighborhood. 

time  step);  these  neighborhoods  are  usually  defined  to  be  local  to  a  center 
reference  ceU,  but  can  extend  beyond  this  conventional  hmit  if  necessary.  This 
last  condition  is  necessary  to  model  the  SA  dissemination  network  topology 
within  a  brigade,  where  neighborhoods  must  span  across  all  lattice  scales  (i.e., 
local  to  global). 

2.2  The  Cellular  Simulation  System 

The  SA  network  topology  model  is  implemented  using  die  Cellular  simulation 
system,  an  easy-to-use  toolkit  for  the  modeling  of  physical  systems  with  cellular 
automata.  The  Cellular  system  consists  of  several  modules: 

•  A  progranuning  language,  Cellang,  and  associated  compiler,  cellc. 

•  An  abstract  virtual  cellular  automata  machine  (avcam)  for  Cellang  code 
execution. 

•  Cell  view,  a  program  for  viewing  simulation  results  in  either  two- 
dimensional  (2-D)  sections  or  in  three  dimensions  with  variable 
perspective. 

Compiled  Cellang  programs  can  be  nm  with  input  provided  at  any  specified  time 
during  the  execution.  The  results  of  an  execution  can  either  be  viewed 
graphically,  as  an  output  stream  of  ceU  locations  and  values,  or  passed  through  a 
custom  filter  before  being  reported.  The  output  stream,  or  a  similar  stream 
produced  by  a  custom  filter,  can  edso  be  directed  into  cellview  for  viewing,  or 
may  be  passed  through  other  programs  that  compile  statistics,  massage  the  data, 
or  merely  act  as  a  valve  to  control  the  flow  of  data  from  the  cellular  automata 
program  to  the  viewer.  For  more  detailed  information  on  the  Cellular  toolkit,  see 
Eckart  (1992a)  and  Eckart  (1992b). 
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2.3  Building  the  Cellular  Automata  Network 

2.3.1  The  Network  Cell 

The  first  step  in  building  the  CA-based  SA  network  is  to  define  the  attributes  of  a 
node  platform  within  the  network.  Within  the  context  of  a  CA  lattice,  a  node  is 
defined  as  a  state  vector  that  represents  a  networked  platform  occupying  a  cell 
within  a  2-D  CA  lattice.  A  nodal  state  vector  maintains  the  following 
information  fields: 

•  node  cell  location  within  the  lattice  (using  Cartesian  coordinates), 

•  nodal  communication  capability  state  (transmission  and  reception 
functionality), 

•  platform  military  echelon, 

•  membership  in  specific  local  SINCGARS  and/or  CSMA  radio  nets, 

•  server  status  (i.e.,  node  is  either  a  SINCGARS  net  client,  a  self  server,  a 
CSMA  server,  or  an  MSG  server), 

•  nodal  SA  connectivity  (see  section  2.4.2),  and 

•  nodal  SA  fitness  (see  section  2.6.3). 

When  appropriate,  the  nodal  state  vector  can  also  keep  running  totals  of  SA 
messages  received  from  the  various  radio  nets  to  which  a  node  belongs.  In  the 
current  version  of  the  network  model,  all  nodes  are  assumed  to  be  groxmd 
vehicles  which  remain  spatially  fixed  relative  to  one  another. 

2.3.2  The  Cellular  Automata  Brigade 

The  next  step  in  building  the  network  is  to  spatially  configure  a  brigade  structure 
relative  to  the  CA  lattice;  an  example  network  is  the  generic  movement  to  contact 
brigade  structure  (U.S.  Department  of  the  Army  1996a)  depicted  in  Figure  3. 
Each  of  the  612  blue  squares  in  the  figure  represents  a  networked  platform  node 
that  occupies  a  unique  (x,  y)  position  within  the  CA  lattice.  Each  of  these 
positions  is  assigned  a  nodal  state  vector  that  is  associated  with  the  platform 
node  occup)ting  the  corresponding  lattice  cell.  Server  nodes  are  mutually 
coordinated  in  that  an  alternate  CSMA  or  MSG  server  will  assume  tiie  server 
identity  if  the  primary  server  becomes  dysfunctional.  This  identity  will  be  used 
only  to  construct  an  SA  network  topology  witiun  the  context  of  a  CA  lattice. 
Also,  the  brigade  layout  depicted  in  Figure  3  serves  only  as  an  illustrative 
example  that  approximates  a  movement  to  contact  structure;  a  real  brigade 
layout  (which  can  also  be  programmed  for  simulation  by  the  CA  model)  would 
likely  be  less  spatially  ordered  than  the  neat  "building  block"  configuration 
shown  in  the  figure  (Wuerz  2000). 
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Figure  3.  Cellular  automata  configuration  of  a  movement  to  contact  brigade  structure. 

2.4  Measuring  Situational  Awareness 
2.4.1  Area  of  Operations  Neighborhoods 

The  next  step  in  model  development  is  defining  the  spatial  scope  of  situational 
awareness  that  is  operationally  required  of  a  platform  in  order  to  formulate  one 
or  more  CA  neighborhoods  relative  to  the  lattice  cell  the  platform  occupies.  In 
the  FBCB2  Operational  Requirements  Document  (ORD),  the  key  performance 
parameter  (KPP)  relevant  to  SA  functionality  states  that  the  system  shall 
"...display  friendly  positions  horizontally  within  a  rmit,  two  echelons  up  and 
down,  and  adjacent  unit  location  one  echelon  down"  (U.S.  Army  Training  and 
Doctrine  Command  1998).  To  allow  platforms  equipped  only  with  SINCGARS 
to  meet  this  KPP,  CSMA  servers  filter  incoming  SA  traffic  based  on  an  area  of 
operations  (AO)  rectangle  drawn  around  a  client  platform  (with  the  client  at  the 
rectangle  center),  where  rectangle  area  increases  with  higher  client  platform 
echelon.  In  this  approach,  a  CSMA  server  will  only  deliver  messages  to  a  client 
which  originate  from  friendly  neighboring  platforms  within  the  client's  AO 
rectangle;  this  is  done  to  prevent  message  saturation  within  the  local  SINCGARS 
net  (U.S.  Army  Communication  Electronics  Command  1999). 

Applying  an  echelon-scaled  AO  rectangle  (which,  in  this  application,  is 
constrained  to  an  equal-sided  square)  to  each  platform  within  the  CA  brigade 
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model  defines  an  AO  neighborhood  relative  to  the  overall  CA  lattice.  Tlie  AO 
neighborhood  is  a  sublattice  surrounding  a  center  node  cell  that  defines  the 
region  of  friendly  neighboring  node  cells  which  the  center  node  must  be 
situationally  aware  of.  Figure  4  displays  the  four  AO  neighborhoods  measuring 
20  km  X  20  km,  15  km  x  15  km,  10  km  x  10  km,  and  5  km  x  5  km,  as  they  are 
defined  for  brigade,  battalion,  company,  and  platoon  echelon  node  platforms, 
respectively  (U.S.  Army  Communication  Electronics  Command  1999);  for  this 
application,  the  CA  lattice  has  been  scaled  to  1  lattice  cell  length  =  125  m.  Note 
that  neighborhood  scaling  is  adjustable  within  the  CA  model  to  reflect  variations 
in  mission,  enemy,  troops,  terrain,  and  time  av'ailable  (METT-T)  conditions. 
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Figure  4.  Area  of  operations  neighborhood  scaling  within  the  CA  brigade. 

2.4.2  Situational  Awareness  Connectivity 

Once  the  spatial  scope  of  situational  awareness  has  been  delineated  and 
implemented  for  each  node  cell  within  the  brigade  (thus  creating  an  AO-based 
network  topology),  a  quantitative  measure  of  the  degree  to  which  a  node  is 
situationally  aware  of  other  friendly  neighbor  nodes  can  now  be  defined.  Within 
the  context  of  the  CA  model,  the  situational  awareness  of  a  T1  node  cell  is  a 
collective  measure  of  the  state  of  the  communication  channels  which  exist 
between  the  reference  node  and  all  neighboring  friendly  nodes  within  the 
reference  node's  AO  neighborhood.  This  is  illustrated  in  Figure  5,  which  depicts 
a  sublattice  structure  of  a  center  receiving  node  (represented  by  an  "R")  and  the 
transmitting  neighbor  nodes  (represented  by  a  "T")  that  lie  within  the  receiving 
node's  AO  neighborhood  (the  "R/T"  notation  represents  server  nodes  which 
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Figure  5.  Local  network  topology  of  a  center  receiving  node's  AO  neighborhood. 


must  be  functionally  capable  of  both  receiving  and  transmitting).  If  the  binary 
states  and  s/{t„)  are  defined  as  the  i*  node's  ability  to  transmit  (l)/not 

transmit  (0)  and  receive  (l)/not  receive  (0)  SA  messages  at  time  tn  (i.e.  the  n* 
discrete  simulation  time  step),  respectively,  then  the  state  of  the  local  network 
topology  relative  to  the  node  can  be  represented  by  a  metric  called  the  SA 
connectivity,  as  in 


^  i  j 


where  N,  =  number  of  friendly  neighbor  nodes  within  the  i*  node's  AO 
neighborhood  and  ]ij{tn)  =  binary  communication  channel  continuity  from  the 
to  node.  This  channel  continuity  metric  is  equal  to  1  only  if  (a)  they*  neighbor 
node  lies  within  the  i*  node's  AO  neighborhood,  and  (b)  the  channel  connecting 
nodes  i  and  j  are  intact  at  time  step  f„.  Thus,  C(t„)  indicates  the  fraction  of 
neighbor  nodes  within  the  i*  node's  AO  neighborhood  that  can  transmit  SA 
messages  at  t„.  The  representation  of  SA  connectivity  expressed  in  equation  (1)  is 
analogous  to  the  local  energy  E,  within  a  two-spin,  variable-range  king 
model  from  statistical  physics,  which  has  also  been  applied  to  the  modeling  and 
simulation  of  diverse  cooperative/interdependent  phenomena  such  as  protein 
folding  (Bar-Yam  1997)  and  political  decision-related  structures  (Galam  1997). 


It  should  be  noted  that  if  communication  channel  continuity  depends  on  the 
availability  of  relay  server  nodes,  then  the  term  Jij{t„)  in  equation  (1)  is  also  a 
function  of  s^(f„)  and  s'^(f„)  for  each  server  making  up  the  channel.  This  is 
addressed  in  Table  1,  which  defines  Jij{t„)  (and  the  converse  continuity  from  the 
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Table  1.  Definitions  of  /y  (f«)  and  /;,  (t„)  for  all  possible  SA  communication  channels. 


Description  of 
Node  i 

Description  of 
Node  / 

Jij 

Jji  (tn) 

Any  client  in  local 
net  Xi  and  CSMA 
net  Yi 

Any  other  client 
in  local  net  Xi 
and  CSMA  net  Yi 

^AO 

^J.i 

SINCGARS-only 
client  in  local  net  X: 
and  CSMA  net  Yi 

SINCGARS-only 
client  in  local  net 
X2  and  CSMA  net 
Yi 

■^CSMA,X,  (^«)  ^ 

'^CSMA,X2  (^n)  ^ 

*^CSMA,X,  (^n)  ^ 

•^CSMA^Xj  i^n)  ^  ^J,i 

SINCGARS-only 
client  in  local  net  Xi 
and  CSMA  net  Yi 

SINCGARS-only 
client  in  local  net 
X3  and  CSMA  net 
Yz 

•^CSMAX,  ^  “^CSMA^  (^n  ) 

'^‘^MSG,Yi(^n)^^y,/ 

SINCGARS-only  j 

client  in  local  net  Xi  | 
and  CSMA  net  Yi 

EPLRS-equipped 
client  in  local  net 
X2  and  CSMA  net 
Yi 

‘^CSMA,X,  ^ 

•^CSMA^Xj 

j 

^CSMA,X,  (^«)  ^  ^j,i 

SINCGARS-only 
client  in  local  net  Xi 
and  CSMA  net  Yi 

EPLRS-equipped 
client  in  local  net 
X3  and  CSMA  net 
Y2 

*^CSMA^  (fn  )  *^CSMA^j  (fn  ) 

'^CSMA.X,  i^n  )  ^ 

•^MSG.Y,  (^n)  ^  ^7,/ 

SINCGARS-only 
client  in  local  net  Xi 
and  CSMA  net  Yi 

CSMA  server  in 
local  net  X2  and 
CSMA  net  Yi 

•^CSMA,X,  i^n)  ^^i,J 

•^CSMA,X,  ^ 

SINCGARS-only 
client  in  local  net  Xi 
and  CSMA  net  Yi 

CSMA  server  in 
local  net  X3  and 
CSMA  net  Y2 

•^CSMA,Xi  (^n  )  ^ 

•^MSaY,  (^n)  ^7,/ 

SINCGARS-only 
client  in  local  net  Xi 
and  CSMA  net  Yi 

Self  server  in 
CSMA  netYi 

•^CSMA,X,  (^n  )  ^ 

•^CSMA,X,  (^n  )  ^  ^J,i 

SINCGARS-only 
client  in  local  net  Xi 
and  CSMA  net  Yi 

Self  server  in 
CSMA  net  Y2 

*^CSMA,X,  (^fj) 

!  ‘^MSG.Yj  (^n)  ^ 

‘^CSMA,Xi  (^n  )  ^ 

-^MSG,  Y,  (^«  )  ^  ^j,i 

SINCGARS-only 
client  in  local  net  Xi 
and  CSMA  net  Yi 

MSG  server  in 
local  net  X4  and 

1  CSMA  net  Yi 

*^CSMA,X, 

'^CSMA,X«(^«)  ^^1,/ 

‘^CSMA,X,  (^n)  ^  ^7,1 

SINCGARS-only 
client  in  local  net  Xi 
and  CSMA  net  Yi 

MSG  server  in 
local  net  X5  and 
CSMA  netY2 

*^CSMA;C,  (^fl  )  ^  *^CSMA^j  (^n  ) 

^  “^MSG.Y, 

•^CSMA,X,  (^n) 

EPLRS-equipped 
client  in  local  net  Xi 
and  CSMA  net  Yi 

EPLRS-equipped 
client  in  local  net 
X2  and  CSMA  net 
Yi 

*^CSMA,X,  (^rt)  ^  ^J,i 

EPLRS-equipped 
client  in  local  net  Xi 
and  CSMA  net  Yi 

EPLRS-equipped 
client  in  local  net 
X3  and  CSMA  net 

^CSMKXy  i^n)  ^ 

■^MSG.Yj  (^n  )  ^  ^i,j 

‘^CSMA.X,  (^/i)  ^ 

“^MSG,  Y,  ^  ^7,/ 

Yz 


Table  1.  Definitions  of  Jij  (f„)  and  Jji  (tn)  for  all  possible  SA  communication  channels 
(continued). 


Description  of 
Node  z 

Description  of  Node 

j 

Jij  (^n) 

EPLRS-equipped 
client  in  local  net 
Xi  and  CSMA  net 
Yi 

CSMA  server  in 
local  net  X2  and 
CSMA  netYi 

‘^CSMA,X,  (4  )  ^  ^y,i 

EPLRS-equipped 
client  in  local  net 
Xi  and  CSMA  net 
Yi 

CSMA  server  in 
local  net  X3  and 
CSMA  netY2 

*^MSG,Y2  (^n)  ^ 

*^CSMA,X,  ^ 

•^CSMA,  Y,  ^ 

EPLRS-equipped 
client  in  local  net 
Xi  and  CSMA  net 
Yi 

Self  server  in 
CSMA  netYi 

*^CSMA,Xi  (^n  )  ^  ^j,i 

EPLRS-equipped 
client  in  local  net 
Xi  and  CSMA  net 
Yi 

Self  server  in 
CSMA  net  Y2 

^MSG.Yj  (^«  )  ^  ^ij 

‘^CSMA.X,  ^ 

*^MSG,Y,  ^  ^j\i 

EPLRS-equipped 
client  in  local  net 
Xi  and  CSMA  net 
Yi 

MSG  server  in 
local  net  X4  and 
CSMA  netYi 

'^CSMA.X,  (^n  )  ^ 

•^CSMA^X,  (Jn)  ^  ^J,i 

EPLRS-equipped 
client  in  local  net 
Xi  and  CSMA  net 
Yi 

MSG  server  in 
local  net  X5  and 
CSMA  netY2 

■^CSMA,  X,  (^<1 ) 

‘^CSMA.X,  (^n)  ^ 

*^MSG,Yj  (^n)  ^  ^j\i 

CSMA  server  in 
local  net  Xi  and 
CSMA  net  Yi 

hj 

5-^ 

CSMA  server  in 
local  net  Xi  and 
CSMA  net  Yi 

CSMA  server  in 
local  net  X3  and 
CSMA  netY2 

^MSG,Y2  ^ 

•^MSaY,  (^«)  ^  ^J,i 

CSMA  server  in 
local  net  Xi  and 
CSMA  netYi 

Self  server  in 
CSMA  netYi 

ZAO 

CSMA  server  in 
local  net  Xi  and 
CSMA  net  Yi 

Self  server  in 
CSMA  net  Y2 

•^MSG,  Yj  (^n  )  ^  ^ij 

*^MSG,  Y,  i^n  )  ^  ^j\i 

CSMA  server  in 
local  net  Xi  and 
CSMA  netYi 

MSG  server  in 
local  net  X4  and 
CSMA  netYi 

'^CSMA,X4  i^n  )  ^ 

s:AO 

5y./ 

CSMA  server  in 
local  net  Xi  and 
CSMA  netYi 

MSG  server  in 
local  net  X5  and 
CSMA  netY2 

■^CSMA,  X5 

^MSG.Y,  (^J  ^ 

*^MSG,Y,  (K)  ^  ^j,i 

Self  server  in 


Self  server  in 


Table  1.  Definitions  of  /,y  (f„)  and  Jji  (tn)  for  all  possible  SA  communication  channel 
(continued). 


Description 
of  Nodei 

Description  of  Node/ 

Jijitn) 

Jjiitn) 

Self  server  in 
CSMA  net  Yi 

Self  server  in 

CSMA  netY2 

^MSG^yS^h) 

•^MSG.Y, 

Self  server  in 
CSMA  net  Yi 

MSG  server  in 
local  net  X4  and  CSMA 
net  Yi 

^CSMA.xS^n)  ^ 

55  AO 

5;.i 

Self  server  in 
CSMA  net  Yi 

MSG  server  in 
local  net  Xs  and  CSMA 
net  Y2 

^CSMA.X^  ^  ^ 

^MSG,Y2  ^  ^  ^ 

•^MSG.Yi  (^rt)  ^ 

MSG  server 
in 

Local  net  X4 
and  CSMA 
net  Yi 

MSG  server  in 
i  local  net  X5  and  CSMA 
net  Y2 

^CSMA,Xs  ^  ^ 

^MSGJz  ^  ^ij  ^ 

‘^CSMA.X^  (^n  )  ^ 

“^MSG.Yi  ^ 

Notes: 

•^CSMA  X  )  •  communication  capability  state  of  the  CSMA  server  in  local  SINCGARS  net  X„. 


^msg,  y  )  •  communication  capability  state  of  the  MSG  server  in  CSMA  net  Y^. 

:  delta  function  representing  the  inclusion  of  the;*  node  within  the  i*  node's  AO  neighborhood. 
a:  binary  AND  operator. 

fth  to  y*  node  Jjiltn])  for  all  possible  SA  communication  channels.  As  indicated  in 
this  table,  an  SA  communication  channel  is  built  using  the  following  generic 
elements. 

•  CSMA  servers.  The  binary  communication  capability  state  of  the  CSMA 
server  in  local  SINCGARS  net  Xn  at  time  step  t„  is 

■^CSMA.X,  ~  “^CSMA.X,  '^'^CSMA.X,  '  WhetO  )  and  ‘S'csMA.X. 

the  transmission  and  reception  capability  states  of  the  server,  respectively. 

•  MSG  servers.  The  communication  capability  state  of  the  MSG  server  in 

CSMA  net  Ym  at  time  step  t„  is  =  where 

‘^MSG,Y„(^/.)^^®  tTansmission  and  reception  states  of  the 
server,  respectively. 

•  AO  neighborhoods.  The  binary  metric  dfj  is  a  delta  function  representing 

the  inclusion  of  the  node  within  the  i*  node's  AO  neighborhood,  and  is 
defined  as 

J'^  within  the  cell’s  AO  neighborhood 

'■'  [  0  otherwise 
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In  these  expressions,  the  symbol  a  is  the  binary  AND  operator. 

It  is  interesting  to  consider  the  various  scales  over  which  the  SA  connectivity  is 
calculated  for  the  nodes  within  the  CA  brigade  (Figure  3).  Given  the  total 
number  of  SA  transmitting  neighbor  nodes  N,  relative  to  the  node  cell  (as 
defined  in  equation  [1]),  one  can  define  the  comiectivity  distribution  P(N),  which 
gives  the  probability  that  a  randomly  selected  reference  node  has  N  neighboring 
nodes  within  its  AO  neighborhood  (where  N  =  0, 1,  2,  611).  Figure  6  displays 

the  P(N)  distribution  for  the  CA  brigade;  this  conveys  the  range  of  scales  over 
which  nodes  receive  SA  messages.  Note  that  the  higher  the  number  of  nodes 
within  a  platform's  AO  neighborhood,  the  more  likely  it  is  to  be  a  higher-echelon 
platform. 

At  this  point,  the  neighbor  nodes  within  a  center  reference  cell's  AO 
neighborhood  do  not  directly  interact  with  the  reference  cell  as  part  of  its  cell 
transition  fimction.  At  it  is  used  here,  the  term  "transition"  explicitly  refers  to  a 
change  in  the  reference  cell's  functional  transmission  and  reception  capabilities. 
In  other  words,  the  fmictional  states  of  node  cells  evolve  independently  of  one 
another;  what  the  CA  model  still  lacks  at  this  point  in  our  discussion  is  a  type  of 
neighbor  cell  that  directly  interacts  with  a  reference  node  cell  to  perturb  its 
commimication  capabilities  through  lO  stress.  This  type  of  threat  cell  is 
introduced  in  the  next  section. 


Figure  6.  Connectivity  distribution  P(N)  for  the  CA  brigade. 
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2.5  Introducing  Information  Operations  Stress 

The  dynamics  of  the  CA  model  are  now  added  through  the  introduction  of  an  lO 
stress,  which  directly  interacts  with  a  node  ceU  in  order  to  perturb  its 
communication  capabilities.  Interaction  of  lO  stress  with  nodes  can  result  in 
either  (a)  temporary  dysfunction  of  a  node  (where  s,’'(f„)  =  0,  Si^(f„)  =  0,  or  s,^(f„) 
=  Si^{t„)  =  0),  (b)  permanent  dysfunction  of  the  node,  or  (c)  "spoofing"  of  the 
global  positioning  system  (GPS)  receiver  operating  onboard  the  node  platform. 
This  last  type  of  effect  forces  the  impacted  node  to  laxmch  an  unscheduled  SA 
update  message  (which  is  then  routed  throughout  die  brigade  via  the  network 
architecture  depicted  in  Figiue  1);  in  this  case,  both  Si^{tn)  and  remain  equal 
to  1.  lO  stress  can  arise  either  from  a  stochastic  parameter  embedded  within  a 
nodal  state  vector  or  an  external  neighboring  cell  (representing  a  localized  lO 
threat  source)  that  generates  a  local  "stress  field"  (these  two  different  tjrpes  of 
stress  will  be  explicitly  modeled  in  section  3). 

The  dynamics  which  emerge  from  node  cell/IO  stress  interactions  are  modeled 
as  perturbations  to  a  node's  communication  capability  states.  Given  that 
equation  (1)  defines  the  i*  node's  SA  connectivity  within  a  stress-free 
environment,  then  the  stress-perturbed  SA  connectivity  of  that  node  at  time  step 
tnis 


In  this  equation,  s,^  (fn-i),  s/(f„-i),  and  /,>  (f„-i)  are  the  i*  node's  prior  reception  and 
node's  prior  transmission  capability  states,  and  the  to  z*  node  channel 
continuity  at  time  step  t„.i,  respectively,  and  As,R  (t„).  As/  (t„),  and  A/y  (f„)  are 
perturbations  (induced  at  the  current  time  step  f„)  to  the  z*  node's  prior  reception 
and  y*  node's  prior  transmission  capability  states,  and  the  to  z*  node  prior 
channel  continuity,  respectively.  The  last  metric  reflects  tiie  introduction  (and 
possible  subsequent  removal)  of  stress-related  perturbations  to  those 
neighboring  server  nodes  connecting  the  z*  andj*  nodes  (as  defined  in  Table  1). 
Each  of  the  metrics  Asjt'^(f„),  Asi^  (f„),  and  A/,)  (f„)  may  assume  values  of  either  -1  or 
1,  reflecting  the  loss  or  return  of  communication  capabilities,  respectively.  Thus, 
equation  (2)  models  a  process  of  dynamic  response,  where  stress-induced 
perturbations  can  affect  nodes  either  directly  (s,^  [f„.i]  +  As,R  [f„]  =  0)  or  indirectly 
lUij  (fn-l)  +  A/y  (f„)]  =  0  or  [S;^(tn.l)  +  As/(fn)]  =  0). 

2.6  Analyzing  the  Emergent  Network  E)5mamics 

There  are  several  data  processing  filters  that  have  been  coded  into  the  CA  model 
for  analyzing  time  series  data.  In  this  section,  these  filters  are  described  in  detail. 
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2.6.1  Normalized  Network  Stress 

The  first  data  filter  calculates  the  normalized  perturbative  stress  that  is 

introduced  into  the  entire  brigade  system  via  localized  lO  stress.  This  can  be 
expressed  as 

where  M  is  the  total  number  of  node  platforms  within  the  digitized  brigade  and 
hiStress  indicates  a  stress-induced  perturbation  of  the  z*  node's  commtmication 
capability  states  at  time  step  t„.  Thus,  (/zs*«ss  (f„)>  indicates  the  fraction  of  stress- 
impacted  nodes  as  a  function  of  time.  For  the  "movement  to  contact"  brigade 
structure  presented  in  section  2.3.2,  M  =  612.  It  should  be  noted  that  while 
(f„))  reflects  interactions  between  lO  stress  sources  and  target  node  cells,  Ci{t„),  on 
the  other  hand,  reflects  inter-ceU  dependencies  which  directly  contribute  to  a 
node's  SA  connectivity  as  defined  in  equation  1.  Thus,  node  cell/ 10  stress 
source  interactions  are  distinguished  apart  from  reference  node/  neighbor  node 
topological  interdependencies. 

2.6.2  Globally  Averaged  Situational  Awareness  Coimectivity  <C(t„)) 

The  second  data  filter  calculates  globally  averaged  SA  connectivity  as  a  function 
of  scenario  time,  where  an  average  value  is  calculated  at  each  point  in  discrete 
simulation  time  across  all  node  platforms  within  the  CA  brigade.  The  dynamic 
globally  averaged  SA  connectivity  for  the  CA  brigade,  (C(fn)>,  is  defined 

r  M 

where  again  M  =  612  is  the  total  number  of  nodes  within  the  brigade. 


2.6.3  Situational  Awareness  Fitness 

Given  that  the  total  SA  connectivity  state  of  the  M-node  CA  brigade  can  be 
represented  as  an  M-tuple  point  within  a  unit  hypercube  (which  can  be  thought 
of  as  representing  the  global  system  "phase  space"),  the  hypercube  axes  (where 
the  z*  axis  represents  the  SA  connectivity  level  C,(fn)  of  the  z'*  node)  can  be 
"coarse-grained"  (Wolfram  1985)  or  partitioned  into  interval-based  fitness  states. 
Figure  7  defines  how  a  five-state  partition  (used  by  Guzie  for  vulnerability  risk 
assessment  [Guzie  2000])  is  applied  for  evaluating  the  SA  connectivity  Ci{t„)  of 
the  z*  node;  this  partition  divides  the  unit  interval  into  the  interval-based  fitness 
states  of  very  high  (blue),  high  (green),  medium  (yellow),  low  (red),  and  very  low 
(orange).  The  assigned  colors  are  used  to  visually  convey  a  node's  fitness  state 
during  a  simulation  run  (which  will  be  discussed  in  the  subsequent  section  of 
this  report).  This  has  the  effect  of  dividing  the  unit  hypercube  continuxim  into  5^ 
discrete  blocks,  where  each  block  represents  a  unique  global  system  state. 
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VERY  HIGH 

CM  >  0.9 

HIGH 

0.6  <  C,(g  <  0.9 

MEDIUM 

0.4  <  C,.(g  <  0.6 

LOW 

0.1  <  c,(g<o.4 

VERY  LOW 

11 

Q(g<o.i 

Figure  7.  Application  of  Guzie's  five-state  partition  to  the  SA  connectivity  of  the  i"’  node. 

For  the  612-node  brigade  structure,  there  are  global  system  states.  Figure  8 
illustrates  how  coarse-graining  the  SA  connectivity  levels  for  each  of  the  612 
brigade  nodes  modifies  the  global  system's  dynamic  state  vector  (where  Ci 
through  C622  indicate  the  SA  connectivity  for  nodes  1-612,  respectively)  as  the 
system  is  perturbed  via  interactions  with  lO  stress  sources.  Here,  continuous 
levels  of  nodal  SA  connectivity  are  mapped  into  one  of  five  fitness  states  at  each 
subsequent  time  step  of  a  simulation  rrm.  Finally,  Pfuness{tn)  is  defined  as  the 
fraction  of  nodes  in  a  particular  fitness  state  at  time  step  tn. 


to  t, 
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Figure  8.  Dynamics  of  the  SA  connectivity  state  vector  for  the  CA  brigade. 
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The  third  data  filter  coded  into  the  CA  model  calculates  the  five  fitness  metrics— 

Pvery  highi^i^/  Phigh(it^/  Pmedium(tt^/  Piow(it^/  and  Pvery  law  (fn)*  These  metrics  provide  a 
different  perspective  relative  to  (C{tn)).  They  illustrate  the  distribution  of  various 
situational  awareness  levels  across  individual  platform  nodes,  rather  than  refer 
to  a  global  average.  Thus,  when  the  five  metrics  Pvery  high{tn),  Phigh{tn),  Pmedium{t„), 
Ptow{tn),  and  Pvery  low  (tn)  (where  Pvery  high{tn)  Phigh(tn)  Pmedium{tn)  Plow(tn)  "t  Pvery  low 

(tn)  =  1  for  all  tn)  are  taken  in  combination,  they  offer  a  concise  way  to  represent 
the  612-tuple  system  state  (one  of  the  discrete  blocks  in  the  "phase  space" 
hypercube)  at  f„.  It  should  be  noted,  however,  that  there  can  be  many  different 
global  system  states  represented  by  the  same  5-tuple  combination  [Pvery  high{tn), 

Phigb(tr^f  Pmediumitj^ f  Plow{tn)/  Pvery  low  (f«)]* 


2.6.4  Course-Grained  Entropy 

The  fitness  states  introduced  in  the  previous  section  partition  the  M-tuple  global 
SA  connectivity  system  state  at  time  step  t„  into  five  mutually  exclusive  sets  of 
states,  where  the  probability  of  occupying  the  /c*  fitness  state  is  Pk{tn)  and 

(5) 

k=l 

These  conditions  define  a  complete  finite  probability  scheme  based  on  the  coarse 
graining  of  nodal  SA  connectivity  measures  into  discrete  states.  Given  such  a 
scheme,  the  global  state  entropy  H{t„)  is  defined  as  a  quantitative  measure  of  die 
disorder/volatility  within  the  global  system  at  time  step  t„,  where 

=  log, (P, (t. >)  .  (6) 

k=I 

The  state  entropy,  which  is  based  on  the  Shannon  entropy  of  an  information 
source,  can  range  from  a  level  of  0  (indicating  that  all  nodes  are  in  the  same 
fitness  state)  up  to  a  level  of  1  (aU  nodes  are  equally  distributed  among  all 
possible  fitness  states). 

The  concept  of  state  entropy  can  be  extended  in  order  to  quantify  the  temporal 
disorder  between  system  states  measured  at  successive  time  steps  Ua  and  t„. 
Thus,  the  temporal  entropy  H{t„,  t„.i),  based  on  the  conditional  Shannon  entropy  as 
measured  between  an  information  source  and  receiver,  is  defined  as 

m..t,^,)  =  -'t'tp,(l.>rP(k.l.lJ.l..,)Aog,(P(k.l.\J.I.-,))  ■  (7) 

J=I  k=I 

In  this  equation,  P(k,  t„  |  j,  f„-i)  is  the  conditional  probability  that  a  node  in  the/* 
fitness  state  at  time  step  f„-i  maps  into  the  k*  fitness  state  at  time  step  U.  As  with 
state  entropy,  the  temporal  entropy  can  range  from  a  level  of  0  (indicating  a  static 
condition  where  all  nodes  remain  in  the  same  fitness  state  from  t„.\  to  t„)  up  to  a 
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level  of  1  (indicating  a  maximally  volatile  condition  where  nodes  are  equally 
likely  to  map  into  any  one  of  the  possible  fitness  states  from  f„-i  to  tn). 

2.6.5  Return  Map 

Once  the  Pfitness{tn)  time  series  have  been  calculated  for  an  lO  scenario,  the  data 
within  each  time  series  can  be  analyzed  for  certain  structural  characteristics.  This 
is  done  through  the  use  of  a  return  map,  which  is  a  2-D  scatter  plot  of  Pfitness{tnn) 
vs.  Pfitness{tn)  for  each  fitness  state.*  A  return  map  illustrates  the  dynamics  of  a  1-D 
time  series  (i.e.,  a  time  series  with  just  one  independent  variable)  in  a  recursive 
rather  tihan  successive  format.  The  utility  of  a  return  map  for  portraying  time 
series  structure  is  illustrated  in  Figure  9,  which  portrays  the  dynamics  of  tire 
chaotic  1-D  logistic  map  x„+i  =  4x„(l-  x„)  in  both  standard  time  series  (Figure  9[a]) 
and  return  map  (Figure  9[b])  formats.  In  this  figure,  the  first  plot  appears  to 
depict  a  random  time  series  with  values  in  the  interval  [0, 1],  while  the  second 
plot  (displaying  x„+i  vs.  x„)  clearly  indicates  the  presence  of  an  attractor  (a  subset 
of  one  or  more  phase-space  points  of  the  form  [x„,  x„+i]  towards  which  a 
deterministic  dynamical  process  wiU  evolve  or  is  "attracted"). 


3.  Information  Operations  Simulations 


In  this  section,  two  different  models  of  lO  stress  are  introduced,  and  simulation 
results  of  CA  model  runs  incorporating  these  stresses  are  presented  and 
discussed.  Specifically,  these  stress  models  address  (a)  a  "mean  field"  stress 
embedded  within  a  node  cell  and  (b)  localized  stress  fields  generated  by  source 
cells  external  to  node  cells. 

3.1  Mean  Field  Stress 

The  first  type  of  lO  stress  was  modeled  as  a  stochastic  process  which  acts  to 
perturb  the  functional  communication  states  of  each  node  within  the  SA 
network.  This  lO  stress  arises  from  a  random  binary  variable  representing  nodal 
digital  coirununication  capabiHty  (i.e.,  a  node  either  can  [1]  or  cannot  [0] 
communicate  with  other  nodes,  where  "communicate"  implies  both  transmission 
and  reception  subcapabilities).  Such  a  random  variable  is  embedded  within  each 
nodal  state  vector,  which  are  synchronously  updated  for  each  node  in  the 
network.  Since  these  distributed  random  variables  evolve  independently  from 
one  another,  there  is  no  spatial  correlation  between  the  stress  at  different  nodes. 


*  An  equivalent  version  of  the  return  map  would  plot  Pfitness{in)  vs.  Pfitness{tnA)- 
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This  t5^e  of  lO  stress  is  similar  to  the  mean  field  approximation  from  statistical 
physics,  where  every  individual  element  within  a  multi-element  system  is 
subjected  to  the  average  influence  of  a  global  field  interacting  with  the  rest  of  the 
system  (Bar-Yam  1997).  In  this  sense,  the  "mean  field"  lO  stress  is  at  best  a  first- 
order  approximation  to  a  real  multi-threat  lO  environment,  serving  to 
demonstrate  the  emergent  behavior  of  the  brigade's  SA  network  topology  under 
a  coherent  hypothetical  rather  than  realistic  stress. 

The  dynamics  of  this  mean  field  stress  are  described  by  a  stationary  Markov 
chain  which  is  applied  uniformly  and  independently  to  each  of  the  nodes  within 
the  digitized  brigade.  This  Markov  chain  can  be  written  as  the  iterative  map 


^Si  (K  )  ~  ^s,  (K-l  ’^n-I  )’ 


(8) 


where 

p(s;(>.))  p(s!(<.))  -  p(sf<t.))]  (9) 

and 

p„(t.-,MpU(t.-,))  p(^;(‘.-,))  p^d..,))  -  p(s!-(t,_, ))]  (10) 


are  row  vectors  denoting  state  probabilities  for  the  network  node  at  time  steps 
t„  and  f„-i,  respectively,  and  X(tn,  U-i)  is  a  stationary  one-step  transition  matrix. 
The  notation  P(sHtn))  indicates  the  probability  that  the  i*  node  occupies  the  A*  of 
L  possible  states  at  time  step  tn,  and  s,Ktn)  =  sA^(fn)]  is  the  i*  node's 

communication  state  vector  consisting  of  transmission  and  reception  capability 
states,  respectively. 


The  stationary  transition  matrix  A.(f«,  U-i)  introduced  in  equation  (8)  can  reflect 
either  transient/ permanent  nodal  communication  dysfunction  or  spoofing  of  the 
GPS  receiver  mounted  within  a  node,  causing  the  node  to  launch  an  unscheduled 
SA  update  message  across  the  network.  The  former  form  of  the  transition  matrix 
is  expressed  as 


p(s^(tJ\s^(t„_j)) 

J I  s'^( ))  p(s^( )  I  s^( t^j  )}  p(s^^( t„  )  I  s‘^( j) 


,  (11) 


while  the  latter  form  is  expressed  as 


pisUij^s'd.-,)) 


p{s^(i.>\sUt..,))\  ' 


(12) 
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where  s  s  s  and  s  s(f„)  represent  the  binary  nodal  communication 
states  functional,  transient  dysfunction,  permanent  dysfunction,  and  fimctional 
but  spoofed,  respectively.  Thus,  the  stochastic  iterative  map  in  equation  (8)  can 
describe  the  transitional  probability  for  each  network  node  between  fimctional 
(S(^(tn)  =  =  1)  and  transient  or  permanent  dysfimction  (where  s,^(tn)  =  Si^{tn) 

=  0  in  both  cases),  or  between  functional  and  fimctional  but  spoofed  states. 

Since  the  mean  field  stress  model  considers  the  impact  of  network  flooding  on 
nodal  connectivity,  the  dynamics  of  SA  message  dissemination  must  be 
incorporated  into  the  model.  If  it  is  assumed  that  (1)  a  node  platform  must 
report  its  new  position  every  time  it  moves  a  distance  of  100  m  and  (2)  aU  nodes 
within  the  brigade  are  moving  forward  at  a  constant  speed  of  3  km/hr,  then  each 
node  automatically  sends  out  a  new  SA  message  once  every  120  s.  The 
dynamics  of  SA  message  dissemination  throughout  (he  CA  brigade  is  based  on 
the  process  used  in  the  Next  Generation  Performance  Model  (NGPM)  developed 
by  the  U.S.  Army  Communication  Electronics  Command  (CECOM)  (U.S.  Army 
Communication  Electronics  Command  1999);  see  Appendix  A  for  details  on  this 
process.  Finally,  in  order  to  prevent  flooding  due  to  automatic  SA  message 
reporting,  initial  reporting  times  between  0  and  120  s  are  assigned  to  each  node 
from  a  uniform  random  distribution. 

When  applying  a  mean  field  stress  as  described  in  equation  (8),  the  SA 
connectivity  described  in  equation  (1)  must  be  modified  to  the  form 

J 

where  s,(f„)  and  S;(f„)  are  the  binary  communication  capability  states  (i.e., 
communicate/can't  communicate)  of  the  i*  and  y*  nodes  at  time  step  t„, 
respectively.  In  equation  (13),  the  communication  channel  continuity  between 
the  i*  andy*  nodes  is  now  expressed  as 

=  (14) 

where  the  binary  channel  continuity  metric  6,y<’"'’"“''y(f„)  =  1  only  if  the  charmel 
transmission  latency  L,;  (the  total  amount  of  time  required  for  an  SA  message  to 
travel  between  the  and  y*  nodes)  is  less  than  a  preset  threshold  limit  Lthreshou- 
The  situation  where  8y""i‘"“''y(fn)  =  0  (i.e.,  L,)  >  Lthresboid)  arises  due  to  excessive 
message  queue  loading  of  the  servers  making  up  the  communication  channel. 
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This  process  is  analogous  to  balls  being  thrown  into  a  funnel  with  a  trap  door 
(Figure  10).  The  commuiiication  channel  connecting  the  to  the  node  is 
represented  as  a  wide-mouthed  fuiuiel,  which  can  collect  and  disseminate  SA 
messages  from  multiple  transmitting  nodes  (including  the  node)  in  order  to 
route  the  messages  to  the  node.  Each  incidence  of  an  incoming  ball  (i.e.,  an  SA 
message)  contributes  an  additional  member  to  a  stack  of  balls  (i.e.,  queued  SA 
messages)  waiting  to  be  released  through  a  trap  door  at  the  bottom  of  the  fumael. 
This  trap  door  opens  (at  a  fixed  frequency)  just  long  enough  to  allow  one  ball  to 
exit  before  closing  again,  representing  a  constant  outgoing  message  baud  rate  for 
the  commimication  channel.  As  the  funnel  fills  up,  the  top  of  the  stack  of  balls 
reaches  and  then  surpasses  a  threshold  level  (i.e.,  Uhreshou)-  At  this  point,  the 
channel  continuity  metric  =  0,  meaning  that  any  new  SA  message 

coming  in  from  the  node  (i.e.,  the  red  ball  at  the  top  of  the  figure)  will  not 
arrive  at  the  f*  node  in  time  to  be  useful.  Finally,  =  1  again  only 

when  the  stack  level  in  the  funnel  falls  below  the  threshold  latency  level. 


Incoming  SA 
messages  from 
other  nodes. 


^threshold 


Figure  10.  Representation  of  an  inter-nodal  communication  channel  as  a  funnel 
collecting  and  dispersing  balls  through  a  trap  door. 
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Once  the  form  of  the  mean  field  stress  and  associated  SA  coimectivity  have  been 
established,  the  CA  model  simulating  this  type  of  lO  stress  can  be  run.  Figure  11 
depicts  globally  averaged  SA  connectivity  (averaged  over  10  simulation  nms 
using  different  random  number  seeds)  emerging  across  the  brigade  network 
under  exposure  to  the  form  of  mean  field  stress  that  induces 
transient/permanent  nodal  commtmication  dysfunction.  For  these 
simulations,  the  transitional  probabilities  in  equation  (11)  are  set  to 


and 


■)^dysfmcfy  .  ; 


0.99899  0.00100  0.00001 
0.29999  0.70000  0.00001 
0  0  1.00000 


0.9899  0.0100  0.0001 
0.2999  0.7000  0.0001 
0  0  1.0000 


(15) 


(16) 


where  (i„ ,  ^„_, )  and  (^„ ,  )  are  transition  matrices  associated  with  low 

and  high  levels  of  nodal  dysfunction  stress,  respectively.  Also,  tiie  channel 
latency  threshold  Lthresiwu  is  set  to  300  s  (5  min),  so  that  a  communication  channel 
is  considered  broken  if  the  difference  between  the  transmitted  and  actual 
positions  of  the  reporting)*  node  is  >  500  m.  The  figure  readily  shows  that  both 
low  and  high  levels  of  this  first  type  of  mean  field  stress  results  in  a  roughly 
linear  decrease  in  average  SA  connectivity.  This  is  due  to  the  coherent  nature  of 
the  mean  field  stress  as  it  is  applied  throughout  the  brigade  structure,  and  also  to 
the  lack  of  spatial  correlation  between  the  functional  states  of  network  nodes. 

Next,  simulations  are  run  using  the  second  type  of  mean  field  stress,  which 
induces  GPS  spoofing  and  subsequent  network  flooding.  Figure  12  depicts 
globally  averaged  SA  connectivity  (again  averaged  over  10  simulation  nms  using 
different  random  number  seeds)  emerging  across  the  brigade  network  under 
exposure  to  GPS  spoofing  stress.  For  these  simulations,  the  transitional 
probabilities  in  equation  (12)  are  set  to 


0.99  0.01 
0.30  0.70 


(17) 
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steady  State  Mean  Field  Stress  —  Modal  Dysfunction 


Simulation  Time  (seconds) 

Figure  11.  Average  SA  connectivity  for  the  brigade  network  exposed  to  a  mean  field 
stress  which  induces  nodal  transient/ permanent  dysfunction. 


steady  State  Mean  Field  Stress  —  GPS  Spoofing 


Simulation  Time  (seconds) 

Figure  12.  Average  SA  connectivity  for  the  brigade  network  exposed  to  a  mean  field 
stress  which  induces  nodal  GPS  spoofing. 


and 


0.9  0.1 
0.3  0.7 


(18) 


where  and  ,t^  are  transition  matrices  associated  with  low 

and  high  levels  of  GPS  spoofing  stress,  respectively.  Again,  the  channel  latency 
threshold  LthresJwu  is  set  to  300  s.  Comparing  this  figure  with  Figure  11,  it  is  clear 
that  the  message  flooding  which  results  from  spoofing  is  tihe  dominant  cause  of 
SA  connectivity  degradation  when  considering  both  types  of  mean  field  stress. 
Both  the  low  and  high  levels  of  spoofing  stress  cause  the  average  SA  connectivity 
to  decrease  in  staggered  bursts  of  activity  (a  phenomenon  known  within 
complexity  science  as  "punctuated  equilibrium"  [Bak  and  Sneppen  1993;  Bak 
1996;  Newman  and  Sneppen  1996])  xmtil  a  quasi-steady-state  condition  is 
reached  at  t„  =  1850  -  1900  s.  At  this  point,  server  message  queues  are 
sufficiently  loaded  so  that  the  transmission  latency  for  all  susceptible 
communication  channels  achieves  or  surpasses  LmreshoM,  and  will  then  likely 
remain  in  this  state  if  the  simulation  timeframe  is  extended  indefinitely. 

Although  the  mean  field  model  serves  to  comparatively  demonstrate  different 
forms  of  lO  stress,  it  has  shortcomings  in  two  different  areas.  First,  the  model 
can  significantly  underpredict  the  impact  of  GPS  spoofing  on  nodal  SA 
connectivity.  Since  the  brigade  C2  network  structure  has  deliberately  been 
decoupled  from  the  SA  network  structure  as  coded  within  the  CA  model,  the 
cascading  decrements  in  average  SA  connectivity  depicted  in  Figure  12  can  only 
reflect  SA-message-induced  server  queue  saturation.  It  is  likely  that  the  addition 
of  a  superimposed  C2  structure  wiU  act  to  accelerate  server  queue  saturation 
(where,  for  the  majority  of  nodes,  the  chaimel  continuity  metric  will 

jump  to  a  value  of  zero  much  sooner)  due  to  C2  message  contention.  Thus,  the 
network  response  shown  in  Figurel2  can  at  best  serve  only  as  a  lower  bound 
when  trying  to  predict  the  negative  impact  of  GPS  spoofing  on  nodal  SA 
connectivity. 

A  more  serious  shortcoming  of  the  mean  field  model  is  its  failure  to  realistically 
consider  correlations  between  spatially  dispersed  stress  sources,  as  well  as  the 
resultant  correlations  between  SA  connectivity  levels  associated  with  impacted 
nodes.  In  the  mean  field  approximation  of  lO  stress,  stress  sources  are  modeled 
as  random  variables  embedded  within  nodal  state  vectors.  These  random 
variables  represent  stochastic  perturbations  to  a  node's  communication 
capability  state  which  are  locally  constrained  to  that  node.  Within  this  context, 
there  is  no  imposed  spatial  correlation  between  the  perturbations  induced  in 
different  network  nodes.  In  other  words,  there  is  no  stress  source  existing 
external  to  and  apart  from  a  node.  In  order  to  better  analyze  the  emergent 
network  behavior  induced  by  a  more  realistic  lO  threat,  stress  sources  must 
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be  introduced  which  can  model  an  external  perturbative  field  that  might  impact 
multiple  nodes  simultaneously.  This  enhanced  CA  model  is  discussed  in  the 
next  section. 

3.2  Localized  External  Stress  Fields 
3.2.1  Jammer  Bomb  Scenario 

A  more  reahstic  form  of  lO  stress  arises  from  perturbations  to  nodal 
communication  due  to  localized  external  fields  that  are  generated  by 
electromagnetic  (EM)  sources.  In  this  second  stress  model,  the  CA  dynamics  are 
added  by  introducing  spatially  dispersed  source  cells  representing  generic  radio 
frequency  (RF)  "jammer  bombs"  that  act  to  locally  jam  platform-mounted  radios 
passing  by  (Bothner  2000).  In  this  case,  the  stress  local  to  the  i*  node  at  time  U  is 

hr“(l.)  =  SUp(,B,(IJ,  (19) 

k 

where  "sup"  is  the  maximum  operator  appMed  to  aU  neighboring  lattice  cells  k 
relative  to  the  i*  node,  the  binary  metric  Bk(t„)  =  1  only  when  a  jammer  bomb 
occupies  the  fc*  cell  relative  to  cell  i  at  time  t„,  and  is  a  binary  coupling 
coefficient  between  the  i*  and  fc*  lattice  sites  and  is  equal  to  1  only  when  the  fc* 
site  is  within  an  isotropic  jamming  neighborhood  relative  to  tiie  i*  site  (i.e.,  the 
Euclidean  distance  from  the  centers  of  ceU  i  to  cell  k  dik  <  3.2  cell  lengths). 
Figure  13  illustrates  the  36-cell  jamming  neighborhood  (i.e.,  a  center  cell  plus  36 
neighboring  cells)  as  implemented  within  the  CA  model.  Then,  the  functional 
receive  state  of  the  i*  node  is 

=  (20) 

Thus,  a  node  is  jammed  if  it  passes  within  the  jamming  neighborhood  of  a 
jammer  bomb;  alternately,  a  node  is  jammed  if  and  only  if  there  is  at  least  one 
jammer  bomb  within  an  identical  jamming  neighborhood  relative  to  the  target 
node  (where  the  target  node  is  at  the  center  of  the  translated  neighborhood). 
The  latter  condition  is  necessary  since  the  jammer  bombs  are  designed  to 
advance  (downward)  toward  the  brigade  at  a  uniform  steady  rate  of  one  cell  per 
time  step;  this  simulates  the  same  relative  movement  as  if  the  vehicle  nodes  were 
advancing  (upward)  towards  the  stationary  bombs.  Since  the  localized  jammer 
bomb/node  interactions  act  to  perturb  nodal  reception  capabiHty  states,  this  in 
turn  perturbs  nodal  SA  connectivity  levels  (as  defined  in  equation  [2]) 
throughout  the  CA  brigade. 
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Figure  13.  Jamming  neighborhood  for  a  generic  RF  jammer  bomb. 

3.2.2  Simulation  Results 

Once  the  CA  model  (which  includes  the  cellular  brigade  configuration,  rules 
defining  cellular  interdependencies  required  for  nodal  SA  connectivity 
calculations,  and  additional  rules  defining  local  jammer  bomb/node  cell 
interactions)  has  been  defined  and  coded,*  the  parameters  that  delineate  an  lO 
simulation  can  be  specified.  First,  a  simulation  timescale  must  be  defined. 
Given  that  the  CA  lattice  has  already  been  scaled  to  125  m/cell  length,  only  the 
uniform  speed  at  which  the  brigade  advances  (remember  that  nodes  do  not 
move  relative  to  each  other,  implying  a  imiform  rate  of  forward  movement)  must 
be  specified.  If  it  is  assumed  that  the  uniform  platform  speed  is  10  km/hr,  then 
1  simulation  time  step  =  125  m  /  (10,000  m  /  3,600  s)  =  45  s. 

Next,  the  process  of  introducing  jammer  bombs  into  the  scenario  is  described. 
In  the  current  application  of  the  model,  we  define  a  two-phase  scenario,  which 
entails  (1)  the  digitized  brigade  encountering  and  moving  through  a  field  of 
jammer  bombs  previously  planted  by  hostile  forces,  and  (2)  hostile  unmanned 
aerial  vehicles  flying  over  the  advancing  brigade  to  deliver  a  barrage  of 
additional  jammer  bombs.  Figure  14  illustrates  the  pattern  of  cells  within  the  CA 
lattice  which  are  designated  as  either  potential  preplanted  bomb  sites  or  aerial 
barrage  targets. 


*  See  Appendix  B  for  a  listing  of  the  jammer  bomb  scenario  Cellang  source  code,  and  Appendix 
C  for  a  listing  of  the  data  filter  C  source  code  that  was  co-compiled  with  the  Cellang  code. 
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Figure  14.  Potential  jammer  bomb  site/target  cells  within  the  CA  brigade. 

The  first  phase  of  the  lO  scenario  commences  at  fn  =  0  (i.e.,  the  start  of  the 
simulation),  and  then  continues  to  unfold  until  tn  =  60  min.  During  this  time, 
jammer  bombs  are  inserted  into  the  row  of  dark  gray  cells  labeled  "preplanted 
bomb  sites"  at  the  top  of  Figure  14,  where  the  probability  of  finding  a  jammer 
bomb  at  a  cell,  given  that  the  cell  is  a  potential  bomb  plant  site,  is  P(bomb  I  plant 
site)  =  0.2 /time  step,  and  the  average  preplanted  bomb  population  within  the  CA 
lattice  is  480  bombs/simulation  trial.  Once  a  jammer  bomb  has  been  inserted 
into  a  cell,  the  bomb  moves  forward  into  the  brigade  at  a  constant  rate  of  1  cell 
length/time  step  (125  m/45  s),  and  continues  to  jam  friendly  node  platforms 
which  penetrate  its  jamming  neighborhood.*  This  is  continued  until  either  the 
bomb  is  destroyed  by  an  advancing  platform  within  the  brigade  (where  P destroyed 
=  1  given  that  the  bomb  and  platform  vehicle  meet  head-on),  or  the  bomb's 
battery  is  drained  of  voltage;  this  last  event  can  occur  anywhere  from  37.5  to 
52.5  min  (reflecting  a  preset  bomb  lifetime)  after  the  initial  bomb  is  inserted  into 
a  plant  cell.  Thus,  the  net  effect  of  this  first  scenario  phase  essentially  results  in  a 
preexisting  bomb  field  that  evolves  deterministically  once  it  moves  into  the  CA 
simulation  window. 


*  It  should  be  noted,  however,  that  the  potential  bomb  plant  sites  remain  fixed  relative  to  the 
brigade  node  cells. 
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The  second  phase  of  the  lO  scenario  commences  at  t„  =  48  min,  and  then 
continues  to  unfold  until  t„  =  96  min  (thus  overlapping  with  a  portion  of  the  first 
phase).  During  this  time,  jammer  bombs  are  inserted  into  the  light  gray  cells 
labeled  "aerial  barrage  targets,"  which  are  positioned  in  the  open  areas  between 
node  cells  (see  Figure  14).  The  probability  of  finding  a  jammer  bomb  at  a  cell, 
given  that  the  ceU  is  a  potential  aerial  barrage  target,  is  P(bomb  |  target  site) 
=  0.005/ time  step,  and  the  average  air-dropped  bomb  population  within  the  CA 
lattice  is  163  bombs  per  simulation  trial.  As  with  the  jammer  bombs  in  the  first 
phase  of  the  scenario,  bombs  move  forward  at  a  constant  rate  of  1  cell 
length/ time  step,*  and  they  continue  to  jam  friendly  nodes  xmtil  the  bomb  is 
either  destroyed  or  runs  out  of  battery-generated  power;  in  this  phase,  the  latter 
event  consistently  occurs  75  min  (reflecting  the  maximum  bomb  lifetime)  after 
the  bomb  is  inserted  into  a  target  cell.  Thus,  the  second  scenario  phase 
essentially  results  in  an  evolving  bomb  field  superimposed  over  a  portion  of  the 
preplanted  bomb  field.  Note  that  even  though  the  bombs  are  stochastically 
introduced  into  the  CA  simulation  window,  they  behave  deterministically  once 
introduced. 

Twenty  simulation  trials  of  the  two-phase  jammer  bomb  scenario  are  rim  using 
the  CA  model,  where  each  trial  is  seeded  with  a  different  random  number 
(needed  to  stochastically  vary  the  insertion  of  jammer  bombs  into  botii  plant  and 
target  cells).  Figure  15  displays  six  different  "snapshots"  taken  at  various  time 
steps  throughout  one  of  the  simulation  trials.  The  color  of  a  platform  node  cell 
indicates  one  of  five  fitness  levels  (as  defined  in  Figure  7),  and  a  black  ceU 
indicates  tiie  presence  of  a  jammer  bomb.  The  time  reported  in  the  black  text  box 
at  the  bottom  of  a  snapshot  reflects  a  dimensionless  simulation  time  step,  rather 
than  the  scaled  time  in  units  of  minutes.  In  these  snapshots,  the  following 
progression  of  events  unfolds: 

•  Time  =  3  min.  The  moving  brigade  approaches  a  preplanted  field  of 
jammer  bombs. 

•  Time  =  14  min.  All  of  the  platforms  in  the  lead  security  force  are  jammed  as 
this  unit  penetrates  the  bomb  field,  resulting  in  very  low  fitness.  Even 
though  they  are  not  directly  jammed,  many  of  the  platforms  in  the  front 
half  of  the  brigade  are  reduced  from  very  high  to  high  fitness  as  a  result  of 
jamming  in  the  lead  unit. 

•  Time  =  43  min.  The  brigade  continues  to  penetrate  the  bomb  field.  The 
lead  security  force  prevents  most  of  the  bombs  from  reaching  the  center 
platforms  within  the  brigade;  however,  fitness  levels  of  these  inner 
platforms  are  still  reduced  because  of  jamming  in  the  front  and  sides 
of  the  brigade. 


*  As  with  the  potential  bomb  plant  sites,  the  aerial  barrage  target  sites  also  remain  fixed  relative 
to  brigade  node  cells. 
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Figure  15.  CA  jammer  bomb  simulation  snapshots:  (e)  time  =  98  min;  (f)  time  =  131  min 
(continued). 

•  Time  =  68  min.  As  the  lead  security  force  reaches  the  end  of  the  preplanted 
bomb  field,  an  aerial  barrage  has  recently  commenced,  resulting  in  an 
increased  number  of  jammer  bombs  within  the  brigade.  Many  platforms 
operate  under  reduced  fitness  levels  (high,  medium,  and  very  low)  because 
they  are  either  directly  jammed  or  they  lose  commrmication  channel 
continuity  with  neighboring  platforms. 

•  Time  =  98  min.  Although  most  of  the  preplanted  jammer  bombs  have 
expired  or  been  destroyed  (reverting  most  of  the  platforms  in  the  front  half 
of  the  brigade  back  to  their  original  very  high  fitness  states),  the  air-dropped 
bombs  continue  to  degrade  fitness  levels  of  platforms  in  the  back  half  of  the 
brigade. 

•  Time  =  131  min.  Most  of  the  brigade  has  moved  beyond  the  air-dropped 
jammer  bombs;  fitness  levels  for  these  platforms  have  reverted  back  to  very 
high,  leaving  only  a  few  jammed  or  reduced-fitness  platforms  in  the  rear. 

These  snapshots  were  selected  from  a  total  of  225  time  samples.  They  were 
generated  within  a  single  simulation  run,  where  the  uniform  sampling  rate  is 
1  time  sample/45  s. 
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Figure  16  displays  a  plot  of  the  time-dependent  fraction  of  jammed  network 
nodes,  resultant  from  the  jammer  bomb  lO  scenario.  This  scenario  was 

run  over  multiple  instances  (where  an  instance  refers  to  a  single-trial  stochastic 
initial  placement  of  jammer  bombs  within  the  2-D  CA  lattice)  using  the  set  of 
20  sequential  random  seeds  {0, 1,  2,...,  19}.*  The  time  series  in  this  plot  rises  and 
then  falls  in  a  quasi-linear  fashion  (except  for  an  exponential-like  tail),  with 
numerous  stair-step  regions  indicating  short-term  equilibrium  jamming 
conditions.  In  this  time  series,  {/?’''■“*(!«))  seems  to  stabilize  for  finite  intervals  of 
time  until  advancing/retreating  jammer  bombs  force  the  fraction  of  jammed 
nodes  to  adjust  accordingly.  Finally,  the  maximal  value  of  (/z®‘''“*(f«))  (0.21,  which 
occurs  at  about  68  min  into  the  scenario)  roughly  corresponds  to  the  point  of  the 
maximal  jammer  bomb  population  within  the  CA  brigade.  This  is  also  the  point 
where  preplanted  bombs  have  completed  their  entry  into  the  brigade  from  the 
front;  thus,  all  new  incoming  bombs  are  stochastically  introduced  via  the  aerial 
barrage  phase  of  the  scenario. 

Fraction  of  Jammed  Modes  Uithin  the  Brigade  --  20  Run  flug 


Figure  16.  Fraction  of  jammed  nodes  as  a  function  of  scenario  time  t„,  for  the  CA  brigade 
in  the  jammer  bomb  scenario. 

Figure  17  is  a  plot  of  the  d}mamic  globally  averaged  SA  connectivity  for  the  CA 
brigade,  (C(fn)),  which  is  averaged  over  20  simulation  trials  using  the  standard 
set  of  random  seeds,  where 

7  612 

=  (21) 


*  Since  this  sequence  of  20  random  seeds  will  also  be  used  in  subsequent  jammer  bomb 
simulation  runs  for  consistency  purposes,  it  will  henceforth  be  referred  to  as  the  standard  set  of 
random  seeds. 
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SA  Connectivity  for  Jemiier  Bomb  Scenario  —  £0  Run  Avg 


Simulation  Time  <Minutes> 

Figure  17.  Dynamic  globally  averaged  SA  connectivity  as  a  function  of  scenario  time  f„, 
for  the  CA  brigade  in  the  jammer  bomb  scenario. 

for  each  trial.  Interestingly,  the  time  series  in  this  plot  approximately  resembles  a 
mirror-image  reflection  of  the  time  series  depicted  in  Figure  16,  where 

the  minimal  value  of  {C{t„))  (0.633)  occurs  concurrently  with  the  maximal  value 
of  As  with  the  time  series,  the  SA  connectivity  stabilizes  for 

finite  intervals  of  time  tmtil  advancing/retreating  local  jamming  environments 
force  (C(f«))  to  adjust  accordingly. 

Figure  18  shows  plots  of  Pv&tyJiigh(^^n)/  Pmediumitn')/  Plowi^^n)/  and  P very  low  (tn)  as  a 

function  of  scenario  time  f„,  where  again  Pfuness  state  (U  is  the  fraction  of  nodes  in  a 
fitness  state  at  time  tn-  Each  of  the  time  series  plots  show  how  platform  nodes 
simultaneously  transition  from  the  unperturbed  very  high  fitness  state  down  to 
either  the  high  or  very  low  fitness  states  for  most  nodes  (with  the  remaining  nodes 
transitioning  to  either  the  medium  or  low  states).  Then,  there  is  an  exponential 
return  to  the  original  very  high  state  as  the  brigade  moves  beyond  the  jammer 
bombs.  In  particular,  the  Pvery  low  (tn)  time  series,  which  rises  and  then  falls  in  a 
fairly  linear  fashion  (except  for  the  exponential  tail),  accounts  for  those  nodes 
whose  SA  connectivity  is  mainly  impacted  by  local  jamming.  On  the  other  hand, 
the  Phigh(tn)  time  series  accoimts  for  nodes  which  are  mainly  impacted  by  the  loss 
of  communication  channel  continuity  with  nonlocal  neighbor  nodes.  Finally,  for 
comparison  purposes.  Figure  19  displays  the  data  from  all  of  the  20  simulation 
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trials  for  both  the  Pvenj  high{tn)  (green  points)  and  Pvery  low  {tn)  (red  points)  time  series, 
whose  average  values  were  plotted  in  Figure  18.  Thus,  in  Figure  19,  each  of  the 
green  and  red  points  represent,  for  one  simulation  instance,  the  fraction  of  nodes 
in  the  very  high  and  very  low  fitness  states,  respectively,  measured  at  a  discrete 
simulation  time  step. 

Figures  20  and  21  illustrate  the  coarse-grained  state  and  temporal  entropies 

H{t„)  and  H{t„,  f„-i),  respectively,  for  the  jammer  bomb  scenario.  As  introduced  in 
section  2.6.4,  "entropy"  is  a  measure  of  the  disorder/volatility  within  a  multi¬ 
element  system  which  can  range  from  0  up  to  1.  In  both  figures,  the  entropy  is 
sampled  according  to  three  different  node  categories. 

•  Global.  The  sample  population  consists  of  all  nodes  within  the  brigade 
structure. 

•  CSMA  net.  The  independent  sample  populations  are  separate  battalion- 
level  CSMA  net  communities  (including  the  brigade-area  CSMA  net);  after 
calculation,  the  net-specific  entropies  are  averaged  to  come  up  with  a  single 
value. 

•  Local  net.  The  independent  sample  populations  are  separate  local 
SINCGARS  net  communities  (including  individual  EPLRS-equipped  self¬ 
servers,  which  are  treated  as  local  nets  consisting  of  one  member);  after 
calculation,  the  net-specific  entropies  are  again  averaged  to  come  up  with  a 
single  value. 

Both  state  and  temporal  entropies  are  sampled  as  a  function  of  discrete 
simulation  time  step  and  averaged  over  20  simulation  nms  using  the  standard 
set  of  random  seeds. 

The  state  entropy  time  series  depicted  in  Figure  20  illustrate  the  progressive 
variation  in  successive  nodal  SA  fitness  levels  as  distributed  among  the  five 
fitness  states  for  global,  CSMA  net,  and  local  net  populations.  In  this  context, 
H{t„)  =  0  indicates  that  aU  nodes  within  a  sampling  population  occupy  the  same 
fitness  state,  while  H{t„)  =  1  indicates  that  the  nodes  are  equally  distributed 
among  the  five  possible  fitness  states.  For  the  time  series  plots  shown  in  the 
figure,  the  state  entropy  of  the  global  population  first  jumps  to  a  quasi-stable 
level  between  0.4  and  0.5  tmtil  f„  =  78-80  min,  at  which  point  tihe  state  entropy 
slightly  decreases  to  a  level  between  0.35  and  0.4,  and  finally  decreases  to  a  level 
of  0  once  all  node  vehicles  have  moved  sufficiently  beyond  the  last  jammer 
bombs.  On  the  other  hand,  the  time  series  reflecting  averaged  state  entropy  for 
CSMA  and  local  net  populations  both  steadily  reach  peak  levels  of  0.35  and  0.28 
at  the  transition  time  step  t„  =  78-80  min,  at  which  point  they  steadily  decrease 
down  to  a  level  of  0.  Note  that  although  both  time  series  share  very  similar  time 
profiles,  the  local  net  state  entropy  is  always  slightly  less  than  tiiat  for  the  CSMA 
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Temporal  Entropy 


nets.  Taken  together,  the  three  time  series  describe  a  situation  where  SA  fitness 
levels  distributed  across  the  nodes  within  a  local  SINCGARS  net  (and,  to  a  lesser 
extent,  across  the  nodes  within  a  CSMA  net  community)  are  roughly  coherent  at 
each  point  in  time,  while  SA  levels  distributed  across  the  entire  brigade  structure 
are  moderately  varied. 

The  temporal  entropy  time  series  depicted  in  Figure  21  illustrates  the  persistence 
(or  lack  of)  of  nodal  fitness  levels  across  single  time  step  transitions  throughout 
the  scenario  run  time.  In  this  context,  H{tn,  fn-i)  =  0  indicates  that  all  nodes  within 
a  sampling  population  remain  in  the  same  fitness  state  from  t„-i  to  t„,  while  H{tn, 
tn-i)  =  1  indicates  that  the  nodes  are  equally  likely  to  transition  to  any  of  the  five 
fitness  states  over  this  single-step  time  interval.  As  with  the  state  entropy  time 
series,  the  temporal  entropy  time  series  for  each  sampling  population  evolve  to 
maximum  levels,  and  then  subsequently  decrease  back  to  a  level  of  0.  However, 
in  this  case,  the  time  series  for  the  different  sample  populations  are  fairly  spread 
out,  with  global,  CSMA  net,  and  local  net  temporal  entropies  achieving  maximal 
levels  of  0.38,  0.23,  and  0.06,  respectively.  Taken  together,  the  three  temporal 
entropy  time  series  illustrate  an  increasing  trend  in  SA  fitness  level  variability  as 
a  function  of  increasing  sampling  population  size,  from  very  stable  SA  levels 
across  local  nets  up  to  moderately  rmstable  SA  fitness  levels  across  the  global 
brigade  nodal  population. 

In  the  case  of  SA  connectivity,  stress-induced  perturbations  act  to  drive  the 
global  network  state  vector  to  various  points  within  the  SA  connectivity  phase 
space  (as  previously  introduced  in  section  2.6.3).  However,  it  is  very  difficult  to 
explore  this  high-dimensional  phase  space  for  attractors;  instead,  we  look  for  a 
collective  structure  within  the  coarse-grained  fitness  space  by  studying  the  various 
Pfitness{tn)  time  series  that  are  generated  during  a  simulation.  This  structure  is  said 
to  demonstrate  non-trivial  collective  behavior  (NTCB),  which  was  first  reported  by 
Chat6  and  Manneville  as  observed  in  the  global  dynamics  of  four-dimensional 
cellular  automata  (Chat6  and  Marmeville  1991;  Chat6  and  Manneville  1992). 

Within  the  context  of  a  return  map,  collective  structure  can  be  defined  as  a  set  of 
one  or  more  points  representing  stationary  Pfitness  levels  towards  which 
perturbative  stress  sources  collectively  drive  the  entire  population  of  networked 
nodes.  Since  an  lO  threat  scenario  can  posit  an  open  system,  where  perturbative 
energy  is  stochastically  introduced  into  a  simulation  from  locations  outside  of  the 
CA  lattice,  the  interaction  dynamics  between  network  nodes  and  external  stress 
sources  can  at  best  be  deterministic  and  discontinuous.  It  should  thus  be  noted 
that  the  discontinuous  nature  of  the  Pfitness{tn)  time  series  will  result  in  collective 
structures  which  can  be  "smeared"  across  variable-sized  regions  of  fitness  state 
space. 

Figure  22  displays  a  return  map  associated  with  tire  jammer  bomb  scenario;  this 
map  was  created  with  the  Pvery  high{tn)  (green  points)  and  Pvenj  low  (tn)  (red  points) 


36 


Jdfftftier  Bomb  Return  Hap  —  Fitness  ■  Very  High  8.  Very  Lou 
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Figure  22.  Return  map  for  the  Pveryhigh{tn)  (green  points)  and  P^ery  low  {tn)  (red  points)  time 
series. 

data  previously  shown  in  Figure  19.  This  return  map  illustrates  the  time  series 
dynamics  in  a  recursive  rather  than  successive  format.  While  the  red  points 
(clustered  along  a  diagonal  line  with  slope  =  1)  reinforce  the  linear  behavior  of 
the  Pvery  low  (tn)  time  series,  the  green  points  display  the  complex  emergent 
behavior  of  the  Pvay  high{tn)  time  series  in  what  is  clearly  a  collective  structure. 
There  are  indications  of  what  appears  to  be  multimodal  dynamics  (symmetric 
branching  along  the  diagonal  line*)  within  certain  intervals  of  Pvery  highitn)- 

Although  the  collective  structure  associated  with  the  Pvay  high{tn)  time  series 
resembles  a  period  doubling  bifurcation  arising  within  a  1-D  chaotic  map  (Ott 
1993),  the  jammer  bomb  scenario  clearly  defines  an  open  system  (illustrated  by 
the  dispersive  "smearing"  of  Pvery  high(tn)  data  points).  Thus,  the  collective 
phenomena  depicted  in  the  figure  cannot  strictly  qualify  as  deterministic  chaos 
because  the  global  system  state  at  to  (i.e.,  initial  positions  of  preplanted  jammer 
bombs  relative  to  the  brigade)  does  not  alone  guarantee  a  determinable  global 
SA  cormectivity  state  vector  for  tn  >  to. 

3.2.3  Echelon-Specific  Nodal  Sampling 

By  adjusting  the  scope  of  the  node  sampling  populations  within  the  CA  model,  it 
is  possible  to  separately  measure  the  time-series  response  of  network  nodes  at 
the  brigade,  battalion,  and  combined  company /platoon  echelon  levels 


*  This  symmetric  return  map  structure  is  very  likely  resultant  from  the  vertical-axis  spatial 
symmetry  in  the  simple  movement  to  contact  brigade  structure  (Figure  3). 
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throughout  the  CA  brigade.  For  example.  Figures  23  and  24  display  the  dynamic 
fraction  of  jammed  nodes  and  the  globally  averaged  SA  onnectivity,  respectively, 
specifically  sampled  from  brigade,  battalion,  and  combined  company/platoon 
echelon  node  populations  using  the  standard  set  of  random  seeds.  In  Figure  23, 
the  combined  company/platoon  echelon  time  series  response  is  very  similar  to 
that  of  the  global  node  population  (see  Figure  16),  while  battalion  and  brigade 
echelon  time  series  display  lower  levels  and/ or  later  onset  of  significant  levels  of 
Qistress{^l^yf  This  latter  response  results  from  many  of  the  battalion  and  brigade 
nodes  being  positioned  apart  from  the  front  and  flanks  of  the  advancing  brigade 
structure.  Figure  24  depicts  similar  relative  behavior  amongst  the  node 
populations,  with  brigade  and  combined  company/  platoon  nodes  exhibiting 
very  similar  dynamic  levels  of  SA  connectivity  as  that  calculated  for  the  entire 
network  node  population  (see  Figure  17),  while  battalion  nodes  exhibit  slightly 
higher  levels  of  SA  comiectivity  throughout  the  jammer  bomb  scenario. 


Simulation  Time  <Mlnutes> 

Figure  23.  Fraction  of  jammed  nodes  as  a  function  of  scenario  time  f„,  for  brigade  (Bde), 
battalion  (Bn),  and  combined  company/ platoon  (Co-Pit)  echelon  nodes. 

3.2.4  Sensitivity  Investigation  of  Model  Parameters 

The  CA  model  is  constructed  such  that  various  sensitivity  analyses  can  be 
performed  on  parameters,  such  as  AO  neighborhood  size,  jammer  bomb  residual 
lifetime,  and  jamming  neighborhood  size.  In  this  section,  the  sensitivity  of  model 
parameters  is  demonstrated  by  varying  the  conditions  of  the  jammer  bomb  lO 
scenario  discussed  in  sections  3.2.1  and  3.2.2.  Note  that  all  simulation  results 
presented  here  are  generated  by  using  the  standard  set  of  random  seeds 
previously  introduced  in  section  3.2.2. 
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Sinulation  Time  <Minutes> 

Figure  24.  Globally  averaged  SA  connectivity  as  a  function  of  scenario  time  tn/  for 
brigade  (Bde),  battalion  (Bn),  and  combined  company/platoon  (Co-Pit) 
echelon  nodes. 

3.2.4. 1  Reduced  AO  Neighborhood  Scenario 

In  the  first  modified  scenario,  the  AO  neighborhoods  originally  associated  with 
brigade,  battalion,  company,  and  platoon  echelon  nodes  within  the  CA  brigade 
(Figure  4)  are  each  decreased  by  a  factor  of  0.5  in  both  x-  and  y-directions  (thus 
decreasing  total  neighborhood  areas  by  a  factor  of  0.25).  These  re-scaled  AO 
neighborhoods  are  illustrated  in  Figure  25.  As  a  result  of  these  reduced 
neighborhoods,  the  connectivity  distribution  P(N)  for  the  CA  brigade  has  also 
adjusted  accordingly,  as  illustrated  in  Figure  26.  It  is  interesting  to  note  that  this 
new  connectivity  distribution  approximately  resembles  a  power-law  distribution 
(except  for  the  large  number  of  nodes  with  N  =  495)  characteristic  of  scale-free 
networks.  In  these  networks,  the  coimectivity  distribution  is  described  by  P(N) 
~  N'“,  where  2.1  <  a  <  4.0  (Barabasi  and  Albert  1999). 

Figure  27  displays  a  plot  comparing  the  time-dependent,  globally  averaged  SA 
connectivity  resultant  from  both  the  original  jammer  bomb  and  reduced  AO 
neighborhood  scenarios  (a  comparative  plot  of  (/z®'’'®®®(fM))  vs.  tn  is  not  presented 
since  jammer  bomb  parameters  remain  imaltered  from  the  original  scenario). 
Asis  clearly  demonstrated  in  the  figure,  uniformly  reducing  all  AO 
neighborhood  areas  by  a  factor  of  0.25  only  serves  to  increase  average  SA 
connectivity  levels  by  10%  or  less.  This  result  demonstrates  the  nonlinear 
response  of  overall  SA  network  topology  to  changes  in  echelon-specific 
connectivity  ranges. 
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Figure  25.  Reduced  AO  neighborhood  scaling  within  the  CA  brigade. 


Reduced  AO  Neighborhoods 


Figure  26.  Connectivity  distribution  P(N)  for  the  CA  brigade  with  reduced  AO 
neighborhoods. 


SA  Connectivity  Hith  Reduced  AO  Neighborhoods  —  20  Run  Avg 


Figure  27.  Globally  averaged  SA  connectivity  as  a  function  of  scenario  time  for  both 
the  original  jammer  bomb  and  reduced  AO  neighborhood  scenarios. 

Figure  28  displays  time  series  plots  of  the  fitness  metrics  Pvay  high(tn),  Phighitn), 
Pmediumitn),  Pioiv(t„),  and  Pvery  low  (tn)  from  the  reduced  AO  neighborhood  scenario. 
As  with  the  SA  connectivity  time  series  depicted  in  Figure  27,  the  Pvery  higk(tn)  and 
Phigh(f«)  time  series  display  a  slight  increase  and  decrease  (respectively)  in 
dynamic  levels  by  factors  between  0.05  and  0.10  when  compared  with  similar 
measures  from  the  original  jammer  bomb  scenario  (see  Figure  18).  The 
remaining  three  time  series,  however,  essentially  remain  unaltered  from  the 
previous  scenario  results.  This  result  demonstrates  that  the  most  significant 
response  to  modifying  SA  network  connectivity  ranges  manifests  in  nonlocal 
interactions  (i.e.,  those  nodes  which  are  not  directly  jammed  but  nevertheless 
suffer  connectivity  losses  due  to  jammed  neighbor  nodes). 

Figures  29  and  30  display  state  and  temporal  entropies,  respectively,  as  a 
function  of  time  step  f«  for  the  reduced  AO  neighborhood  scenario.  In  the  case  of 
measurements  made  with  respect  to  the  global  nodal  population,  dynamic  levels 
of  both  state  and  temporal  entropies  have  significantly  increased  from  those 
resultant  from  the  original  jammer  bomb  scenario,  especially  during  the  first 
60  min  of  the  scenario  time  window.  This  is  likely  due  to  a  reduction  in 
correlation  between  intemodal  SA  connectivity  levels,  which  results  from 
decreased  overlapping  of  AO  neighborhoods  associated  with  neighboring  nodes. 
State  and  temporal  entropies  measured  from  CSMA  and  local  net  populations, 
however,  essentially  remain  unchanged  from  the  previous  scenario. 
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Temporal  Entropy  lor  Reduced  AO  Neighborhoods  —  20  Run  Avg 


Figure  30.  Temporal  entropy  as  a  function  of  scenario  time  for  the  reduced  AO 
neighborhood  scenario. 

Figures  31  and  32  display  the  data  from  all  of  the  20  simulation  trials  of  the 
reduced  AO  neighborhood  scenario  for  both  Pvery  highitn)  (green  points)  and  P very  low 
(tn)  (red  points)  in  standard  time  series  and  return  map  formats,  respectively.  In 
Figure  31,  each  of  the  time  series  resembles  its  counterpart  from  the  original 
jammer  bomb  scenario  (Figure  19),  except  that  the  branching  observed  in  the  Pvery 
high{tn)  time  series  from  the  original  scenario  has  disappeared.  Similarly,  the 
return  map  associated  with  the  Pvery  high(tn)  data  depicted  in  Figure  32  has  lost  the 
definitive  multimodal  branching  structure  observed  in  the  Pvery  highitn)  return  map 
from  the  original  scenario  (Figure  22).  These  results  demonstrate  that  the  overall 
dynamic  response  of  the  stressed  global  networked  system  becomes  more  linear 
as  the  number  of  internodal  connections  is  reduced. 

3.2.4.2  Reduced  Jammer  Bomb  Lifetime  Scenario 

In  the  second  modified  scenario,  the  jammer  bomb  scenario  is  next  modified  as 
follows; 

•  During  the  first  phase  of  the  scenario  (i.e.,  preplanted  bomb  field),  jammer 
bomb  functional  lifetimes  are  reduced  by  a  factor  of  0.5,  now  ranging  from 
18.75  to  26.25  min  after  a  bomb  is  initially  inserted  into  a  plant  cell. 

•  During  the  second  phase  of  the  scenario  (i.e.,  aerial  barrage),  jammer  bomb 
functional  lifetimes  are  also  reduced  by  a  factor  of  0.5,  so  that  bombs  now 
uniformly  expire  37.5  min  after  being  inserted  into  a  target  cell. 

Note  that  the  AO  neighborhoods  for  all  network  nodes  are  restored  to  the 
original  configuration,  as  depicted  in  Figure  4. 
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P<Very  High>CGREEM]  or  P<Very  LowJCRE 


Figure  32.  Return  map  for  the  Pveryhighitn)  (green  points)  and  Pven/iou>(tn)  (red  points)  time 
series  from  the  reduced  AO  neighborhood  scenario. 


Figure  33  displays  a  plot  of  the  time-dependent  fraction  of  jammed  network 
nodes  resultant  from  the  reduced  bomb  lifetime  scenario.  The  impact  of  shorter- 
lived  jammer  bombs  is  clearly  illustrated  in  this  figure,  which  compares 
resultant  from  the  current  scenario  with  the  time  series  response  from  the 
original  scenario.  Here,  it  is  seen  that  the  response  is  identical  out  until  t„ 
=  18.75  min,  the  point  where  preplanted  bombs  begin  to  expire.  From  that  point 
in  time  until  about  t„  =  100  min,  the  difference  between  the  two  time  series 
ranges  from  about  25%  to  50%,  where  reduced  bomb  lifetime  significantly 
reduces  Of  particular  interest  is  an  interval  of  apparently  quasi¬ 

equilibrium  behavior  within  the  modified  scenario  time  series  (where 
seems  to  temporarily  stabilize),  commencing  at  about  tn  =  20  min  (right  after 
preplanted  bombs  begin  to  expire)  and  continuing  until  =  48  min  (the 
commencement  of  the  aerial  barrage  phase).  Finally,  the  improved  network 
performance  is  less  significant  from  t„  =  100  min  to  the  end  of  the  scenario,  at 
which  point  most  of  the  preplanted  bombs  in  the  original  scenario  have  expired 
and  the  brigade  gradually  moves  beyond  the  remaining  bombs  delivered  during 
the  recently  concluded  aerial  barrage. 


Traction  of  Jammed  Nodes  Within  the  Brigade  —  20  Run  Aug 


Figure  33.  Fraction  of  jammed  nodes  as  a  function  of  scenario  time  tn,  for  both  the 
original  and  reduced  jammer  bomb  lifetime  scenarios. 

Figures  34-37  depict  time  series  of  the  globally  averaged  SA  connectivity,  five- 
state  Pfiiness{tn)  mctrics,  state  entropy,  and  temporal  entropy,  respectively, 
associated  with  the  reduced  bomb  lifetime  scenario.  Each  of  these  time  series 
behaves  very  similarly  to  its  counterpart  from  the  original  scenario,  with  the 
notable  ubiquitous  exception  of  the  interval  of  quasi-equilibrium  behavior  first 
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Jammer  Bomb  SA  Connectivity  Hz’  Shortened  Bomb  Lifetime  —  20  Run  Aug 


Figure  34.  Globally  averaged  SA  connectivity  as  a  function  of  scenario  time  tn,  for  the 
reduced  jammer  bomb  lifetime  scenario. 


Simulation  Time  (minutes) 


FigUfG  35.  Pvery  high{tn^/  Phigh{tn)/  Pmedhimi^tr^/  Plowi^iji)/  Sfld,  P very  low  (fn)  ^.S  cl  funCtiOH  of  SCGnS-fiO 

time  tn,  for  the  reduced  jammer  bomb  lifetime  scenario. 


state  Entropy  for  Shortened  Bomb  Lifetime  —  20  Run  Aug 


Figure  36.  State  entropy  as  a  function  of  scenario  time  for  the  reduced  jammer  bomb 
lifetime  scenario. 


Temporal  Entropy  for  Shortened  Bomb  Lifetime  —  20  Run  Aug 


Figure  37.  Temporal  entropy  as  a  function  of  scenario  time  fn,  for  the  reduced  jammer 
bomb  lifetime  scenario. 
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pointed  out  in  the  time  series  from  tihe  modified  scenario  (Figure  33). 

In  particular,  the  Pvenj  high{tn),  Phigh{tn),  and  Pvery  i<m>  (tn)  time  series  in  Figure  35,  the 
CSMA  and  local  net  state  entropy  time  series  in  Figure  36,  and  the  local  net 
temporal  entropy  time  series  in  Figure  37  are  all  remarkably  level  tihroughout 
this  time  interval,  with  the  two  state  entropy  time  series  essentially  overlapping 
(indicating  near-identical  distributions  of  SA  connectivity  levels  throughout  both 
CSMA  and  local  net  nodal  populations).  Taken  togetiher,  these  plots  reinforce 
the  conclusion  that  SA  cormectivity  levels  have  achieved  a  near-equilibrium 
condition  by  =  20  min  (especially  within  battalion-level  and  local  SINCGARS 
net  communities),  which  is  perpetuated  until  the  commencement  of  the  aerial 
barrage  phase  of  the  jammer  bomb  scenario  at  tn  =  48  min. 

Figures  38  and  39  display  all  data  from  the  20  simulation  trials  of  the  reduced 
jammer  bomb  lifetime  scenario  for  both  Pvery  high{t„)  (green  points)  and  Pvery  low  (tn) 
(red  points)  in  standard  time  series  and  return  map  formats,  respectively.  Again, 
each  of  the  time  series  in  Figure  38  resembles  its  coxmterpart  from  the  original 
jammer  bomb  scenario,  except  that  the  branching  observed  in  the  Pvery  high{tn)  time 
series  from  the  original  scenario  has  now  sharply  differentiated  into  two 
relatively  level  branches  (with  some  additional  distributed  points  resulting  from 
variable  bomb  lifetimes)  within  the  tn  =  20  min  to  48  min  time  interval  discussed 
previously.  Similarly,  the  return  map  associated  with  the  Pvery  high{tn)  data 
depicted  in  Figure  39  retains  tine  definitive  multimodal  branching  structure  first 
observed  in  the  Pvery  high{tn)  return  map  from  the  original  scenario.  In  addition,  the 
Pvery  high{tn)  time  scries  in  Figure  38  exhibits  more  significant  point  dispersion  from 
about  tn  =  70  min  to  =  110  min  tiian  is  seen  in  the  original  scenario  response,  as 
well  as  a  more  abbreviated  tail  from  f„  =  110  min  to  the  end  of  the  scenario.  This 
behavior  is  illustrated  from  another  perspective  in  the  return  map  (Figure  39), 
where  the  lower  branching  structure  appears  to  exhibit  a  wider  spread  than  that 
observed  in  the  return  map  associated  with  the  original  scenario.  These  results 
seem  to  indicate  that  the  overall  dynamic  network  response  within  the 
deterministic  first  phase  of  the  jammer  bomb  scenario  becomes  more  stable  as 
jammer  bomb  lifetimes  are  reduced;  this  equilibrium  then  vanishes  as  the  second 
phase  commences  (which  is  a  characteristic  of  the  discontinuous  dynamics  of  this 
phase).  Due  to  the  emergent  nature  of  this  behavior,  however,  it  is  not 
appropriate  to  conclude  that  the  first  phase  equilibrium  would  continue  to 
manifest  itself  given  even  shorter  bomb  lifetimes  without  rxmning  further 
simulations. 

3.2.4.3  Combined  Moore  Neighborhood/Increased  Bomb  Density  Scenario 

In  the  third  modified  scenario,  tiie  original  jammer  bomb  scenario  presented  in 
sections  3.2.1  and  3.2.2  is  modified  as  follows; 

•  The  jamming  neighborhood  of  all  jammer  bombs  is  reduced  from  36  down 
to  8  neighboring  cells. 
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Shortened  Somb  Lifetine  —  Fitness  =  Very  High  and  Very  Lou  --  20  Run  Aggregate 
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lifetime  scenario. 
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Figure  39.  Return  map  for  the  Pvery  high(tn)  (green  points)  and  Pvery  low  (tn)  (red  points)  time 
series  from  the  reduced  jammer  bomb  lifetime  scenario. 
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•  The  average  density  of  preplanted  jammer  bombs  is  increased  by  a  factor  of 
five  (i.e.,  from  480  to  2,400  bombs  per  simulation  trial). 

•  The  average  density  of  air-dropped  jammer  bombs  is  increased  by  a  factor 
of  10  (i.e.,  from  163  to  1,630  bombs  per  simulation  trial). 

Figure  40  depicts  the  modified  jamming  neighborhood  as  defined  in  the  previous 
paragraph.  As  is  clearly  illustrated  in  the  figure,  the  modified  CA  neighborhood 
includes  only  those  neighboring  cells  that  share  either  an  edge  or  vertex  with  the 
center  cell.  This  type  of  nine-cell  configuration  is  conventionally  knovm  as  a 
Moore  neighborhood  (see  Figure  2[b]). 


Figure  40.  Reduced  jamming  neighborhood  consisting  of  eight  neighbor  cells. 

Figures  41  and  42  display  plots  of  the  time-dependent  fraction  of  jammed 
network  nodes  and  globally  averaged  SA  connectivity  resultant  from  the 
modified  lO  scenario,  respectively.  This  modified  scenario  was  run  using  the 
standard  set  of  random  seeds.  Although  the  average  net  jamming  area  coverage 
per  simulation  trial  has  increased  from  the  original  to  modified  scenarios  by 
factors  of  1.2  and  2.4  for  preplanted  and  air-dropped  bombs,  respectively,  the 
associated  changes  in  the  fraction  of  jammed  nodes  and  average  SA  connectivity 
are  clearly  nonlinear.  Figure  41  illustrates  this  point  by  comparing  the  time- 
series  response  of  the  fraction  of  jammed  nodes  for  both  the  original  and 
modified  jammer  bomb  scenarios.  The  high-frequency  structure  which  emerges 
(in  both  plots)  from  tn  =  8  min  to  t„  =  68  min  is  due  to  a  combination  of  a 
maximally  allowable  preplanted  bomb  population  (where  every  one  of  the  dark 
gray  cells  originally  shown  in  Figure  14  is  occupied  by  a  jammer  bomb  during 
each  advancing  time  step  in  the  first  scenario  phase)  and  a  jamming  range  which 
doesn't  necessarily  jam  both  the  primary  and  alternate  servers  in  a  radio  net 
(which  was  the  case  in  the  original  scenario). 
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Traction  of  Jasined  Nodes  Uithin  the  Brigade  —  &0  Run  Avg 


Figure  41.  Comparison  of  fraction  of  jammed  nodes  as  a  function  of  time  tn  for  both  the 
original  jammer  bomb  and  combined  Moore  neighborhood/increased  bomb 
density  scenarios. 


Jammer  Bomb  SA  Connectiuity  U/  Incr,  Bomb  Dens.  &  Red.  Jan  Neigh.  —  £0  Run  Aug 


Figure  42.  Dynamic  globally  averaged  SA  connectivity  for  the  combined  Moore 
neighborhood /increased  bomb  density  scenario. 
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Figure  43  displays  time  series  plots  of  the  fitness  metrics  Pve>y  highitn),  Phigh{tn), 
Pmediumitn),  Piowitn),  and  Pven,  low  {tn)  from  the  Combined  Moore 
neighborhood/increased  bomb  density  scenario.  Here,  the  high-frequency 
structure  observed  in  Figures  41  and  42  reemerges  in  the  Pvety  highitn)  and  P very  low 
(tn)  time  series  out  to  t„  =  68  min.  Then,  it  abruptly  disappears  when  the  PhighiU) 
time  series  emerges,  and  from  this  time  onward,  the  dynamics  appear  to  follow  a 
trend  similar  to  that  exhibited  in  the  original  jammer  bomb  scenario  (see  Figure 
18).  The  Pvery  high{tn)  and  Phighitn)  time  series  exhibit  an  approximate  reflective 
symmetry,  while  the  Pvery  low  (tn)  time  series  decreases  amplitude  in  a  quasi-linear 
fashion.  Because  this  abrupt  transition  in  the  time  series  dynamics  is  not  evident 
in  the  original  scenario,  it  reinforces  the  fact  that  global  emergent  behavior  in  a 
distributed  interdependent  system  is  usually  very  sensitive  to  any  changes  in 
model  parameters. 


Jarnmer  Bomb  Sfl  Connectivity  14/  Incr.  Bomb  Dens.  &  Red.  Jam  iieigh.  —  £0  Run  Rug 


Figure  43.  PveryUghitn),  Phighitn),  Pnwdiumitn),  Piowitn),  and  Pvery  law  (tn)  aS  a  fuRCtion  of  time  tn, 
for  the  combined  Moore  neighborhood/increased  bomb  density  scenario. 

Figures  44  and  45  display  state  and  temporal  entropies,  respectively,  as  a 
fimction  of  time  step  tn  for  the  combined  Moore  neighborhood /increased  bomb 
density  scenario.  The  high-frequency  structure  characteristic  of  this  scenario 
again  manifests  itself  within  the  plots  in  both  figures  from  tn  =  8  min  to  tn 
=  68  min.  However,  there  is  also  an  abrupt  increase  in  amplitude  in  both  state 
and  temporal  entropy  time  series  plots  occurring  at  about  f«  =  55  min.  This 
phenomenon  is  very  likely  catalyzed  by  the  commencement  of  the  aerial  barrage 
phase  of  the  scenario,  where  the  steady-state  stress  supplied  by  the  preplanted 
jammer  bombs  transitions  (within  the  time  interval  tn  =  48  min  to  tn  =  55  min) 
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Figure  44.  State  entropy  as  a  function  of  time  tn,  for  the  combined  Moore 
neighborhood /increased  bomb  density  scenario. 


Figure  45.  Temporal  entropy  as  a  function  of  time  tn,  for  the  combined  Moore 
neighborhood /increased  bomb  density  scenario. 
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into  a  form  stochastic  in  nature.  This  is  especially  noticeable  in  Figure  45,  where 
the  temporal  entropies  jump  from  low-amplitude  (^0.0.5  for  brigade,  CSMA  net, 
and  local  net  sampling  populations)  periodic  time  series  to  the  higher  amplitude 
and  noisy  format  previously  observed  in  the  other  jammer  bomb  scenario 
variants. 

Figures  46  and  47  display  the  data  from  all  of  the  20  simulation  trials  for  both 
Pvery  highlit)  (green  points)  and  Pveryiow  (tn)  (red  points)  in  standard  time  series  and 
return  map  formats,  respectively.  In  Figure  46,  both  time  series  exhibit  a 
gradually  increasing  two-state  periodic  behavior  from  t„  =  8  min  out  to  about  tn 
=  40  min  within  the  scenario.  This  behavior  then  transitions  into  a  dynamics 
exhibiting  data  clustering  for  both  Pveryhigh{tn)  and  Pvenjimo  (tn)  (similar  to  what  was 
observed  in  the  original  scenario  and  plotted  in  Figure  19),  plus  two  additional 
new  clusters  for  Pvery  high{tn)  in  the  probability  range  from  0.3  to  0.5.  As  with  the 
return  maps  associated  with  the  previously  discussed  jammer  bomb  scenario 
variants,  the  return  map  shown  in  Figure  47  displays  symmetry  relative  to  the 
diagonal,  with  slope  =  1  for  both  Pvery  high(tn)  and  Pvery  low  (tn);  in  particular,  the 
symmetric  behavior  of  Pvery  high(tn)  indicates  that  this  particular  time  series  exhibits 
quasi-periodic  dynamics,  where  levels  jump  back  and  forth  between  the  upper 
and  middle  data  clusters.  Thus,  platform  nodes  principally  impacted  by  local 
jamming  (represented  by  the  Pvery  iow(tn)  time  series)  consistently  exhibit  very  low 
situational  awareness.  Other  nodes,  whose  situational  awareness  is  mainly 
impacted  by  communication  loss  with  nonlocal  neighbor  nodes  (represented  by 
the  Pvery  high(tn)  and  Phigh(tn)  time  series),  oscillate  between  very  high  and  high  levels 
of  situational  awareness. 

3.2.5  Jammer  Bomb  Susceptibility  Gradient 

Up  to  this  point,  it  has  been  assumed  that  the  jamming  neighborhood  of  a 
jammer  bomb  is  isotropic.  In  other  words,  a  node  platform  is  assumed  jammed  if 
the  Euclidian  distance  between  the  centers  of  the  fc*  jammer  bomb  and  i^  node 
cells  dik  <  3.2  cell  lengths.  This  assumption  is  predicated  upon  another 
simplifying  assumption:  that  each  node  platform  is  equally  susceptible  to 
jamming  from  the  notional  jammer  bomb  stress  source.  A  more  reahstic 
assumption  might  be  that  different  types  of  node  platforms  each  have  their  own 
specific  susceptibility  threshold  to  jamming,  which  is  based  on  the  incident  power 
received  by  a  node  radiating  from  a  jamming  source.  Thus,  the  original  jammer 
bomb  lO  scenario  is  modified  to  accoimt  for  variable  platform  susceptibiHty 
thresholds. 

In  this  modification  to  the  jammer  bomb  scenario,  each  bomb  is  equipped  with  a 
set  of  concentric,  embedded  jamming  neighborhoods,  as  depicted  in  Figure  48. 
Each  of  the  seven  neighborhoods  shown  in  this  figure  is  defined  by  a  unique  dik 
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Figure  46.  Superposition  of  Pvenj  high{tn)  (green  points)  and  Pveiy  imu  (fn)  (red  points)  time 
series  data  from  the  combined  Moore  neighborhood /increased  bomb  density 
scenario. 
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Figure  47.  Return  map  for  the  Pvayhigh(tn)  (green  points)  and  Pveryiow  (tn)  (red  points)  time 
series  from  the  combined  Moore  neighborhood/increased  bomb  density 
scenario. 
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Figure  48.  The  seven  concentric  jamming  neighborhoods  used  to  define  the  discrete 
susceptibility  gradient. 

(where  the  node  can  reside  in  any  of  the  36  lattice  cells  surrounding  the  fc* 
jammer  bomb  cell  at  the  center)  indicated  by  the  numbers  written  within  each 
cell,  where  d*  =  1.0,  1.4,  2.0,  2.2,  2.8,  3.0,  3.2.  In  addition,  each  node  within  the 
brigade  is  assigned  (based  on  sampling  from  a  uniform  random  distribution)  one 
of  seven  possible  jamming  susceptibility  thresholds  Si"’"'”*”?,  where  Si'""'"'”?  =  dik 
(thus  associating  a  susceptibility  threshold  with  a  specific  jamming 
neighborhood).  Once  assigned  a  susceptibility  threshold,  each  node  then  also 
has  an  associated  jamming  resistance  or  "hardness"  threshold  where 

^{jamming  =  ijgjammmg  0.3125  <  ^  1.000.  Then,  a  node  is  jammed  if  and 

only  if  ^  dik  (or  equivalently,  ^  dik).  The  collection  of  concentric 

jamming  neighborhoods,  taken  in  conjimction  with  the  set  of  threshold  values 
for  nodal  susceptibility  5/'”'””'%  defines  a  discrete  susceptibility  gradient.  It  should 
be  noted  that  this  gradient  is  a  notional  construct  that  qualitatively  represents  a 
form  of  nodal  jamming  susceptibility  which  monotonically  decreases  as  a 
function  of  increasing  dik  (where  the  received  jammer  power  oc  l/d“^ ). 

Figures  49-51  display  the  fraction  of  jammed  nodes,  globally  averaged  SA 
connectivity,  and  five-state  fitness  time  series,  respectively,  averaged  across 
multiple  simulation  rxms  (using  the  standard  set  of  random  seeds).  The  results 
depicted  in  these  three  plots  are  very  similar  in  amplitude  to  those  associated 
with  the  combined  Moore  neighborhood/increased  bomb  density  scenario 
(see  Figures  41-43  in  section  3.2.4.3).  The  most  obvious  difference  between  the 
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Figure  49.  Comparison  of  fraction  of  jammed  nodes  as  a  function  of  time  within  the 
susceptibility  gradient  and  original  jammer  bomb  scenarios. 
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Figure  50.  Dynamic  globally  averaged  SA  connectivity  for  the  susceptibility  gradient 
scenario. 
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Figure  51.  Pveryhigh{tn),  Phigh{tn),  Pmedmm{tn),  Piow{tn),  and  Pvery  iow{t„)  as  a  function  of  time  t„  for 
the  susceptibility  gradient  scenario. 

the  two  sets  of  time  series  plots  is  the  lack  in  the  susceptibility  gradient  plots  of 
the  periodic  high-frequency  structure  observed  in  the  first  half  of  the  previous 
scenario  plots.  However,  from  about  t„  =  68  min  to  the  termination  of  the 
scenario,  plots  from  both  scenarios  demonstrate  similar  aperiodic  low-amplitude 
structure.  In  the  case  of  the  current  scenario,  this  aperiodic  structure  is  likely  due 
to  the  variation  in  susceptibility  threshold  levels  amongst  the  jammed  nodes, 
where  locally  clustered  nodes  within  a  jammer  bomb  neighborhood  can  no 
longer  be  assumed  to  be  concurrently  jammed. 

Figures  52  and  53  display  the  state  and  temporal  entropy  time  series, 
respectively,  for  the  susceptibility  gradient  scenario,  again  averaged  across 
simulation  runs  using  the  standard  set  of  random  seeds.  As  with  similar  plots 
associated  with  the  original,  reduced  AO  neighborhood,  and  reduced  jammer 
bomb  lifetime  scenarios,  both  state  and  temporal  entropies  associated  with  the 
current  scenario  demonstrate  dintinct  separation  between  brigade,  CSMA  net 
and  local  net  sampling  populations.  Entropy  amplitude  levels  are  also 
comparable  to  those  observed  in  each  of  the  previous  scenarios  (including  the 
combined  Moore  neighborhood/increased  bomb  density  scenario).  One, 
difference,  however,  is  the  tendency  for  the  susceptibility  gradient  scenario 
entropies  (especially  those  associated  with  the  global  brigade  sampling 
population)  to  suddenly  drop  significantly  in  magnitude  at  about  f„  =  68  to 
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70  min.  This  sharp  dropoff  is  clearly  associated  with  the  previously  noted  point 
in  time  within  each  jammer  bomb  scenario  variant  where  maximal  bomb  density 
within  the  brigade  structure  is  reached  and  then  is  persistently  reduced, 
indicating  a  slight  increase  in  ordered  behavior  as  bombs  switch  from  entering  to 
exiting  the  brigade  territory. 

Figures  54  and  55  display  all  susceptibility  gradient  simulation  trial  data  for  both 
Pvery  high{tn)  (green  points)  and  Pvery  low  (tn)  (red  points)  in  standard  time  series  and 
return  map  formats,  respectively.  There  is  some  general  resemblance  between 
these  plots  and  corresponding  plots  associated  with  the  combined  Moore 
neighborhood/ increased  bomb  density  scenario  (Figures  46  and  47).  In  both 
scenarios,  most  of  the  data  within  the  time  series  plots  demonstrate  higher  and 
lower  levels  of  Pvery  high{tn)  and  Pvery  iow{t„),  respectively,  than  were  observed  within 
the  other  jammer  bomb  scenario  variants.  Also,  tire  return  maps  associated  with 
both  the  current  and  combined  Moore  neighborhood/increased  bomb  density 
scenarios  demonstrate  the  metamorphosis  of  the  diagonal  cluster  (and  branching 
subclusters)  of  Pvery  ugh{tt)  data  points  observed  in  other  scenario  return  maps  into 
smaller  clusters  symmetric  with  respect  to  the  slope  =  1  diagonal.  However,  the 
crisply  defined  periodic  behavior  that  emerged  in  the  early  portion  of  the  time 
series  plot  from  the  previous  scenario  has  been  replaced  by  two  data  clusters 
extending  from  f„  =  8  min  out  to  about  t„  =  68  (i.e.,  the  point  of  maximmn 
average  jammer  bomb  density  within  tire  brigade).  This  change  is  also  reflected 
in  the  collective  structure  depicted  in  Figure  55,  where  the  three  enlongated 
clusters  associated  with  the  Pvery  high{tn)  time  series  from  the  previous  scenario 
(Figure  47)  have  metamorphosized  into  five  compressed,  ovoid-like  clusters, 
which  are  again  symmetric  relative  to  the  slope  =  1  diagonal.  The  new  data- 
smeared  structure  that  emerges  within  the  susceptibility  gradient  scenario  plots 
is  likely  resultant  from  the  previously  mentioned  variable  susceptibility 
threshold  levels  assigned  across  the  nodes  within  the  brigade. 

3.2.6  Disussion  of  Results 

Results  from  each  of  the  jammer  bomb  scenarios  examined  in  this  section 
indicate  flie  following  set  of  shared  characteristics. 

•  The  dynamic  fraction  of  jammed  node  within  the  brigade,  </is''^"“(tn)>, 
increases  (with  small  non-monotonic  fluctuations)  to  a  maximum  value 
occurring  at  t„  =  68  min,  remains  in  a  state  of  quasi-equihbrium  until  about 
tn  =  95  min,  and  then  steadily  decreases  back  down  to  a  level  of  0. 

•  The  dynamic  SA  connectivity  averaged  across  aU  nodes  within  the  brigade, 
(C(tn)),  decreases  (again  with  small  non-monotonic  fluctuations)  down  to  a 
minimum  value  again  occurring  at  tn  =  68  min,  again  remains  in  a  state  of 
quasi-equihbrium  tmtil  about  tn  =  95  min,  and  then  steadily  increases  back 
up  to  a  level  of  1. 


60 


Figure  54.  Superposition  of  Pvm/  ughitn)  (green  points)  and  Pvm/  law  (tn)  (red  points)  time 
series  data  from  the  susceptibility  gradient  scenario. 
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Figure  55.  Return  map  for  the  Pveryhigh{tn)  (green  points)  and  Pven/iow  (tn)  (red  points)  time 
series  from  the  susceptibility  gradient  scenario. 


•  The  SA  fitness  Pvenj  iom{t„)  time  series  roughly  track  the  scenario-respective 
(h^tress^tn))  time  series  in  terms  of  metric  amplitude. 

•  The  SA  fitness  Pvenj  high{tn)  time  series  roughly  track  the  scenario-respective 
(C{t„))  time  series  in  terms  of  relative  time  series  profile  (with  differences  in 
metric  amplitude). 

•  The  SA  fitness  Phigh{tn)  time  series  profile  is  a  mirror-image  reflection  of  the 
respective  SA  fitness  Pvery  high{tn)  time  series  profile. 

•  The  dynamic  state  entropy  H{tn)  and  temporal  entropy  H(t„,  f„.i) 
demonstrate  separation  between  global  brigade,  CSMA  net,  and  local  net 
node  populations  (although  H{t„)  for  CSMA  and  local  net  populations 
within  a  scenario  often  will  track  one  another  very  closely). 

Thus,  in  all  scenario  variants,  the  dynamic  response  of  the  brigade  SA  network 
tracks  the  invasive  lO  stress,  as  jammer  bombs  diffuse  into  and  then  out  of  the 
brigade  structure.  Also,  this  dynamic  response  is  clearly  partitioned  between 
nodes  which  are  directly  jammed  and  other  nodes  which  suffer  SA  connectivity 
degradation  through  jammed  neighbors. 

However,  the  jammer  bomb  scenario  variants  are  distinguishable  by  their 
distinctive  emergent  behaviors  as  displayed  within  their  respective  Pvery  Wg/i(fn) 
time  series.  Although  variations  in  AO  neighborhoods  and  jammer  bomb 
lifetimes  lead  to  differences  in  collective  structure  as  displayed  within  the  return 
maps,  the  most  noticeable  factor  observed  within  jammer  bomb  simulations 
appears  to  be  the  size  of  the  jammer  bomb  neighborhood  and  its  impact  on  Pvery 
bigh{tn)-  The  original  36-cell  jamming  neighborhood  tends  to  drive  the  network  to 
respond  in  a  roughly  linear  fashion,  with  variable  degrees  of  data  spreading 
during  the  stochastic  second  phase  of  a  simulation  (as  witnessed  in  both  the  Pvenj 
to[f„]  and  Pvery  high[tn\  data).  The  smaller  8-cell  Moore  neighborhood,  on  the  other 
hand,  produces  a  more  abrupt  and  discontinuous  response  in  Pvery  high{tn)  (but 
with  little  change  in  Pvery  iow[tn]).  Finally,  the  susceptibility  gradient  scenario, 
which  effectively  applies  a  spectrum  of  jamming  neighborhood  sizes,  results  in 
network  response  almost  as  equally  abrupt  and  discontinuous  as  seen  with  the 
8-cell  jamming  neighborhoods.  This  suggests  that  the  greater  the  jamming  range 
of  a  stress  source  and  the  more  ubiquitous  this  type  of  source  (i.e.,  the  higher  the 
average  susceptibility  of  network  nodes),  the  more  coordinated  the  overall 
network  response  (i.e.,  the  less  resistant  is  the  network  to  associated 
perturbations  to  SA  coimectivity). 

It  is  interesting  to  speculate  on  the  possible  existence  of  a  specific  jamming 
neighborhood  size  that  results  in  a  phase  transition  relative  to  Pvery  high{t„)  (and 
Phigh[tn]  by  association).  This  type  of  phenomenon,  which  has  been  observed  in 
network  models  of  complex  systems  (Sawhill  and  Kauffman  1996;  Luque  and 
Sole  1997),  is  analogous  to  a  thermodynamic  phase  transition.  In 
thermodynamics,  a  material's  characteristic  properties  can  suddenly  change 
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when  its  temperature  is  raised  or  lowered  by  a  small  amoimt;  the  temperature  at 
which  the  change  in  properties  occurs  is  called  the  critical  temperature  Tc. 
Similarly,  within  the  context  of  an  SA  network  operating  within  a  scenario  with 
uniform  jammer  bombs  (i.e.,  all  bombs  consistently  radiate  at  the  same  power 
level  for  a  fixed  lifetime)  and  susceptibility  gradients,  one  could  define  ttie  critical 
jamming  susceptibility  threshold  relative  to  the  average  network 

susceptibility  threshold  where 

=  (22) 

N  =  total  number  of  network  nodes  within  die  brigade,  and  S,  =  susceptibility 
threshold  of  the  i*  node.  Then,  given  the  existence  of  a  phase  transition  relative 
to  Pvery  high{tn),  two  network  SA  connectivity  phases  might  emerge:  (1)  an 
uncorrelated  decoherent  collective  structure  arising  when  <  Sci'”"'"'"?,  and 

(2)  a  quasi-Unear  coherent  collective  structure  arising  when  ^  Sc 

Since  nodal  susceptibility  thresholds  are  direcdy  correlated  to  a  range  of 
embedded  januning  neighborhoods.  Sc  would  be  associated  witii  an 

average  jamming  neighborhood  size  somewhere  between  4  and  36  cells  (where 
1.0  <  <  3.2;  see  section  3.2.5). 


4.  Conclusions 


A  CA  model  has  been  developed  for  the  analysis  of  the  dynamics  of  situational 
awareness  network  topology  (relative  to  a  receiving  platform)  which  emerges  on 
a  global  scale  within  a  digitized  brigade  exposed  to  lO  stress.  Two  different 
forms  of  stress  were  modeled,  including  (1)  a  mean  field  stress  applied  uniformly 
to  all  brigade  nodes  and  (2)  RF  jammer  bombs  exhibiting  local  fields  within  the 
brigade  structure.  In  particular,  analysis  of  the  fitness  state  data  generated  by  the 
jammer  bomb  lO  scenarios  has  demonstrated  that  local  threat/target  interactions 
affect  inter-platform  dependencies  simultaneously  on  various  scales.  If  a  portion 
of  the  CA  brigade  (such  as  a  battalion)  were  to  be  geographically  separated  from 
the  rest  of  the  brigade,  levels  of  SA  connectivity  throughout  both  of  the  separated 
parts  would  be  impacted  as  a  function  of  node  platform  echelon,  so  that  jamming 
in  one  of  the  separated  parts  would  affect  nodal  situational  awareness  (especially 
in  higher  echelon  platforms)  in  the  other  part  differently  than  if  the  parts  were 
unified.  Thus,  the  emergent  behavior  of  the  stressed  brigade  network  must  be 
analyzed  at  the  holistic  macroscale  level. 

The  spatio-temporal  dynamics  of  other  types  of  dispersed  lO  stress  sources 
which  modify  nodal  transmit/receive  functionality  are  hkely  to  create  similar 
modes  of  emergent  global  behavior  in  situational  awareness  levels  which  cannot 
be  anticipated  by  studying  isolated  platforms  or  even  local  area  networks 
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(LANs).  Given  that  the  threat/ target  interaction  dynamics  of  an  lO  stress  source 
can  be  discretized  in  space  and  time  at  the  mesoscale  level,  it  is  a  fairly  straight¬ 
forward  process  to  program  the  stress  source  into  the  CA  model  within  the 
context  of  a  scenario.  There  is  also  no  limitation  on  the  number  of  different  types 
of  lO  stress  sources  which  can  be  programmed  into  a  single  scenario  given  that 
the  discrete  spatio-temporal  interaction  dynamics  are  known,  although  the 
complexity  of  multi-threat  scenario  sensitivity  analyses  can  rapidly  increase  as  a 
function  of  inter-threat  synergy. 
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Appendix  A.  Dynamics  of  Situational  Awareness  Message 
Servers 


The  dissemination  of  situational  awareness  (SA)  messages  throughout  a  digitized 
brigade  involves  the  coordinated  interaction  of  several  different  types  of  servers, 
including  (a)  carrier  sense  multiple  access  (CSMA)  servers,  (b)  multisource  group 
(MSG)  servers,  and  (c)  self  servers  (see  section  1.2.1  in  the  main  report  for  more 
details  on  the  operational  roles  of  these  servers).  In  order  to  properly  simulate 
the  process  of  SA  message  dissemination  throughout  the  brigade,  the  dynamics 
of  both  CSMA  and  MSG  servers  must  be  modeled  to  accoimt  for  both  message 
size  and  transmission  channel  baud  rates.  This  includes  botii  messages 
broadcast  upward  out  of  a  particular  net  as  well  as  other  incoming  messages 
which  must  be  broadcast  across  a  local  net  once  received  by  a  server.  In  all  cases, 
the  disseminated  SA  messages  will  represent  K5.1  variable  message  format 
(VMF)  messages,  which  report  the  positions  of  friendly  platform  nodes  within 
the  brigade  network.  Finally,  an  SA  message  is  limited  to  a  single  bit  in  order  to 
simulate  server  queue  loading  only  (as  opposed  to  measuring  message  network 
transit  times). 

Figure  A-1  illustrates  the  dynamics  of  outbound  SA  message  processing  witiiin  a 
CSMA  server.  First,  the  server  receives  a  message  over  its  local  single  channel 
groimd  airborne  radio  system  (SINCGARS)  net  (this  message  is  also  received  by 
all  other  members  of  the  net).  Next,  the  server  determines  a  course  of  action 
based  on  the  total  number  of  messages  currently  stored  within  its  queue  (this 
number  can  include  a  self-generated  SA  message  in  addition  to  received 
messages). 

•  If  only  one  message  is  queued,  the  server  will  wait  4  s,  transmit  the 
message  over  its  CSMA  net  via  broadcast,  and  then  reset  an  embedded 
timer  which  measures  message  waiting  time  within  the  queue. 

•  If  two  messages  are  queued,  the  server  will  combine  the  messages  into  a 
CSMA  net  communication  transmission  unit  (CTU)  containing  two  bits  and 
transmit  the  CTU  immediately,  and  then  reset  its  message  timer. 

•  If  more  than  two  messages  are  queued,  the  server  will  combine  two 
messages  into  a  CSMA  net  CTU  and  transmit  the  CTU  immediately. 

(1)  If  there  is  now  one  message  left  in  the  queue,  the  server  will  wait  four 
seconds,  transmit  the  message,  and  then  reset  its  message  timer. 

(2)  If  there  are  two  or  more  messages  left  in  the  queue,  the  server  will  wait  one 
second,  and  then  recheck  the  number  of  queued  messages  (thus  restarting 
the  whole  process). 
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Figure  A-1.  Dynamics  of  outbound  SA  message  processing  within  a  CSMA  server. 

This  process  is  restarted  every  time  the  CSMA  server  receives  one  or  more  new 
SA  messages. 

Figure  A-2  illustrates  the  dynamics  of  SA  message  processing  within  an  MSG 
server  (which  only  sends  messages  outbotmd  across  the  brigade-wide  MSG  net). 
First,  the  server  receives  a  message  over  its  CSMA  net  (this  message  is  also 
received  by  all  other  enhanced  position  location  reporting  system  (EPLRS)- 
equipped  members  of  the  net).  As  with  the  CSMA  server  described  previously 
the  MSG  server  determines  a  course  of  action  based  on  the  total  number  of 
currently  queued  messages. 

•  If  one  to  four  messages  are  queued,  the  server  will  wait  4  s,  transmit  the 
messages  over  the  MSG  net  via  broadcast,  and  then  reset  its  message  timer. 

•  If  five  messages  are  queued,  the  server  will  combine  the  messages  into  an 
MSG  net  CTU  containing  five  bits,  transmit  the  CTU  immediately,  and  then 
reset  its  message  timer. 
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•  If  more  than  five  messages  are  queued,  the  server  will  combine  five 
messages  into  a  MSG  net  CTU  and  transmit  the  CTU  immediately. 

(1)  If  there  are  now  one  to  four  messages  left  in  the  queue,  the  server  will  wait 
4  s,  transmit  the  messages,  and  then  reset  its  message  timer. 

(2)  If  there  are  five  or  more  messages  left  in  the  queue,  the  server  wiU  wait  1  s, 
and  then  recheck  the  number  of  queued  messages  (thus  restarting  the 
whole  process). 

As  with  the  CSMA  server,  this  process  is  restarted  every  time  the  MSG  server 
receives  one  or  more  new  SA  messages. 

There  is  a  simplifying  assumption  made  within  this  MSG  server  model. 
Basically,  it  is  assumed  that  an  MSG  net  CTU  can  be  transmitted  within  a  1-s 
time  interval,  which  is  the  dtuation  of  a  simulation  time  step  using  the  mean 
field  model  (see  section  3.1).  In  reality,  given  that  (1)  such  a  CTU  is  600  bits  in 
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length  and  (2)  the  transmission  baud  rate  over  the  MSG  net  is  480  bits  per  second 
(bps),  only  0.8  of  a  CTU  can  actually  be  transmitted  during  one  simulation  time 
step.  Thus,  the  server  queue  saturation  effects  induced  by  GPS  spoofing  will 
actually  evolve  at  a  faster  rate  than  predicted  by  the  mean  field  model; 
i.e.,.staggered  reduction  in  average  SA  connectivity  levels;  see  Figure  12  in 
section  3.1  of  the  main  report). 

Lastly,  Figure  A-3  illustrates  the  dynamics  of  incoming  SA  message  processing 
within  a  CSMA  server.  In  this  instance,  messages  are  received  over  the  server's 
CSMA  net  and  the  brigade-wide  MSG  net,  and  are  subsequently  rebroadcast 
over  the  server's  local  SINCGARS  net.  First,  the  server  receives  a  message 
(which  is  also  received  by  aU  other  EPLRS-equipped  members  of  the  net  over 
which  the  message  was  broadcast).  Again,  the  server  determines  a  course  of 
action  based  on  the  total  number  of  currently  queued  messages. 


Figure  A-3.  D)mamics  of  incoming  SA  message  processing  within  a  CSMA  server. 

•  If  1-50  messages  are  queued,  the  server  wiU  immediately  transmit  all  of  the 
messages  over  the  local  SINCGARS  net  via  broadcast. 

•  If  more  than  50  messages  are  queued,  the  server  will  transmit  50  of  the 
messages,  wait  1  s,  and  then  recheck  the  number  of  queued  messages  (thus 
restarting  the  whole  process). 
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The  50  message/ time  step  limit  is  due  to  (1)  the  data  transmission  baud  rate 
associated  with  a  SINCGARS  net  is  assumed  to  be  16,000  bps  and  (2)  the  size  of 
an  SA  message  broadcast  over  a  SINCGARS  net  is  320  bits.  Thus,  16,000  bps  / 
320  bits  per  message  =  50  messages  transmitted  per  second.  Of  course,  since 
SINCGARS  net  voice/  data  C2  messages  are  not  considered  in  this  model,  the 
actual  SA  message  transmission  rate  would  likely  be  considerably  lower  than 
this  estimate.  Finally,  this  process  is  restarted  every  time  the  CSMA  server 
receives  one  or  more  new  incoming  SA  messages. 
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Intentionally  left  blank. 


Appendix  B.  Cellang  Source  Code  for  the  Situational 
Awareness  Network  Model  With  Jammer 
Bomb  Stress  Sources 


This  appendix  appears  in  its  original  form,  without  editorial  change. 
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Intentionally  left  blank. 
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sa_net_jain  — 

Situational  awareness  dissemination  network  with  jammer  bombs. 


This  CA  models  a  tactical  digital  communication  network  at  the  multiple 
BRIGADE  level,  where  the  network  architecture  is  based  on  a  generic 
Tactical  Internet  (TI)  implementation.  The  model  thus  allows  for  data 
communication  within  a  single  tactical  BRIGADE.  The  current  version  of 
this  code  models  a  spatial/geographical  network  architecture  where  all 
nodes  move  forward  at  the  same  fixed  speed,  so  that  network  spatial 
configuration  is  essentially  static  (actually,  nodes  —  mounted  on  vehicles 
—  do  not  move  forward,  but  rather  invasive  enemy  "jammer”  bombs  move 
downward;  see  below) .  Cells  represent  different  types  of  TI  nodes  within 
the  architecture  (i.e.  clients  and  servers);  these  nodes  represent 
networked  Army  systems  at  either  the  BRIGADE,  BATTALION,  COMPANY,  or 
PLATOON  military  echelon  levels. 


The  model  runs  as  a  cyclic  process  with  *FOUR*  phases  per  simulation 

cycle  (i.e.  1  four-phase  cycle  =  1  realtime  increment): 

Phase  1:  New  or  "advancing"  jammer  bombs  move  into  empty  cells. 

Phase  2:  Perturb  the  state  of  a  TI  node  cell  from  functional  to 
dysfunctional  if  a  jammer  bomb  lies  within  a  symmetric 
jamming  neighborhood  relative  to  the  cell. 

Phase  3:  Determine  the  operational  status  of  the  cell  * s  local  and  CSMA 
servers  (this  status  will  be  self-referential  if  the  reference 
cell  is  itself  a  local  or  CSMA  server) .  The  status  report 
includes  only  the  current  functionality  of  both  servers .  This 
information  will  be  used  in  Phase  4  in  order  to  determine  the 
reference  cell’s  tactical  internet  connectivity. 


Phase  4:  Evaluate  the  operational  fitness  (fit/degraded/unfit )  of  all  TI 
node  cells  based  on  (i)  self-function  and  (ii)  neighbor  cell 
function  within  broadcast  neighborhood  (i.e.  tactical  internet 
connectivity) .  Determine  whether  jammer  bombs  advance  into 
adjacent  empty  cells,  cease  to  operate  (due  to  battery  drainage), 
or  are  destroyed. 


dimensions  of 

#  Output  either 

# 

#  (a)  the  SA  connectivity  (0  to  1000)  which  is  then  divided  by  1000 

#  within  the  data  filter  program,  or 

# 

#  (b)  the  SA  communication  fitness  as  either  VERY  HIGH  (5) , 

#  HIGH  (4),  MEDIUM  (3),  LOW  (2),  or  VERY  LOW  (1). 

# 

output  of  0..1000  #  NOTE:  currently  set  up  to  output  SA  connectivity 

#output  of  0..5 
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#  The  xpos  and  ypos  fields  record  the  absolute  x  and  y  positions, 

#  respectively,  of  a  cell  relative  to  the  2D  CA  lattice. 

# 

xpos  of  0 . . 128 
ypos  of  0 . . 128 

#  0  =  Empty  cell 

#  1  =  TI  node  system  with  SINCGARS  only  (no  EPLRS) 

#  2  =  TI  node  system  with  SINCGARS  AND  EPLRS  —  capable  of 

#  functioning  as  local  stub  net  server 

#  3  =  Self  server  (i.e.  SINCGARS/EPLRS-eguipt  system  with 

#  no  other  systems  in  local  stub  net;  must  act  as  its 

#  own  local  stub  net  server) . 

#  4  =  CSMA  Primary  SA  server  (pre-planned) 

#  5  =  CSMA  Alternate  SA  server  (pre-planned) 

#  6  =  Jammer  bomb 

#  7  =  Empty  cell  which  absorbs  entering  jammer  bombs. 

# 

type  of  0 . . 7 

#  The  echelon  field  identifies  the  military  echelon  of  a  TI  node. 

#  Possible  echelons  are 

# 

#  4  — >  BRIGADE-level  (e.g.  Brigade  Commander) 

#  3  — >  BATTALION-level  (e.g.  Battalion  TOC  S3) 

#  2  — >  COMPANY-level  (e.g.  Company  First  Sargent) 

#  1  — >  PLATOON-level  (e.g.  Platoon  Wing  Man) 

#  0  — >  not  a  TI  node  (i.e.  an  empty  cell) 

# 

echelon  of  0 . . 4 

#  fitness  =  5  (VERY  HIGH) ; 

#  fitness  =  4  (HIGH) ; 

#  fitness  =  3  (MEDIUM) ; 

#  fitness  =  2  (LOW) ; 

#  fitness  =  1  (VERY  LOW) . 

# 

fitness  of  0..10 

#  The  function  field  indicates  the  current  RECEIVE  functionality 

#  of  the  system  hardware/software  required  for  TI  connectivity, 

#  where 

#  1  =  system  is  currently  able  to  transmit  AND  receive 

#  SA  messages; 

#  0  =  system  is  currently  able  to  transmit  BUT  unable  to 

#  receive  SA  messages. 

#  Thus,  only  message  reception  is  adversely  affected  by  a  jammer 

#  bomb  in  the  vicinity  of  a  cell. 

# 

function  of  0..1 

#  The  server_local  field  indicates  a  server  that  is  dynamically 

#  selected  to  function  as  a  local  stub  net  server.  Any  SINCGARS/ 

#  EPLRS-equipt  system  may  act  as  a  local  stub  net  server  (i.e. 

#  cell  types  2,  4,  or  5  —  type  3  is  already  a  server  by 

#  definition  and  doesn't  require  this  field) .  The  local  stub  net 

#  server  is  a  TI  node  that  maintains  and  disseminates  SA  messages 


#  to  and  from  other  client  nodes  within  a  local  SA  routing 

#  neighborhood,  and  also  broadcasts  SA  data  generated  within  its 

#  local  stub  net  outward  across  the  associated  CSMA  net. 

# 

#  0  =  non-server 

#  1  =  server. 

# 

server_local  of  0..1 

#  The  server_csma  field  indicates  a  pre-determined  active  CSMA 

#  net  server;  this  function  may  be  assumed  by  either  the  primary  or 

#  alternate  CSMA  servers  {cell  types  4  and  5,  respectively) .  The 

#  CSMA  server  receives  SA  data  from  its  local  CSMA  net  and  then 

#  re-broadcasts  this  data  across  the  BRIGADE-wide  MSG  net. 

# 

#  .0  =  non-server 

#  1  =  server. 

# 

#  Since  CSMA  servers  across  different  BATTALIONS  (including  the 

#  BRIGADE-AREA  unit;  see  below)  must  exchange  "broadcast  shares" 

#  between  themselves  for  the  privilege  of  broadcasting,  the  token 

#  field  is  used  to  identify  the  CSMA  server  which  currently  holds 

#  the  authority  to  broadcast. 

# 

server_csma  of  0..1 

#  The  bn  field  identifies  the  BATTALION  or  CSMA  community  to 

#  which  a  TI  node  belongs.  A  BRIGADE-AREA  unit,  which  maintains 

#  its  own  CSMA  SA  net,  is  undistinguishable  from  a  subordinate 

#  BATTALION  within  the  context  of  this  model. 

# 

bn  of  0. . 10 

#  The  stub  field  identifies  the  local  STUB  NET  to  which  a  TI  node 

#  belongs.  A  local  stub  net  is  itself  contained  within  an  BATTALION, 

#  which  is  in  turn  contained  within  a  BRIGADE.  Thus,  the  5th  STUB  NET 

#  within  the  3rd  BATTALION  is  identified  by  the  cell  fields 

# 

#  cell. stub  =  5 

#  cell. bn  =  3. 

# 

stub  of  0. . 50 

#  The  target  field  designates  an  empty  cell  (type  =  0)  which  can  either 

#  (a)  become  occupied  by  a  previously-positioned  jammer  bomb  (with 

#  probability  PI)  or  (b)  serve  as  a  target  drop  point  for  a  jammer  bomb 

#  dropped  by  an  over-flying  aircraft  (with  probability  P2) .  In  these 

#  cases,  cell. target  =  1  or  cell. target  =  2,  respectively. 

# 

target  of  0 . . 2 

#  The  j ammer___timer  field  keeps  track  of  the  lifetime  of  a  jammer  bomb. 

#  This  field  is  passed  as  an  agent  to  represent  a  bomb  moving  downward 

#  relative  to  the  BRIGADE  structure. 

# 

jammer_timer  of  0..300 
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#  The  connect  field  represents  a  fraction  which  measures,  for  a 

#  reference  cell,  the  number  of  currently  functional  cells  which  can 

#  potentially  provide  INCOMING  SA  data  to  the  reference  cell  divided 

#  by  the  total  number  of  such  cells  (including  both  functional  and 

#  dysfunctional  cells) .  Since  this  code  is  limited  to  integer  opera- 

#  tions,  the  above  fraction  is  represented  as  the  integer 

#  connect  =  (fraction  *  1000)modulo  0  (i.e.  rounded  down  to  next 

#  integer),  where  0  <=  connect  <=  1000. 

# 

connect  of  0..1000 

#  The  server_status  field  encodes  the  current  operational  status  of  a 

#  cell’s  local  stub  net  server  and/or  CSMA  server  (if  the  cell  itself 

#  is  a  local/CSMA  server,  then  the  report  status  is  self-referential) . 

#  This  field  is  set/reset  for  each  TI  node  cell  during  phase  2  of  a 

#  simulation  cycle  and  then  queried  during  phase  3  of  the  same  cycle 

#  in  order  to  calculate  a  cell’s  connectivity  measure.  The  array 

#  elements  of  the  server_status  field  are 

# 

#  cell.server_status [0] :  current  local  server  function  (0/1); 

# 

#  cell. server_status [1] :  current  CSMA  server  function  (0/1). 

# 

server_status []  for  2  of  0..1 
agent  of 

#  The  bomb  and  timer  fields  indicates  a  message  sent  from  a  bomb- 

#  occupied  cell  to  an  adjacent  empty  cell  into  which  the  bomb  will 

#  "move.”  The  bomb  agent  field  indicates  the  presence  of  a  jammer 

#  bomb;  the  timer  agent  field  contains  the  same  information  as 

#  cell .  j  ammer_timer . 

# 

bomb  of  0 .  .  1 
timer  of  0 . . 300 

#  The  transfer_local  field  indicates  a  message  sent  from  a 

#  dysfunctional  local  stub  net  server  to  a  dynamically-selected 

#  alternate  server  within  the  local  stub  net.  A  field  value  of  "1" 

#  instructs  the  alternate  server  to  assume  the  role  of  active  stub  net 

#  server. 

# 

transfer_local  of  0..1 

#  The  transfer_csma  field  indicates  a  message  sent  from  a  dysfunctional 

#  primary  CSMA  server  to  a  pre-determined  alternate/standby  server 

#  within  the  same  BAT TALI ON /CSMA  community.  A  field  value 

#  of  "1”  instructs  the  alternate  server  to  assume  the  role  of  active 

#  CSMA  server.  A  subsequent  message  with  a  field  value  of  ”0" 

#  (sent  from  the  primary  to  alternate  server  upon  functional  recovery 

#  of  the  primary)  instructs  the  alternate  server  to  relinquish  its 

#  server  role. 

# 

transfer  csma  of  0..1 


end 
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#  Set  the  maximum  neighborhood  search  ranges  to  cover  the  entire  cell 

#  universe  based  on  the  pre-set  constants  x_range  and  y_range.  This  takes 

#  advantage  of  the  toroidal  curvature  of  the  cell  universe. 

# 


const 

imin  : 

=  -47 

# 

1  -  (x  range  of  automata/2) 

const 

imax  : 

=  48 

# 

X  range 

of  automata  /  2 

const 

jmin  : 

-  -63 

# 

1  -  (y  range  of  automata/2) 

const 

jmax  : 

=  64 

# 

y  range 

of  automata  /  2 

#  Set 

#  the 

up  generic  BATTALION-wide 
entire  cell  universe. 

neighborhood  search 

const 

ilow  : 

=  -24 

# 

(imin/2) 

-  1 

const 

ihigh 

:=  25 

# 

(imax/2) 

+  1 

const 

j  low  : 

=  -32 

# 

( jmin/2) 

-  1 

const 

jhigh 

:=  33 

# 

( jmax/2) 

+  1 

#  Operational  fitness  thresholds: 

# 

#  fitness  >  high  — >  VERY  HIGH 

#  fitness  <=  high  &  fitness  >  medium  — >  HIGH 

#  fitness  <=  medium  &  fitness  >  low  — >  MEDIUM 

#  fitness  <=  low  &  fitness  >  very_low  — >  LOW 

#  fitness  <=  very__low  — >  VERY  LOW 

# 

# 

const  high  :=  900 
const  medium  :=  600 
const  low  :=  400 
const  very_low  :=  100 

#  Apply  a  geographic  distance  filter  to  all  TI  node  cells  within  the 

#  BRIGADE. 

# 

const  bde_filter  :=  80 
const  bn__filter  :=  60 
const  co_filter  :=  40 
const  plt_filter  :=  20 

#  Set  the  lifetime  of  a  jammer  bomb  (all  bombs  are  assumed  to  have  the  same 

#  lifetime) . 

# 

const  lifetime  :=  100 

#  Define  a  symmetric  jamming  neighborhood  relative  to  a  center  cell.  By 

#  symmetry  arguements,  the  center  cell  may  either  be  a  jamming  source  (where 

#  neighbors  are  jammed  targets)  or  a  jammed  target  (with  at  least  one 

#  neighbor  as  a  jamming  source) . 

# 

jamming[]  for  36  :=  [0,  1],  [1,  1],  [1,  0],  [1,  -I],  [0,  -1],  [-1,  -1], 

[-1/  0],  [-1,  1],  [0,  2],  [1,  2],  [2,  2],  [2,  1],  [2r  0]  , 
[2,  -1],  [2,  -2],  [1,  -2],  [0,  -2],  [-1,  -2],  [-2,  -2], 
[-2,  -1],  [-2,  0],  [-2,  1],  [-2,  2],  [-1,  2],  [0,  3], 

[1,  3],  [3,  1],  [3,  0],  [3,  -1],  [1,  -3],  [0,  -3], 

[-1,  -3],  [-3,  -1],  [-3,  0],  [-3,  1],  [-1.  3] 

if  time%4  =  0  &  cell. type  =  0  then 


79 


#  Cell  is  a  potential  site  for  a  jammer  bomb  and 

#  time  =  0,  4,  8,  12,  16,  ... 

# 

if  cell. target  =  1  then 

#  A  potential  pre-plant ed  jammer  bomb. 

# 

P__planted  random%1000 
if  P_planted  <  67  then 

#  There  is  a  one-in-fifteen  chance  that  a  targeted 

#  cell  will  contain  a  pre-planted  jammer  bomb  with 

#  a  remaining  lifetime  of  70  simulation  cycles  (i.e. 

#  280  timesteps) . 

# 

cell. type  :=  6 
cell . jammer_timer  :=  30 
elsif  P_planted  >=  67  &  P_planted  <  133  then 

#  There  is  a  one-in-fifteen  chance  that  a  targeted 

#  cell  will  contain  a  pre-planted  jammer  bomb  with 

#  a  remaining  lifetime  of  60  simulation  cycles  (i.e. 

#  240  timesteps) . 

# 

cell. type  :=  6 
cell . jammer_timer  :=  40 
elsif  P_planted  >=  133  &  P_planted  <  200  then 

#  There  is  a  one-in-fifteen  chance  that  a  targeted 

#  cell  will  contain  a  pre-planted  jammer  bomb  with 

#  a  remaining  lifetime  of  50  simulation  cycles  (i.e. 

#  200  timesteps) . 

# 

cell. type  :=  6 

cell . jammer_timer  :=  50 

end 

elsif  cell. target  =  2  then 

#  A  potential  air-dropped  jammer  bomb. 

# 

P^dropped  :=  random%1000 
if  P_dropped  <  5  then 

#  There  is  a  one-in-two-hundred  chance  that  a  targeted 

#  cell  will  receive  an  air-dropped  jammer  bomb, 
cell. type  :=  6 

cell . jammer_timer  :=  0 

end 

end 

type  cell. type 

#  Check  for  a  jammer  bomb  moving  into  the  cell  from  a  previously- 

#  occupied  ground  position. 

# 

new_bomb_timer  :=  cell . j ammer_timer 
forall  bomb: agent 

#  Move  a  bomb  into  the  cell  if  appropriate. 

# 

if  bomb. bomb  then 

if  type  =  0  then 

cell. type  :=  6 

cell . j ammer_timer  :=  bomb. timer  +  1 
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elsif  type  =  6  then 

cell . jammer_timer  :=  bomb. timer  +  1 
when  bomb. timer  <  new_bomb_timer 

end 

end 

end 

elsif  time%4  0  &  cell. type  =  6  then 

#  Move  a  jammer  bomb  out  of  its  previous  location. 

# 

cell. type  :=  0 

elsif  time%4  =  1  &  cell. type  >  0  &  cell. type  <  6  then 

#  Cell  is  a  tactical  internet  network  node  and 

#  time  =  1,  5,  9,  13,  17,  ... 

# 

#  Jam  the  cell’s  radios  (cutting  off  tactical  internet  data  reception) 

#  if  there  is  a  jammer  bomb  present  within  a  jamming  neighborhood. 

# 

function  :=  1 
forall  i 

if  j amming [i] . type  =  6  then 
function  :=  0 
exit 

end 

end 

cell . function  :=  function 

elsif  time%4  =  2  &  cell. type  >  0  &  cell. type  <  6  then 

#  Cell  is  a  TI  node  and  time  =  2,  6,  10,  14,  18,  ... 

# 

cell_bn  :=  cell. bn 
cell_stub  :=  cell. stub 
type  :==  cell. type 

#  Establish  the  current  status  of  the  cell’s  local  stub  net  and  CSMA 

#  servers  (this  information  is  self-referential  if  the  cell  is  itself 

#  either  an  active  local  or  CSMA  server) . 

# 

forall  i : ilow. . ihigh 

forall  j : jlow. . jhigh 

neighbor  :=  [i,  j] 

neighbor__f unction  :  =  neighbor .  function 
neighbor_bn  neighbor. bn 
neighbor_stub  :=  neighbor , stub 
server_local  : =  neighbor . server_local 
server_csma  :  =  neighbor .  server__csma 
if  neighbor_bn  =  cell_bn  &  neighbor's tub  = 
cell_stub  &  server_local  then 

#  Neighbor  is  the  cell’s  local  server. 
server_status_0  :==  neighbor_f unction 

elsif  neighbor_bn  =  cell_bn  &  server_csma  then 

#  Neighbor  is  the  cell’s  CSMA  server. 
server_status_l  :=  neighbor_f unction 

end 

end 

end 
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if  cell. type  =  2  &  cell . server_local  & 

(cell . function  =  0  |  cell . function  =  2)  then 

#  Transfer  the  local  stub  net  server  role  to  the 

#  first  available/functional  EPLRS^equipt  system  (i.e. 

#  a  type  2  cell)  found  in  a  search  of  the  stub  net/ 

#  local  SA  community. 

# 

cell .  server__local  0 
if  cell. echelon  =  1  then 

agent  (0,  0,  1,  0)  ->  [0,  1]  when  [0,  1] -type  =  2  & 

[0,  1] .echelon  =  1  & 

[0,  1]. function 

“>  [1,  0]  when  [1,  0] -type  =  2  & 

[1,  0] .echelon  =  1  & 

[1,  0]. function 

->  [0,  -1]  when  [0,  -ll.type  =  2  & 
[0,  -I], echelon  =  1  & 

[0,  -1]. function 
[-1,  0]  when  [-1,  0]  .type  =  2  & 
[-1,  0]. echelon  ==  1  & 
[-1,  0]. function 

“>  cell  otherwise 
elsif  cell. echelon  =  2  then 

agent  (0,  0,  1,  0)  ->  [0,  2]  when  [0,  2] .type  =  2  & 

[0,  2] . echelon  =  2  & 

[0,  2]. function 

->  [2r  0]  when  [2,  0] .type  =  2  & 

[2,  0] .echelon  ~  2  & 

[2,  0]. function 

->  [0,  “2]  when  [0,  -2]. type  =  2  & 
[0,  -2]. echelon  =  2  & 

[0,  ~2]. function 

->  [-2,  0]  when  [-2,  0] .type  =  2  & 
[-2,  0]. echelon  =  2  & 
[-2,  0]. function 

->  cell  otherwise 
elsif  cell. echelon  =  3  then 

agent  (0,  0,  1,  0)  ->  [0,  1]  when  [0,  1]  .type  =  2  & 

[0,  1] .echelon  =  3  & 

[0,  1]. function 

->  [1,  1]  when  [1,  1] .type  =  2  & 

[1,  1] .echelon  =  3  & 

[1,  1]. function 

->  [1,  0]  when  [1,  0] .type  =  2  & 

[1,  0] .echelon  =  3  & 

[1,  0]. function 

->  [1,  --1]  when  [1,  -1]  .type  =  2  & 
[1,  -1] .echelon  =  3  & 

[1,  ~1] . function 

->  [0,  “1]  when  [0,  -1] .type  =  2  & 
[0,  -1] .echelon  =  3  & 

[0,  -1]. function 

“>  [-1,  -1]  when  [-1,  -1] .type  =  2  & 
[-1,  -1]. echelon  =  3  & 
[-*1,  -1]  .  function 

->  [-1,  0]  when  [-1,  0].type  =  2  & 
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elsif 


[-1,  0] .echelon  =  3  & 
[-1,  0]. function 

“>  [-1,  1]  when  [-1,  IJ.type  =  2  & 
[-1,  1] .echelon  =  3  & 
[-1,  1]. function 

“>  cell  otherwise 
elsif  cell. echelon  =  4  then 

agent(0^  0,  1,  0)  ->  [0,  1]  when  [0,  1] .type  =  2  & 

[0,  1] . echelon  =  4  & 

[0,  1] . function 

->  [1,  1]  when  [1,  1] .type  =  2  & 

[1,  1] .echelon  =  4  & 

[1,  1]. function 

->  [1,  0]  when  [1,  0].type  =  2  & 

[1,  0] .echelon  =  4  & 

[1,  0]. function 

->  [1,  -1]  when  [1,  -1] .type  =  2  & 
[1,  -1] .echelon  =  4  & 

[1,  -1] . function 

->  [0,  -1]  when  [0,  -I] .type  =  2  & 
[0,  -I]. echelon  =  4  & 

[0,  -1]. function 

->  [-1,  -1]  when  [-1^  -IJ-type  =  2  & 
[-1,  -1] .echelon  =  4  & 
[-1,  -1]. function 

->  [-1,  0]  when  [-1,  OJ.type  =  2  & 
[-1^  0]. echelon  ==  4  & 
[-1,  0] , function 

->  [-1,  1]  when  [-1,  l].type  =  2  & 
[-1,  1]. echelon  =  4  & 
[-1,  1] . function 

“>  cell  otherwise 

end 

(cell. type  =  4  |  cell. type  =  5)  &  cell . server_csma  & 

(cell . function  =  0  |  cell . function  =  2)  then 

#  Transfer  the  CSMA  server  role  to  the  standby  CSMA 

#  server  if  it  is  currently  functional. 

# 

cell . server_csma  :=  0 
if  cell. echelon  =  2  then 

agent (0,  0,  0,  1)  “>  [2,  0]  when  {[2,  0] .type  =  4  | 

[2f  0] . type  =5)  & 

[2,  0]. echelon  =  2  & 

[2,  0]. function 

->  [-2,  0]  when  ([-“2,  0]  .type  =  4  | 
[“2,  0] .type  =5)  & 

[“2,  0]. echelon  =  2  & 

[-2f  0]. function 

-*>  cell  otherwise 


elsif  cell. echelon  =  3  then 


agent (0,  0,  0,  1)  ->  [1,  0]  when  ([1,  0] .type  =  4  1 

[1,  0] .type  =  5)  & 

[1,  0] . echelon  =  3  & 

[1,  0] . function 

->  [-1,  0]  when  ([-1,  0].type  =  4  | 
[-1,  0] .type  =  5)  & 

[-1,  0] .echelon  =  3  & 
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end 


elsif  cell. echelon  = 
agent (0,  0,  0, 


end 


[~1,  0]. function 
->  cell  otherwise 
4  then 

1)  ->  [1,  0]  when  {[1,  0].type  =  4  | 
[1,  0] .type  =  5)  & 

[1,  0] .echelon  =  4  & 

[1,  0]. function 

”>  [-1,  0]  when  ([-1,  0].type  =  4  | 
[-1,  0] .type  =5)  & 

[-1,  0] .echelon  =  4  & 
[“1,  0]. function 
“>  cell  otherwise 


cell.server^status  [0]  server__status_0 
cell.server_status [1]  :=  server_status_l 

elsif  time%4  =  3  &  cell. type  >  0  &  cell. type  <  6  then 

#  Cell  is  a  TI  node  and  time  =  3,  1,  11,  15,  19,  ... 

# 

cell_xpos  :=  cell.xpos 
cell_ypos  :=  cell.ypos 
cell_function  :=  cell . function 
type  :=  cell, type 
bn  :=  cell. bn 
stub  :=  cell. stub 

server_local  :=  cell . server_local 
server_csma  :=  cell . server_csma 
cell_local_function  :=  cell . server_status  [0] 
cell_csma_function  ;=  cell . server_status [1] 

#  Calculate  a  cell*s  inbound  connectivity  to  the  tactical  internet 

#  (i.e,  examine  the  cell  as  a  RECEIVER  of  SA  messages)  and  the 

#  correlated  SA  fitness. 

# 

echelon  :=  cell. echelon 
number_total  :=  0 
number_f unction  :=  0 

forall  i:imin..imax  forall  j:jmin..jmax 
neighbor  :=  [i,  j] 
xpos  :=  neighbor .xpos 
ypos  :=  neighbor .ypos 

neighbor_local_f unction  :=  neighbor . server_status [0] 
neighbor_csma_f unction  neighbor . server_status [1] 

delta_x_l  :=  cell_xpos  -  xpos 
delta_x_2  :=  xpos  -  cell_xpos 
delta_y_l  :=  cell_ypos  -  ypos 
delta_y_2  :=  ypos  -  cell_ypos 
delta_x  :=  delta_x_l  when  delta_x_l  >  0 
:=  delta__x_2  when  delta_x_2  >  0 
:=  0  otherwise 

delta_y  delta_y_l  when  delta_y_l  >  0 
:=  delta_y_2  when  delta_y_2  >  0 
:=  0  otherwise 
if  echelon  =  1  then 

within  :=  1  when  delta_x  <  plt_filter 
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&  delta_y  <  plt_filter 
:=  0  otherwise 
elsif  echelon  =  2  then 

within  :=  1  when  delta_x  <  co_filter 
&  delta_y  <  co_filter 
0  otherwise 

elsif  echelon  =  3  then 

within  :=  1  when  delta_x  <  bn_filter 
&  delta^y  <  bn__filter 
:=  0  otherwise 
elsif  echelon  =  4  then 

within  :=  1  when  delta_x  <  bde_filter 
&  delta_y  <  bde__filter 
:=  0  otherwise 

else 

within  :=  0 

end 

nuinber_total  number^total  +  1  when  within  &  (i  !=0  |  j  !=0) 

&  neighbor . type  >  0  &  neighbor . type  <  6 

if  within  then  if  type  =  1  then 

#  Cell  is  a  SINCGARS-only-equipt  member  of  a  stub  net. 

# 

if  neighbor . type  =  1  then 

#  Neighbor  only  has  a  SINCGARS  radio  on  board  (no  \ 

#  EPLRS) . 

# 

if  neighbor. bn  =  bn  &  neighbor . stub  =  stub  then 

#  Neighbor  is  in  the  same  local  stub  as  the 

#  reference  cell. 

# 

numb er_funct ion  :=  number__function  +  1 
when  (i  !=  0  1  j  !=  0) 

elsif  neighbor. bn  =  bn  &  neighbor. stub  !=  stub  then 

#  Neighbor  is  in  the  same  CSMA  net  as  the  reference 

#  cell,  but  in  a  different  local  stub  net.  Thus,  the 

#  communication  channel  between  the  neighbor  and  cell 

#  also  depends  on  functionality  of  both  the 

#  neighbor’s  and  cell’s  local  servers. 

# 

number^function  :=  number_function  +  1 
when  (i  !=0  |  j  1=0)  & 
cell_local_f unction 
&  neighbor_local_function 

else 

#  Neighbor  and  reference  cell  are  in  different  CSMA 

#  nets.  Thus,  the  communication  channel  between  the 

#  neighbor  and  cell  also  depends  on  functionality  of 

#  (i)  the  neighbor’s  local  server,  (ii)  the 

#  neighbor’s  CSMA  server,  and  (iii)  the  cell’s  local 

#  server. 

# 

number_f unction  :=  number_function  +  1 
when  (i  !=0  |  j  !=0)  & 
cell  local  function  & 
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end 


neighbor_local_f unction  & 
neighbor_csma_f unction 


elsif  neighbor .type  =  2  |  neighbor . type  =  3  | 
neighbor . type  =4  |  neighbor . type  =  5  then 

#  Neighbor  is  equipt  with  an  EPLRS  data  radio. 

# 

if  neighbor. bn  =  bn  &  neighbor . stub  =  stub  & 
neighbor . type  !=  3  then 

#  Neighbor  is  in  the  same  local  stub  as  the 

#  reference  cell. 

# 

number_f unction  :=  nuinber_function  +  1 
when  (i  !=  0  |  j  !=  0) 

elsif  neighbor. bn  =  bn  &  neighbor . stub  !=  stub  then 

#  Neighbor  is  in  the  same  CSMA  net  as  the  reference 

#  cell,  but  in  a  different  local  stub  net. 

# 

if  (neighbor .type  =  2  &  ! neighbor . server_local)  | 
neighbor . type  =  4  |  neighbor . type  =  5  then 

#  The  communication  channel  between  the 

#  neighbor  and  cell  depends  on 

#  functionality  of  the  cell’s  local 

#  server  and  the  neighbor’s  local  server. 

# 

number_function  :=  number___function  +  1 
when  (i  !=0  !  j  !=0)  & 
cell_local_f unction 
&  neighbor_local_function 

elsif  (neighbor,  type  ==  2  &  neighbor .  server_local)  | 
neighbor . type  =  3  then 

#  The  communication  channel  between  the  neighbor  and 

#  cell  depends  on  functionality  of  the  cell’s  local 

#  server. 

# 

nuinber_f unction  :=  number_function  +  1 

when  (i  !=0  |  j  !=0)  & 
cell_local_f unction 

end 


else 

#  Neighbor  and  reference  cell  are  in  different  CSMA  nets. 

# 

if  (neighbor .type  =  2  &  ! neighbor . server_local)  | 

neighbor .type  =  4  |  neighbor . type  =  5  then 

#  The  communication  channel  between  the  neighbor  and 

#  cell  depends  on  functionality  of  (i)  the  cell’s 

#  local  server,  (ii)  the  neighbor’s  local  server,  and 

#  (iii)  the  neighbor’s  CSMA  server. 

# 

number_f unction  :=  number^function  +  1 

when  (i  !=0  |  j  !=0)  & 
cell_local_f unction 
&  neighbor_local_function  & 
neighbor_csma_f unction 
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elsif  (neighbor .type  =  2  &  neighbor . server_local)  | 
neighbor .type  =  3  then 

#  The  communication  channel  between  the  neighbor  and 

#  cell  depends  on  functionality  of  the  cell’s  local 

#  server  and  neighbor’s  CSMA  server. 

# 

number_f unction  :=  number^function  +  1 

when  (i  !=0  |  j  !=0)  & 
cel l_local_f unction 
&  neighbor_csma_function 

end 

end 

end 

elsif  type  =  2  |  type  =  3  |  type  =  4  |  type  =  5  then 

#  Cell  is  any  EPLRS-equipt  TI  node. 

# 

if  neighbor . type  =  1  then 

#  Neighbor  only  has  a  SINCGARS  radio  on  board  (no  EPLRS) . 

# 

if  neighbor. bn  =  bn  &  neighbor . stub  =  stub  then 

#  Neighbor  is  in  the  same  local  stub  as  the  reference 

#  cell. 

# 

number_function  :=  number_f unction  +  1 
when  (i  !=  0  |  j  !=  0) 


elsif  neighbor. bn  =  bn  &  neighbor . stub  !=  stub  then 

#  Neighbor  is  in  the  same  CSMA  net  as  the  reference 

#  cell,  but  in  a  different  local  stub  net.  Thus,  the 

#  communication  channel  between  the  neighbor  and  cell 

#  also  depends  on  functionality  of  the  neighbor’s 

#  local  server. 

# 


number_function  :=  number_f unction  +  1 
when  (i  1=0  |  j  !=0)  & 
neighbor_local_f unction 

else 

#  Neighbor  and  reference  cell  are  in  different  CSMA 

#  nets.  Thus,  the  communication  channel  between  the 

#  neighbor  and  cell  also  depends  on  functionality  of 

#  (i)  the  neighbor’s  local  server  and  (ii)  the 

#  neighbor’s  CSMA  server. 

# 

number_function  :=  number_function  +  1 
when  (i  !=0  |  j  !=0)  & 
neighbor_local_f unction  & 
neighbor_csma_f unction 

end 


elsif  neighbor .  type  =  2  I  neighbor .  type  =  3  |  neighbor .  type 
=  4  I  neighbor .type  =  5  then 

#  Neighbor  is  equipt  with  an  EPLRS  data  radio. 

# 

if  neighbor. bn  =  bn  &  neighbor . stub  =  stub  then 

#  Neighbor  is  in  the  same  local  stub  as  the  reference 

#  cell. 


87 


# 


nuinber_function  :=  nuinber_function  +  1 
when  (i  !=  0  |  j  1=  0) 

elsif  neighbor. bn  =  bn  &  neighbor . stub  !=  stub  then 

#  Neighbor  is  in  the  same  CSMA  net  as  the  reference 

#  cell,  but  in  a  different  local  stub  net. 

# 

if  (neighbor . type  =  2  &  ! neighbor . server_local)  | 
neighbor .type  =  4  |  neighbor . type  =  5  then 

#  The  communication  channel  between  the  neighbor  and 

#  cell  depends  on  functionality  of  the  neighbor *s 

#  local  server. 

# 

number_f unction  :=  number_f unction  +  1 
when  (i  !=0  |  j  !==0)  & 
neighbor_local_f unction 

elsif  (neighbor. type  =  2  &  neighbor . server_local)  | 
neighbor . type  =  3  then 

#  The  communication  channel  between  the  neighbor  and 

#  cell  depends  only  on  cell  functionality. 

# 

number^f unction  :=  number_function  +  1 
when  (i  !=  0  |  j  !=  0) 

end 


else 

#  Neighbor  and  reference  cell  are  in  different  CSMA. 

#  nets. 

# 

if  (neighbor . type  =  2  &  Ineighbor . server_local)  I 
neighbor . type  =  4  1  neighbor . type  =  5  then 

#  The  communication  channel  between  the  neighbor  and 

#  cell  depends  on  functionality  of  (i)  the  neighbor’s 

#  local  server  and  (ii)  the  neighbor’s  CSMA  server. 

# 

number^f unction  :=  nuinber_function  +  1 
when  (i  i=0  |  j  !=0)  & 
neighbor_local_function  & 

neighbor_csma_f unction 

elsif  (neighbor .type  =  2  &  neighbor . server_local)  | 
neighbor . type  =  3  then 

#  The  communication  channel  between  the  neighbor  and 

#  cell  depends  on  functionality  of  the  neighbor’s 

#  CSMA  server. 

# 

number_f unction  :=  number^f unction  +  1 
when  (i  !=0  |  j  !=0)  & 
neighbor_csma_f  unction 

end 

end 

end 

end  end 
end  end 
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connect  :=  (number_function*1000)  /  number^total  when 
cel l_f unction 
:=  0  otherwise 


fitness 


5  when  connect 
4  when  connect 
connect 
3  when  connect 
connect 
2  when  connect 
connect 
1  when  connect 
0  otherwise 


>  high 

<=  high  & 

>  medium 

<=  medium  & 

>  low 

<=  low  & 

>  very_low 
<=  very_lcw 


#  Re-assign  server  roles  as  appropriate.  This  is  done  at  the  end  of 

#  the  phase  3  timestep  so  as  not  to  effect  current  cell  message 

#  reception  fitness  metrics.  Thus,  the  actual  role  transfer  is 

#  effectively  delayed  until  the  next  simulation  cycle  (i.e.  the  next 

#  time  step) . 

# 

local_transfer  :=  0 
csma  transfer  :=  0 


forall  message : agent 

#  Acknowledge  receipt  of  server  role  transfer  authorization 

#  messages.  The  actual  role  transfer  is  delayed  until  the 

#  next  simulation  cycle  (i.e.  the  next  timestep). 
local_transfer  :=  1  when  type  =  2  &  message . transfer__local  -  1 
csma_transfer  1  when  (type  =  4  |  type  =  5)  & 

message. transfer_csma  =  1 


end 


#  Assign  server  roles  to  designated  cells  based  on  reception  of 

#  transfer  authorization  messages.  The  receiver  of  the  server 

#  authorization  message  SHOULD  be  currently  functional. 

# 

server__local  :=  1  when  type  =  2  &  local_transf er 
server__csma  :=  1  when  (type  =  4  |  type  =  5)  &  csma_transfer 

cell . server_local  :=  server_local 
cell . server_csma  :=  server_csma 
cell. connect  :=  connect 
cell. fitness  :=  fitness 

cell, output  :=  connect 

#  cell. output  fitness  when  cell. type  >  0  &  cell. type  <  6 

#  :=  0  otherwise 

elsif  time%4  =  3  &  cell. type  =  6  then 

#  Cell  is  a  jammer  bomb  and  time  =  3,  7,  11,  15,  19,  ... 

# 

jammer_timer  :=  cell .  jammer_timer 
new_position  :=  [0,  -1] 

if  jammer^timer  <  lifetime  &  new^position. type  =  0  then 

#  Send  a  message  to  move  the  jammer  bomb  downward 

#  one  lattice  cell. 
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# 

agent  (1,  jammer_timer,  0,  0)  -“>  [0,  -1] 

end 

end 
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Appendix  C.  Data  Filter  Source  Code 


This  appendix  appears  in  its  original  form,  widiout  editorial  change. 
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INTENTIONALLY  LEFT  BLANK. 
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/*  f liter J am. c 

This  data  filtering  program  should  ONLY  be  used  with  the  Cellang  code 
modeling  the  various  jammer  bomb  10  scenarios. 

*/ 


#include  <stdio.h> 

^**^*^**^**************vr*****  Model  Constants  ★**************^**************/ 

/*  Is  the  number  of  time  steps  in  each  complete  phase  of  the  computation.  */ 
#define  TIME_PHASES  4 

/*  The  total  maximum  number  of  tactical  internet  nodes  in  the  simulation.  */ 
#define  MAX__N0DES  612 

/*  Determines  the  consecutive  number  of  complete  time  phase  computations  to 
keep. 

#define  TIME_PERIOD  1 

/*  The  total  number  of  tactical  internet  echelon  levels  in  the  simulation.  */ 
#define  NUM_0F_ECHEL0NS  4 

/*  The  total  number  of  possible  nodal  fitness  states  in  the  simulation.  */ 
#define  NUM_OF_NODE_STATES  5 

/*  Various  constants  for  dealing  with  different  sampling  windows /partitions . 

#define  NUM__OF_BN_NETS  7 

♦define  NUM_OF_STUB_NETS  26 

/*  This  is  the  maximum  of  NUM_OF_BN_NETS  and  NUM_OF_STUB_NETS .  */ 

♦define  MAX_NUM_OF_SUBNETS  26 

y*^*vr*******-,t****************  Data  Structures  ***★*★★★*★****★***★***********/ 

typedef  struct  { 

int  vhigh_nodes [TIME_PERIOD] ; 
int  high_nodes [TIME_PERIOD] ; 
int  med_nodes [TIME_PERIOD] ; 
int  low_nodes [TIME_PERIOD] ; 
int  vlow_nodes [TIME_PERIOD] ; 

}  echelon; 

static  echelon  fitness_count [NUM_0F_ECHEL0NS] ; 

/*  This  structure  maintains  the  total  amount  of  time  that  each  node  has  thus 
far  spent  in  a  particular  state. 

typedef  struct  { 

int  vhigh_time,  high_time,  med_time,  low_time,  vlow_time; 
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}  node__fitness_tiines; 

static  node__f itness__times  node_times  [MAX_NODES]  ; 

/'^  This  structure  maintains  the  current  fitness  state  of  each  node  and  the 
consecutive  time  to  present  that  the  node  has  been  in  that  state.  E.g. 

time_in__state  [3]  .  state  ==  2  and  time_in_state  [3]  .  time  ==  5 

means  that  the  4th  node  has  been  in  state  2  for  the  past  5  time  periods. 


*/ 

typedef  struct  { 
int  state; 
int  time; 

}  node_t ime_in_s t at e ; 

static  node_time_in_state  time_in_state  [MAX_NODES]  ; 

/*  This  structure  maintains  the  state  that  each  node  was  in  during  the  most 
recent  TIME_PERIOD  previous  time  steps /phases .  The  zeroth  (0th)  entry 
is  the  most  recent  (i.e.  current)  state. 

*/ 

static  int  last_states  [MAX___NODES]  [TIME_PERIOD]  ; 

/*  This  structure  indicates  the  position  within  the  CA  universe  of  each 

node  as  well  as  its  echelon  level,  the  battalion  to  which  it  belongs,  and 
the  stub  net  it  is  a  part  of. 

*/ 

typedef  struct  { 
int  X,  y; 
int  echelon; 
int  bn_net; 
int  stub_net; 

}  inf o_for_node; 

static  info_for_node  node_info [MAX_NODES] ; 

/*  Counts  the  number  of  nodes  belonging  to  each  echelon  level.  Where 

echelon_count [3]  =  #  Brigade  nodes 
2  =  #  Battalion  nodes 

1  =  #  Company  nodes 

0  =  #  Platoon  nodes 

static  int  echelon_count  [NUM_OF_ECHELONS]  ; 

/*^*******^**************^**:^vr  Data  Gathering  **************************^***/ 
extern  int  dim_size(int  dim)  ; 

extern  long  int  get__field  (int  field,  int  dim_0_index,  int  dim_l_index)  ; 
extern  unsigned  long  int  get_time  (void) ; 

static  unsigned  long  int  time; 

static  int  x_dim_si2e,  y__dim_size; 
static  int  number  of  nodes  =  0; 
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ttdefine  OUTPUT_FIELD  0 

ttdefine  TYPE_FIELD  3 

#define  ECHELON_FIELD  4 

#define  BN_FIELD  11 

#define  STUB_FIELD  12 

/*  Input  array  files  that  the  CA  takes  the  x  and  y  coordinates  for  the 
nodes  from. 

*/ 

#define  X_NODE_COORD  "nodes_x” 

#define  Y_NODE_COORD  "nodes_y" 

/*  Initializes  the  data  structures  in  preparation  for  getting  the  nodal 
data  at  the  end  of  each  time  phase.  0  is  returned  if  operations  are 
okay  and  should  continue,  1  is  returned  if  there  is  a  problem. 

*/ 

static  int  init_data (void)  { 
int  i,  j; 

FILE  *x_coord_f ile,  *y_coord_f ile; 

/*  Initialization  at  time  ===  0  */ 
x_dim_size  =  dim_size(0); 
y_dim_size  =  dim_size(l); 

memset (echelon_count ,  0,  NUM_OF_ECHELONS*sizeof (int) ) ; 
memset (last_states,  0,  MAX_NODES*TIME_PERIOD*sizeof (int) ) ; 

x_coord_file  =  fopen (X_NODE_COORD,  "r”) ; 
y_coord_file  =  fopen ( Y_NODE_COORD,  "r") ; 

/*  Make  sure  the  files  got  opened.  */ 
if  {x_coord_file  ==  NULL)  { 

fprintf (stderr,  "FILTER  ERROR:  Unable  to  open  the  file  ’%s*\n", 
X_NODE_COORD)  ; 
return  1; 

} 

if  {y_coord_file  ==  NULL)  { 

fprintf (stderr,  "FILTER  ERROR:  Unable  to  open  the  file  *%s*\n”, 
Y_NODE_COORD) ; 
return  1; 

} 

/*  Get  position  and  echelon  values  for  each  node.  */ 
while  (1)  { 

int  X,  y,  echelon; 
int  x_scan,  y_scan; 

x_scan  =  fscanf (x_coord_file,  "  %d",  &x) ; 
y_scan  =  fscanf (y_coord_file,  "  %d",  &y) ; 

if  (x_scan  <=  0  |  |  y__scan  <=  0)  break; 
else  if  (x_scan  +  y__scan  ==  1)  { 

/*  A  different  number  of  coordinates  was 
placed  in  each  file. 

*/ 

fprintf (stderr,  "FILTER  ERROR:  Number  of  value" 
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I  ir 


in  coordinate  files 


X  NODE  COORD 


and  Y_NODE_COORD  " *  are  not  the  ” 

^*same .  \n”)  ; 
return  1 ; 

}  else  if  (x  <  0  ||  x  >  dim_size(0))  { 

fprintf (stderr,  "FILTER  ERROR:  The  value  of  the  x” 
"  coordinate  from  X_NODE_COORD  at” 

"  position  %d  is  out  of  bounds. \n", 
number_of_nodes+l ) ; 
return  1; 

}  else  if  (y  <  0  M  y  >  dim_size(l))  { 

fprintf (stderr,  "FILTER  ERROR:  The  value  of  the  y" 
"  coordinate  from  Y_NODE_COORD  ”*  at" 

"  position  %d  is  out  of  bounds. \n”, 
number_of_nodes+l )  ; 
return  1; 

} 


if  {number_of_nodes  ==  MAX_NODES)  { 

fprintf (stderr,  "FILTER  ERROR:  Increase  the  size  of” 

"  MAX_NODES  in  filter. c  and  recompile . \n" )  ; 
return  1; 

} 

node_inf o [number_of_nodes] . X  =  x; 

node_info [number_of_nodes] . y  =  y; 

echelon  =  get_field  (ECHELON^FIELD,  x,  y) ; 

/*  Check  for  a  legal  echelon  value.  */ 

if  (echelon  <  1  I  I  echelon  >  NUM_OF_ECHELONS )  { 

fprintf (stderr,  "FILTER  ERROR:  Echelon  value  (%d)  for 
”  node  at  (%d,  %d)  is  out  of  range. \n", 
echelon,  x,  y) ; 
return  1; 

} 

node_info  [number_of__nodes]  .  echelon  =  echelon-1; 

echelon_count [echelon-1] ++/ 

node_times  [number_of_nodes]  .  vhigh_time 

=  node_times  [number_of_nodes]  . high_time 
=  node_times [number_of_nodes] .med_time 
=  node_times [number_of_nodes] .low_time 
=  node_times [number_of_nodes] . vlow_time  =  0; 

node__info  [nuinber_of_nodes]  .bn_net  =  get_field  (BN_FIELD,  x, 

node_inf  o  [number_of_nodes]  .stub_net  = 

get__field  (STUB^FIELD,  x,  y)  ; 

/*  A  node  cell  should  NEVER  be  in  state  0,  thus 
whatever  state  it  is  in  will  not  match  and  the 
data  gathering  loop  below  will  set  the  correct 
initial  value. 

*/ 

time  in  state [number  of  nodes]. state  =  0; 


t ime_in_s tat e  [numb er__of_nodes]  .  time  =  0; 
number_of__nodes++  ; 

} 

fclose  (x_coord_file)  ; 
fclose (y_coord_file) ; 

/*  Zero  out  all  elements  of  the  echelon  time  phase  counts.  */ 
for  (i  =  0;  i  <  NUM_OF__ECHELONS;  i++) 

for  (j  =  0;  j  <  TIME_PERIOD;  j++)  { 

fitness_count [i] . vhigh_nodes [ j ]  =  0; 
fitness_count [i] .high_nodes [ j ]  =  0; 
fitness__count  [i]  .med_nodes  [  j  ]  -  0; 
fitness__count  [i]  .  low__nodes  [  j  ]  =  0; 
fitness_count  [i]  .  vlow__nodes  [  j  ]  =  0; 

} 

return  0; 

} 

/*  Gets  the  data  at  the  end  of  each  complete  time  phase.  0  is  returned  if 
operations  are  okay  and  should  continue,  1  is  returned  if  there  is  a 
problem. 

*/ 

static  int  get_data (void)  { 
int  i ,  j  ; 

/*  Fitness  counts  for  each  echelon  during  this  time  step.  */ 
int  vhigh_count [NUM_OF_ECHELONS] , 
high_count [NUM_OF_ECHELONS] , 
med_count [NUM_OF_ECHELONS]  , 
low_count  [NUM_OF__ECHELONS]  , 
vlow_count [NUM_OF_ECHELONS] ; 

/*  Initialize  echelon  fitness  counts.  */ 
memset  ( vhigh_count ,  0,  NUM__OF_ECHELONS*sizeof  (int)  )  ; 
memset (high_count,  0,  NUM_OF_ECHELONS*sizeof (int) ) ; 
memset {med_count,  0,  NUM_OF_ECHELONS*sizeof (int) ) ; 
memset (low_count,  0,  NUM_OF_ECHELONS*sizeof (int) ) ; 
memset (vlow_count,  0,  NUM_OF_ECHELONS*sizeof (int) ) ; 

/*  Get  the  data  for  each  node.  */ 
for  (i  =  0;  i  <  number_of_nodes;  i++)  { 
int  X  =  node_info [i] .x; 
int  y  =  node_info [i] . y; 
int  echelon  ==  node_info  [i]  .  echelon; 

int  fitness  =  (int)  get_field  {OUTPUT___FIELD,  x,  y)  ; 

/*  This  works  ONLY  when  the  cell. output  field  (i.e., 
field  0)  in  the  Cellang  code  mirrors  the  cell. fitness 
field. 

*/ 

switch  (fitness)  { 

case  5:  /*  VERY  HIGH  */  ^ 

node_times [i] . vhigh_time++; 
vhigh_count [echelon] ++; 

break; 
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} 


case  4:  /*  HIGH  */ 

node_times [i] .high_time++; 
high_count [echelon] ++; 

break; 

case  3:  /*  MEDIUM  */ 

node_tiines  [i]  .ined_time++; 
ined_count  [echelon]  ++; 

break; 

case  2:  /*  LOW  '^ / 

node^times [i] . low_time++; 
low_count [echelon] ++; 

break; 

case  1:  /*  VERY  LOW  */ 

node_times [i] . vlow_time++; 
vlow_count [echelon] ++; 

break; 


/*  Shift  the  remembered  states  and  add  the  latest  one.  */ 
for  (j  =  TIME__PERI0D“1;  j  >  0;  j  — ) 

last_states  [i]  [ j ]  =  last^states [i]  [ j-1] ; 
last_states [i] [0]  =  time_in_state [i] . state; 

if  {time_in_state [i] . state  ==  fitness) 
time_in_state [i] .time++; 

else  { 

time_in_state  [i]  .  state  ==  fitness; 
time_in_state  [i] . time  =  1; 

} 

} 

/*  Shift  the  fitness  counts  for  the  previous  time  steps  and  then 
include  the  count  for  the  current  time  step. 

*/ 

for  (j  =  0;  j  <  NUM_OF_ECHELONS;  j++)  { 

for  (i  =  TIME_PERIOD-l;  i  >  0;  i  — )  { 

fitness_count  [  j  ]  .  vhigh_nodes  [i]  = 

f itness_count [ j ] . vhigh_nodes [i-1] ; 
fitness_count [ j ] . high_nodes [i]  = 

fitness_count [ j ] . high_nodes [i-1] ; 
fitness_count [ j ] .med_nodes [i]  = 

fitness_count [ j  ]  .med_nodes [i-1] ; 
fitness_count [ j ] . low_nodes [i]  = 

f itness_count [ j ] . low_nodes [i-1] ; 
fitness_count [ j ] . vlow_nodes [i]  = 

fitness_count [ j ] . vlow_nodes [i-1] ; 

} 


fitness_count [j] .vhigh_nodes [0]  =  vhigh_count [j]; 
fitness_count [ j ] .high_nodes [0]  =  high_count [ j ] ; 
fitness_count [ j ] .med_nodes [0]  =  med_count [ j ] ; 
fitness_count [ j ] . low_nodes [0]  =  low_count[j]; 
fitness_count [ j ] . vlow_nodes [0]  =  vlow_count [ j ] ; 


return  0; 

} 
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y*************^**^^******  Filter  Configuration  File  *vr*****^****:;r*  +  ***vr*^:fc-**/ 

/*  This  is  the  file  which  contains  the  specification  for  what  type  of  filter 
to  do  and  the  kind  of  graph  specifications  to  generate  for  feeding  to 
gnuplot . 

*/ 

#define  FILTER_SPEC_FILE  "f ilter . conf ” 

/*  By  default,  the  first  filter  (i.e.  ”filter_l”)  will  be  performed.  */ 
static  int  desired^f ilter  =  1; 

/★*★***★*★***★★**★*★*★★★★****  Specific  Filters  ****★★***★*★*★★**★*★*★*******/ 
#include  <math.h> 


/*  The  filter_level  determines  which  kinds  of  nodes  will  be  considered  in 
building  certain  statistics. 

Level  Meaning 


0  All  nodes 

1  Only  brigade  nodes 

2  Only  battalion  nodes 

3  Only  platoon  and  company  nodes 

4  CSMA  net  communities  are  averaged  together 

5  Local  stub  net  communities  are  averaged  together 


*/ 

static  int  filter_level  =  0; 

/*  Basic  zeroth  filter.  Outputs  the  fraction  of  jammed  nodes  within  a 
simulation  cycle.  NOTE:  This  filter  will  work  correctly  ONLY  when 
the  cell. output  field  in  the  Cellang  code  mirrors  the  complemented 
value  of  the  cell . function  field.  In  other  words, 

cell. output  :=  1  -  cell . function . 


*/ 

static  void  filter_0  (void)  { 

int  i,  total  =  0,  node_count  =  0; 

/*  Add  up  the  #  of  nodes  for  the  current  time  step /phase. 

*/ 

for  (i  =  0;  i  <  numb er_of_n odes;  i++)  { 

if  (filter_level  ==  1  &&  node_info [i] . echelon  !=  3) 
continue; 

else  if  ( filter_level  ==  2  &&  node^info [i] . echelon  !=  2) 
continue; 

else  if  (filter_level  ==  3  &&  node_info [i] . echelon  !=  0 

&&  node_info [i] . echelon  !=  1) 

continue; 

total  +=  time_in_state [i] . state; 
node_count++  ; 

} 
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} 


printf(»»%ci  %lf\n",  time/TIME_PHASES, 

(double)  total  /  (double)  node_count) ; 


/*  Basic  first  filter.  Outputs  the  situational  awareness  connectivity 
within  a  simulation  cycle.  NOTE:  This  filter  will  work  correctly  ONLY 
when  the  cell. output  field  in  the  Cellang  code  mirrors  the  value  of  the 
cell. connect  field. 

*/ 

static  void  filter_l  (void)  { 

int  i,  total  =  0,  node_count  =  0; 

/*  Add  up  the  #  of  nodes  for  the  current  time  step/phase. 

*/ 

for  (i  =  0;  i  <  number_of__nodes;  i+t)  { 

if  ( filter_level  ==  1  &&  node_info [i] . echelon  !=  3) 
continue; 

else  if  ( f ilter_level  ==  2  &&  node_info [i] . echelon  !=  2) 
continue; 

else  if  (filter_level  ”  3  &&  node_info [i] . echelon  !=  0 

&&  node_info [i] . echelon  !=  1) 

continue; 


total  +=  time_in_state  [i]  .  stated- 
node  count ++; 


} 


printf(”%d  %lf\n”,  time/TIME_PHASESd 

(double)  total  /  (double)  node_count  /  1000.0); 


/*  Basic  second  filter.  Outputs  total  and  fraction  of  cells  in  each  of  five 
fitness  states  for  the  current  time  step. 

*/ 

static  void  filter  2  (void)  { 


int  i,  jd  temp; 

int  step[NUM__OF_NODE_STATES]  ,  node_count  =  0; 


/*  Initialize  arrays.  */ 

memset (step,  0,  NUM_OF_NODE_STATES*sizeof (int) ) ; 


/*  Add  up  the  #  of  nodes  in  each  fitness  category  for  each  echelon 
for  the  current  time  step/phase. 

*/ 

for  (i  =  0;  i  <  NUM_OF_ECHELONS;  i+4-)  { 

if  (filter_level  ==  1  &&  i  !=  3) 
continue; 

else  if  (filter_level  ==  2  &&  i  !=  2) 
continue; 

else  if  (filter__level  ==  3  &&  i  !=  0  &&  i  !=  1) 
continue; 


step[4] 
step [3] 
step  [2] 
step  [1] 


+=  fitness_count [i] . vhigh_nodes [0] ; 
+=  f itness_count [i] . high_nodes [0] ; 
+=  f itness_count  [i]  .med__nodes  [0]  ; 

+=  fitness  count [i].low_nodes[0]; 
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} 


step[0]  +=  fitness  count [ij . vlow  nodes [Oj; 


} 

node_count  =  step[0]  +  step[l]  + 

printf("  %d",  time/TIME_PHASES; 
printf(”  %lf",  (double)  step[4] 
printf{”  %lf",  (double)  step[3] 
printf("  Ilf",  (double)  step[2] 
printf("  Ilf",  (double)  step[l] 
printf("  llf\n",  (double)  step[ 


step[2]  +  step[3]  +  step[4]; 


(double)  node_count); 

( double )  node_count ) ; 

( double )  node_count ) ; 
(double)  node_count) ; 

/  (double)  node  count); 


/*  Basic  third  filter.  This  produces  both  the  state  and  temporal  entropies. 
The  helper  function  "calc_entropies"  makes  it  easier  to  do  the  averaging 
over  both  bn  and  stub  networks.  Sub_net  is  the  network  number  to 
calculate  the  entropy  for  if  the  filter^level  is  either  4  or  5. 


static  calc_entropies (double  *state,  double  ^temporal,  int  sub_net); 

static  void  filter_3  (void)  { 

double  statement ropy,  temporal_entropy; 
int  i  ; 


if  (f ilter_level  ==  4)  { 

/*  Average  the  entropies  over  all  bn  nets.  */ 

double  total_state  =  0.0,  total_temporal  =  0.0; 
for  (i  =  0;  i  <  NUM__OF_BN_NETS;  i++)  { 

calc_entropies (&state_entropy,  &temporal_entropy,  i) ; 
total__state  +=  state_entropy; 
total_temporal  +=  temporal_entropy; 

} 

state_entropy  =  totalis tate  /  NU]yi_OF_BN_NETS; 
temporal_entropy  =  total_temporal  /  NU]yi_OF_BN_NETS ; 

}  else  if  ( filter_level  ==  5)  { 

/*  Average  the  entropies  over  all  stub  nets.  */ 


} 


double  total_state  =  0.0,  total__temporal  =  0.0; 
for  (i  =  0;  i  <  NUM_OF_STUB_NETS;  i++)  { 

calc_entropies ( &state_entropy ,  &temporal_entropy ,  i ) ; 
total_state  +=  state_entropy; 
total_temporal  +=  temporal_entropy; 

} 

state_entropy  =  total_state  /  NUM_OF_STUB_NETS; 
temporal_entropy  =  total_temporal  /  NUM_OF_STUB_NETS; 

}  else 

calc_entropies (&state_entropy,  &temporal_entropy,  0) ; 

/*  Print  the  results.  */ 

printf("ld  Ilf  llf\n",  time/TIME_PHASES, 

-’State_entropy,  -temporal_entropy)  ; 


static  double  last_P  [NUM_OF_NODE_STATES]  [MAX_NUM_OF_SUBNETS]  ; 
/*  Helper  function  (and  the  real  guts)  for  filter_3,  */ 
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static  calc_entropies  (double  ’^state,  double  ^temporal,  int  sub_net)  { 

double  state_entropy  =  0,  P [NUM_OF_NODE_STATES] ; 

double  teinporal_entropy  =  0, 

P_cond[NUM_OF_NODE_STATES] [NUM_OF_NODE_S TATES] ; 

int  j,  k,  step[NUM_OF_NODE_STATES]  ,  node__count  =  0; 

/*  Initialize  arrays.  */ 

memset (step,  0,  NUM_OF_NODE_STATES*sizeof (int ) ) ; 
memset  (P,  0,  NUM__OF_NODE_STATES*sizeof  (double)  )  ; 
memset (P_cond,  0, 

NUM  OF  NODE  STATES^NUM  OF  NODE  STATES^sizeof (double)); 


/*  Add  up  the  #  of  nodes  in  each  fitness  category  for  each  echelon, 
but  only  for  the  current  time  step/phase. 

*/ 

for  (k  =  0;  k  <  number_of__nodes ;  k++)  { 

if  (filter_level  ==  1  &&  k  1=  3) 
continue; 

else  if  (f ilter_level  ==  2  &&  k  1=  2) 
continue; 

else  if  (filter_level  ==  3  &&  k  i=  0  &&  k  !=  1) 
continue; 

else  if  ( f ilter_level  ==  4  &&  node_info [ k] .bn_net  !=  sub_net) 
continue; 

else  if  (f ilter^level  ===  5  &&  node_info  [ k]  .  stub_net  !=  sub_net) 
continue; 


} 


step[0]  += 
step[l]  += 
step[2]  += 
step [3]  += 
step[4]  += 


(time_in_state [k] . state  ==  0); 
(time_in_state [k] . state  ==  1)  ; 
(time_in_state [ k] . state  ==  2); 
(time_in_state [k] . state  ==  3) ; 
(time  in  state [k] . state  ==  4); 


/*  VERY  LOW 
/*  LOW  */ 

/*  MEDIUM  */ 

/*  HIGH  */ 

/*  VERY  HIGH  */ 


node_count  =  step[0]  +  step[l]  +  step [2]  +  step [3]  + 


step [4] ; 


/*  Calculate  the  state  entropy.  */ 

for  (k  =  0;  k  <  NUM_OF_NODE_STATES;  k++)  { 


P[)c]  =  (step[k]  ==  0  ?  0  :  (double)  step  [k]  /  (double)  node_count)  ; 
state_entropy  +=  log (pow (P [k] ,  P[k])); 

} 

state  entropy  =  state_entropy/log  (  (double)  NUM_OF__NODE_STATES)  ; 

/*  Calculate  the  conditional  probabilities.  */ 
for  (j  =  0;  j  <  NUM_OF_NODE_STATES;  j++) 
for  (k  =  0;  k  <  NUM_OF_NODE_S TATES;  k++)  { 

Count  the  number  of  nodes  now  in  state  k  that  were  in 
state  j  during  the  last  time  step/phase. 

*/ 

int  i,  k_count  =  0,  last_j_count  =  0; 
for  (i  =  0;  i  <  number_of_nodes ;  i+t)  { 

if  (filter_level  ==  1  &&  node_info [i] . echelon  !=  3) 
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} 


continue; 

else  if  (filter_level  ==  2  &&  node_info [i] . echelon  !=  2) 
continue; 

else  if  (filter_level  ==  3 

&&  node_info [i] . echelon  !=  0 
&&  node___info  [i]  .  echelon  !=  1) 

continue; 

else  if  (*filter_level  ==  4 

&&  node_info [i] .bn_net  !=  sub_net) 
continue; 

else  if  (filter_level  “  5 

&&  node_info [i] .bn_net  1=  sub_net) 
continue; 


} 


if  (time_in_state[i] .state  ==  NUM_OF_NODE_S TATES -k 
&&  last_states [i] [0]  ==  NUM_OF_NODE_STATES- j ) 
k_count++; 

if  (last_states[i]  [0]  -=  NUM_OF_NODE_STATES’- j  ) 
last_j_count++; 


/*  NOTE:  If  last_P[j]  =  0  then  it  MUST  be  the  case  that 
k_count  =0,  in  which  case  the  conditional  probability 
is  simply  set  to  0  (rather  than  attempting  the  faulty 
division  by  0) . 

*/ 

P__cond[k]  [j]  =  (last_j_count  ==  0  ? 

0  : 

(double)  k__count  /  (double)  last_j_count )  ; 

} 

/*  Calculate  the  temporal  entropy.  */ 
for  (j  =  0;  j  <  NUM_OF_NODE_STATES;  j++) 
for  (k  =  0;  k  <  NUM  OF_NODE_S TATES;  k++)  { 


temporal_entropy  +=  log (pow (P_cond [k] [ j ] , 

last_P[j]  [sub_net]  *  P_cond[k]  [j] ) )  ; 

} 

temporal_entropy  =  temporal_ent ropy/log ( (double)  NUM_OF_NODE_STATES) ; 

/*  Report  the  results.  */ 

Estate  =  state_entropy; 

^temporal  =  temporal_entropy; 

/*  Update  last_P  for  the  next  time  step /phase  calculations.  */ 
for  (k  =  0;  k  <  NUM_OF_NODE_S TATES;  k++) 
last  P [k] [sub_net ]  =  P[k]; 


/*  Basic  fourth  filter.  Produces  return  maps  for  all  five  fitness  level 
time  series.  Each  line  of  output  for  a  given  time  step/phase 
represents  the  fraction  of  nodes  within  each  state  for  the  current  time 
step  as  well  as  for  the  previous  time  step.  This  filter  produces  no 
output  for  T=0. 

*/ 

static  void  filter  4  (void)  { 
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int  i,  step [NUM_OF_NODE_S TATES ] ,  node_count  =  0; 
static  double  last_vhigh_frac  =  0, 
last_high_frac  =  0, 
last_rned_frac  =  0, 
last_low_frac  =  0, 
last_vlow_frac  =  0; 

/*  Initialize  arrays.  */ 

memset (step,  0,  NUM_OF_NODE_STATES*sizeof (int) )  ; 

/*  Add  up  the  #  of  nodes  in  each  fitness  category  for  each  echelon, 
but  only  for  the  current  time  step/phase. 

*/ 

for  (i  =  0;  i  <  NUM_OF_ECHELONS ;  i++)  { 

if  (filter_level  ==  1  &&  i  !=  3) 
continue; 

else  if  (f ilter_level  ==  2  &&  i  !=  2) 
continue ; 

else  if  (filter_level  ==  3  &&  i  !=  0  &&  i  !=  1) 
continue; 

step[4]  +=  f itness_count [i] . vhigh_nodes [0]  ; 
step[3]  +=  fitness_count [i] . high_nodes [0]  ; 
step[2]  +=  fitness_count [i] .med_nodes [0] ; 
step[l]  +=  fitness_count [i] . low_nodes [0]  ; 
step[0]  +=  fitness_count [i] . vlow_nodes [0] ; 

} 

node_count  =  step[0]  +  step[l]  +  step[2]  +  step[3]  +  step[4]; 

/*  No  output  for  T=0 .  */ 
if  (time  >0)  { 

printf ("%d",  time/TIME_PHASES) ; 

printfC  %.2f",  (double)  step[4]  /  (double)  node_count)  ; 
printf ("  %f",  last_vhigh_frac) ; 

printfC  %.2f",  (double)  step[3]  /  (double)  node_count)  ; 
printfC  %f",  last_high_frac)  ; 
printfC  %.2f",  (double)  step[2]  /  (double)  node_count); 

printfC  %f",  last_med_frac)  ; 
printfC'  %.2f",  (double)  step[l]  /  (double)  node_count); 

printfC  %f",  last_low_frac)  ; 
printfC  %.2f",  (double)  step[0]  /  (double)  node_count); 
printfC  %f\n",  last_vlow_frac)  ; 

} 

last_vhigh_frac  =  (double)  step [4]  /  (double)  node_count; 

last  high_frac  =  (double)  step [3]  /  (double)  node_count; 
last~med_frac  =  (double)  step [2]  /  (double)  node_count; 
last_low_frac  =  (double)  step[l]  /  (double)  node_count; 
last  vlow  frac  =  (double)  step[0]  /  (double)  node_count; 

/*  This  is  the  main  driver  for  filtering.  */ 
int  do_filter  (void)  { 

time  =  get_time(); 


if  (time  ==  0)  { 

FILE  *config; 


} 


} 


if  (init^data  ( )  )  return  1;  Quit  now  if  there  is  a  problem.  */ 

/*  Read  the  desired  filter  from  the  FILTER_SPEC_FILE.  */ 

if  ((config  =  fopen  (FILTER_SPEC_FILE,  "r") )  =-  NULL)  { 
fprintf (stderr, 

"FILTER  ERROR:  Cannot  open  file  *%s*\n”, 
FILTER_SPEC_FILE)  ; 
return  1; 

} 

if  (1  !=  fscanf (config,  "  %d",  &desired_filter) )  { 
fprintf (stderr, 

"FILTER  ERROR;  No  int  at  beginning  of  file  ’%s’\n", 
FILTER_SPEC_FILE)  ; 
return  1; 

} 

if  (1  !=  fscanf (config,  "  %d",  &f ilter_level ) )  { 

filter_level  =  0; 

} 

close (config) ; 

printf  ("#  FILTER  %d\n#  LEVEL  %d\n", 
desired  filter,  filter  level) ; 


if  (time%TIME_PHASES  ==  0)  { 

if  (get_data())  return  1;  /*  Quit  now  if  there  is  a  problem.  */ 


switch  (desired_filter)  { 

case  0:  filter_0();  break; 
case  1:  filter_l{);  break; 
case  2:  filter_2{);  break; 
case  3:  filter_3();  break; 
case  4:  filter_4();  break; 
default : 

fprintf (stderr,  "FILTER  ERROR;  Unknown  "  \ 
"filter  kind  ’%d*  encountered. \n", 
desired__filter )  ; 
return  1; 

} 

} 


/*  Since  nothing  is  done  before  the  end  of  a  time  phase,  do  nothing.  */ 


return  0; 


105 


Intentionally  left  blank. 


106 


NO.  OF 

COPIES  ORGANIZATION 


NO.  OF 

COPIES  ORGANIZATION 


2  DEFENSE  TECHNICAL 
INFORMATION  CENTER 
DTIC  OCA 

8725  JOHN  J  KINGMAN  RD 
STE0944 

FT  BELVOIR  VA  22060-6218 

1  HQDA 

DAMOFDT 

400  ARMY  PENTAGON 

WASFHNGTON  DC  20310-0460 

1  OSD 

OUSD(A&T)/ODDR&E(R) 

DRRJTREW 

3800  DEFENSE  PENTAGON 
WASHINGTON  DC  20301-3800 

1  COMMANDING  GENERAL 
US  ARMY  MATERIEL  CMD 
AMCRDATF 
5001  EISENHOWER  AVE 
ALEXANDRIA  VA  22333-0001 

1  INST  FOR  ADVNCD  TCHNLGY 
THE  UNIV  OF  TEXAS  AT  AUSHN 
3925  W  BRAKER  LN  STE  400 
AUSHN  TX  78759-5316 

1  US  MILITARY  ACADEMY 

MATH  SCI  CTR  EXCELLENCE 
MADN  MATH 
THAYER  HALL 
WEST  POINT  NY  10996-1786 

1  DIRECTOR 

US  ARMY  RESEARCH  LAB 

AMSRLD 

DR  D  SMITH 

2800  POWDER  MILL  RD 

ADELPHI MD  20783-1197 

1  DIRECTOR 

US  ARMY  RESEARCH  LAB 
AMSRLCIAIR 
2800  POWDER  MILL  RD 
ADELPFD  MD  20783-1197 


3  DIRECTOR 

US  ARMY  RESEARCH  LAB 
AMSRL  Cl  LL 
2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1197 

3  DIRECTOR 

US  ARMY  RESEARCH  LAB 
AMSRL  Cl  1ST 
2800 POWDER  MILLED 
ADELPHI  MD  20783-1197 


ABERDEEN  PROVING  GROUND 

2  DIRUSARL 

AMSRL  Cl  LP  (BLDG  305) 


107 


NO.  OF 

COPIES  ORGANIZATION 


NO.  OF 

COPIES  ORGANIZATION 


1  OASDC3I 

JBUCHHEISTER 
RM  3D174 

6000  DEFENSE  PENTAGON 
WASHINGTON  DC  20301-6000 

1  OUSD  AT  STRT  TAC  SYS 

DRSCHNEITER 
RM3E130 

3090  DEFENSE  PENTAGON 
WASHINGTON  DC  20301-3090 

1  OUSD(A&T)/S&T  AIR  WARFARE 

RMUTZELBURG 
RM3E139 

3090  DEFENSE  PENTAGON 
WASHINGTON  DC  20301-3090 

1  OUSD(A&T)/S&TLAND 
WARFARE 
AVHLU 
RM3B1060 

3090  DEFENSE  PENTAGON 
WASHINGTON  DC  20310-3090 

1  UNDER  SECY  OF  THE  ARMY 

DUSAOR 
ROOM2E660 

102  ARMY  PENTAGON 
WASHINGTON  DC  20310-0102 

1  ASST  SECY  ARMY 

ACQSTN  LOGISTCS  &  TECHLGY 
SAAL  ZD  ROOM  2E673 

103  ARMY  PENTAGON 
WASHINGTON  DC  20301-0103 

1  ASST  SECY  ARMY 

ACQSTN  LOGISTCS  &  TECHLGY 
SAAL  ZP  ROOM  2E661 
103  ARMY  PENTAGON 
WASHINGTON  DC  20310-0103 

1  ASST  SECY  ARMY 

ACQSTN  LOGISTCS  &  TECHLGY 
SAAL  ZS  ROOM  3E448 
103  ARMY  PENTAGON 
WASHINGTON  DC  20310-0103 


1  OADCSOPS  FORCE  DEV  DIR 
DAMOFDZ 
ROOM  31522 
460  ARMY  PENTAGON 
WASHINGTON  DC  20301-0460 

1  HQDA 
ODCSPER 
DAPEMRRM2C733 
300  ARMY  PENTAGON 
WASHINGTON  DC  20301-0300 

1  US  ARMY  ARMAMENT  RDEC 
AMSTAARTD 
MHSETTEBLDGl 
PICATINNY  ARSENAL  NJ 
07806-5000 

1  US  ARMY  MISSILE  RDEC 
AMSMIRD 
DRWMCCORKLE 
REDSTONE  ARSENAL  AL 
35898-5240 

1  NATICK  SOLDIER  CENTER 

SBCNT 
PBRANDLER 
KANSAS  STREET 
NATICK  MA  01760-5056 

1  US  ARMY  TANK  AUTOMTV  RDEC 

AMSTATR 
J  CHAPIN 

WARREN  MI  48397-5000 

1  US  ARMY  INFO  SYS  ENGRGCMD 

AMSELBETD 
DRFJENIA 

FT  HUACHUCA  AZ  85613-5300 

1  US  ARMY  SMTRNG  INST  CMD 
AMSTICG 
DR  M  MACEDONIA 
12350  RESEARCH  PKWY 
ORLANDO  FL  32826-3726 

1  USARMYTRADOC 

BATTLE  LAB  INTEGRATION 
TECH  &  CONCEPTS  DIR 
ATCDB 

FT  MONROE  VA  23651-5850 


108 


NO.  OF 

COPIES  ORGANIZATION 


NO.  OF 

COPIES  ORGANIZATION 


j 


1  USARMYTRADOCANLCTR 
ATRCW 
AKEINTZ 

WSMR  NM  88002-5502 

1  US  ARMY  RESEARCH  OFFICE 

4300  S  MIAMI  BLVD 
RESEARCH  TRIANGLE  PARK  NC 
27709 

1  US  ARMY  RESEARCH  LAB 

AMSRLSL 
CHOPPER 

WSMR  NM  88002-5513 

1  US  ARMY  RESEARCH  LAB 
AMSRLSL  E 
JPALOMO 

WSMR  NM  88002-5513 

1  US  ARMY  RESEARCH  LAB 
AMSRLSLEA 
R  FLORES 

WSMR  NM  88002-5513 

1  US  ARMY  RESEARCH  LAB 
AMSRLSL  El 
J  NOWAK 

FT  MONMOUTH  NJ  07703-5601 


ABERDEEN  PROVING  GROUND 


ABERDEEN  PROVING  GROUND  (COND 

1  US  ARMY  EVALUATION  CENTER 
CSTEAECSV 
LDELATTRE 

4120  SUSQUEHANNA  AVE 
APGMD  21005-3013 

12  DIRUSARL 
AMSRLSL 
DR  WADE 
JBEILFUSS 
AMSRLSL  B 
J  FRANZ 
W  WINNER 
AMSRLSL  BA 
MRITONDO 
AMSRLSL  BD 
J  MORRISSEY 
AMSRLSL  BE 
DBELY 
AMSRLSL  BG 
A  YOUNG 
AMSRLSL  BN 
D  FARENWALD 
AMSRLSL  E 
M  STARKS 
AMSRLSL  EC 
E  P  ANUSKA 
AMSRLSL  EM 
J  FEENEY 


1  SBCCOMRDEC 
AMSSBRTD 
JZARZYCKJ 
5183  BLACKHAWK  RD 
APG  MD  21010-5424 

1  US  ARMY  DEV  TEST  COM 
CSTEDTCTTT 
APG  MD  21005-5055 

1  US  ARMY  EVALUATION  CENTER 
CSTEAEC 
W  HUGHES 

4120  SUSQUEHANNA  AVE 
APG  MD  21005-3013 


109 


NO- OF 

COPIES  ORGANIZATION 


NO.  OF 

COPIES  ORGANIZATION 


3  INFO  SCIENCES  TEAM 
PHYS  SCIENCE  LAB 
DR  R  BERNSTEIN  JR 
DR  R  SMITH 
DR  M  COOMBS 
PO  BOX  30002 
LAS  CRUCES  NM 
88003-8002 

1  BOOZ  ALLEN  &  HAMILTON  INC 
P  MURPHY 
4001  FAIRFAX  DRIVE 
SUITE  750 

ARLINGTON  V  A  22203 

1  CONVERGENT  COMPUTING  INC 
J  D  ECKART 

1260  RADFORD  ST  SUITE  A 
CHRISTIANSBURG  VA  24073 


ABERDEEN  PROVING  GROUND 

1  US  ARMY  DEV  TEST  CMD 

CSTEDTCTTT 
APGMD  21005-5055 

17  DIRUSARL 
AMSRLSLBE 
PTANENBAUM 
RSANDMEYER 
J  ANDERSON 
R  SAUCIER 
AMSRLSLBN 
EHORAVANTE 
AMSRLSLEA 
D  BAYLOR 
B  RUTH  (10  CPS) 

RZUM  BRUNNEN 


V 


r 


5  US  ARMY  RESEARCH  LAB 
AMSRLSLEA 
DR  K  MORRISON 
DRPDJANG 
DLANDIN 
T  READER 
GGUZIE 

WSMRNM  88002-5513 


1  US  ARMY  RESEARCH  LAB 

AMSRLSLEI 
PBOTHNER 
FT  MONMOUTH  NJ 
07703-5601 


♦ 


no 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
0MB  No.  0704-0188 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  Including  the  time  for  reviewing  Instructions,  searching  existing  data  sources, 
gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  Information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this 
collection  of  Information,  including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson 
Davis  Hlphwav.  Suite  1204,  Artinqiton,  VA  22202-4302.  and  to  the  Office  of  Manaqement  and  Budget.  Paperwork  Reduction  Pro|ect(0704-01M).  WashInqtonJiC  20503, _ 


1 .  AGENCY  USE  ONLY  (Leave  b!ank)  2.  REPORT  DATE  I  3.  REPORT  TYPE  AND  DATES  COVERED 


December  200 1  |  Final,  October  1 999-September  2000 


4.  TITLE  AND  SUBTITLE  5.  FUNDING  NUMBERS 

Modeling  the  Emergent  Dynamics  of  a  Tactical  Situational  Awareness  Network  665604D670 
Within  an  Information  Operations  (10)  Environment 


6.  AUTHOR{S) 

Brian  G.  Ruth  and  J.  Dana  Eckart* 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 
U.S.  Army  Research  Laboratory 
ATTN:  AMSRL-SL-EA 
Aberdeen  Proving  Ground,  MD  21005-5068 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

ARL-TR-2641 


9.  SPONSORING/MONITORING  AGENCY  NAMES(S)  AND  ADDRESS(ES) 


10.SPONSOR1NG/MONITORING 
AGENCY  REPORT  NUMBER 


11,  SUPPLEMENTARY  NOTES 


*Convergent  Conqjuting,  Inc.,  1260  Radford  St,,  Suite  A,  Christiansburg,  VA  24073. 


12a.  DISTRIBUTION/AVAILABILITY  STATEMENT 
Approved  for  public  release;  distribution  is  unlimited. 


13.  ABSTRACTfMax/mum  200  words) 

A  model  developed  to  analyze  the  global  emergent  dynamics  of  situational  awareness  (SA)  dissemination  within  a 
digitized  brigade  undergoing  hostile  information  operations  (10)  stress  events  is  presented.  Networked  battlefield 
platforms  are  modeled  as  interconnected  tactical  internet  nodes  through  the  use  of  cellular  automata.  The  cellular 
automata  form  a  discrete  spatially  extended  dynamical  system  consisting  of  a  parallel  networked  lattice  of  computational 
cells  in  two  dimensions.  The  global  impact  of  localized  10  stress  events  on  SA  dissemination  can  then  be  analyzed  as  a 
function  of  time  by  mapping  the  SA  dissemination  architecture  onto  the  cellular  automata  lattice.  The  lattice  cells  can 
represent  tactical  internet  client  or  server  platforms.  This  model  demonstrates  the  fluctuating  patterns  of  situational 
"self-awareness"  (i.e.,  a  friendly  platform*s  situational  awareness  of  all  other  fiiendly  platforms  within  its  area  of 
operations)  which  emerge  across  the  digitized  brigade  on  a  global  scale  and  also  illustrates  the  emergent  dynamics  of  SA 
dissemination  within  a  "movement-to-contact"  brigade  structure. 


14.  SUBJECT  TERMS 


cellular  automata,  tactical  communication  network,  situational  awareness,  information 
operations,  emergent  behavior 


17.  SECURITY  CLASSIFICATION 
OF  REPORT 

UNCLASSIFIED 


NSN  7540-01-280-5500 


18.  SECURITY  CLASSIFICATION 
OF  THIS  PAGE 

UNCLASSIFIED 


19.  SECURITY  CLASSIFICATION 
OF  ABSTRACT 

UNCLASSIFIED 


15.  NUMBER  OF  PAGES 


16.  PRICE  CODE 


20.  LIMITATION  OF  ABSTRACT 


Standard  Form  298  (Rev.  2-89) 
Prescribed  by  ANSI  Std.  239-18 


Intentionally  left  blank. 


112 


